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SUMMARY 

In  2006,  we  published  a  paper  that  is  based  on  the  present  report1.  Because  of  a  restriction  on  the 
number  of  printed  pages,  we  had  to  omit  a  good  part  of  the  analysis.  The  omissions  were  so 
extensive  that  even  the  author  could  not  follow  the  mathematical  trend  when  looking  at  the  paper 
a  few  years  later.  For  this  reason,  we  went  back  to  the  original  notes  and  compiled  this  report  that 
gives  a  detailed  account  of  the  subject  matter.  In  addition  to  aiding  in  understanding  how  the 
results  in  the  paper  were  obtained,  this  report  can  also  guide  the  reader  in  extending  the  results  to 
other  geometries,  such  as  a  circular  or  rectangular  waveguide. 

The  problem  we  deal  with  here  is  the  classical  problem  of  radiation  of  a  semi-infinite  coaxial  line 
into  a  half-space.  With  respect  to  a  rectangular  coordinate  system  xyz,  the  coaxial  line  is  along 
the  z-axis  and  originates  at  minus  infinity.  The  inner  conductor  terminates  at  z  =  0,  while  the 
outer  one  opens  up  into  an  infinite  plane  perpendicular  to  the  coaxial  line.  Details  on  the  history 
and  bibliography  of  this  problem  are  given  in  the  aforementioned  paper  and  will  not  be  repeated 
here.  Our  work  on  this  problem  differs  from  that  of  other  works  in  that  it  is  an  exact  formulation 
of  the  problem  in  tenns  of  boundary  integral  equations  (BIEs).  The  idea  for  this  approach 
originated  in  our  work  on  scattering  by  an  indentation  on  a  ground  plane2.  Its  main  characteristic 
is  that  the  domain  of  the  integral  equations  does  not  involve  the  infinite  plane. 

In  this  report,  we  not  only  provide  the  missing  analysis  for  the  geometry  we  just  described  but 
we  also  present  the  problem  in  a  more  general  setting.  In  Part  1,  we  derive  BIEs  for  a  monopole 
over  a  perfectly  conducting  plane,  driven  by  a  coaxial  line.  The  monopole  may  or  may  not  be  a 
continuation  of  the  center  conductor  of  the  coaxial  line  and  can  be  quite  general  in  shape.  The 
integral  equations  extend  over  the  monopole,  the  opening  to  free  space  of  the  coaxial  line  and 
over  either  the  semi-infinite  walls  of  the  coaxial  line  or  over  a  finite  part  of  them  and  the  inter- 
wall  spacing  at  the  end  of  the  finite  part. 

In  Part  2,  we  employ  the  equations  we  found  in  Part  1  to  the  problem  we  presented  in  the  paper. 
We  take  advantage  of  the  circular  symmetry  of  the  problem  to  reduce  the  vector  integral 
equations  to  three  scalar  equations.  By  expressing  the  unknown  current  densities  in  terms  of  the 
natural  modes  of  the  coaxial  line,  we  show  that  the  problem  can  be  reduced  to  solving  a  single, 
scalar  integral  equation.  We  also  derive  expressions  for  the  far  fields  in  the  half-space  in  which 
the  coaxial  line  radiates.  In  Part  3,  we  take  advantage  of  the  orthogonality  properties  of  the 
modal  functions  to  convert  the  scalar  integral  equation  into  an  infinite  system  of  linear  algebraic 
equations.  We  also  convert  the  coefficients  of  the  system  from  double  to  single  integrals  and 
proceed  to  show  how  to  compute  them.  In  Part  4,  we  provide  detailed  information  on  how  we 
compute  the  system  of  equations  and  consider  four  different  coaxial  lines  for  which  we  compute 
a  number  of  quantities  of  interest.  In  the  last  part,  Conclusions,  we  summarize  the  work  and  offer 
suggestions  for  further  work. 


1  J.  S.  Asvestas,  “Radiation  of  a  coaxial  line  into  a  half-space”,  IEEE  Trans.  Antennas  Propagat.,  Vol.  54,  No.  6,  pp. 
1624-1631,(2006). 

2  J.  S.  Asvestas  and  R.E.  Kleinman,  “Electromagnetic  Scattering  by  Indented  Screens,”  IEEE  Trans.  Antennas 
Propagat.,  Vol.  42,  pp.  22-30,  (1994). 
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PART  1 

INTEGRAL  EQUATIONS  FOR  A  MONOPOLE  OVER  A  GROUND  PLANE 

ABSTRACT 

In  Part  1,  we  present  a  new  approach  to  the  problem  of  monopole  radiation  over  an  infinite, 
perfectly  conducting  plane.  The  monopole  is  fed  by  an  air-filled,  infinite,  coaxial  line.  The 
approach  we  use  is  mathematically  rigorous  and  physically  exact.  It  leads  to  a  system  of 
boundary  integral  equations  (BIEs)  that  can  be  solved  numerically  using  well  established 
methods. 

Part  1  has  three  chapters.  In  the  first  chapter,  we  derive  BIEs  that  extend  over  the  entire  length  of 
the  walls  of  the  infinite  coaxial  line.  The  monopole  is  not  physically  connected  to  the  center 
conductor  of  the  coaxial  line.  The  BIEs  involve  the  electric  surface-current  density  as  the 
unknown. 

In  the  second  chapter,  the  geometry  remains  the  same  as  in  the  first  part  but  the  BIEs  extend  over 
the  part  of  the  walls  of  the  coaxial  line  that  border  its  connection  to  the  infinite  plane.  This  is 
accomplished  by  introducing  an  additional  unknown,  namely,  the  magnetic  surface-current 
density.  We  show  that  the  results  of  the  first  part  can  be  obtained  from  the  results  of  the  second 
by  a  series  of  transformations  and  we  also  derive  formulas  for  the  far  field. 

In  the  third  chapter,  we  make  the  monopole  a  natural  extension  of  the  center  conductor  of  the 
coaxial  line.  We  then  use  the  results  of  the  first  two  parts  to  obtain  BIEs  for  this,  more  realistic, 
case. 

We  point  out  that  the  shape  of  the  monopole  is  quite  arbitrary  but  that  the  BIEs  are  valid  only 
when  all  materials  are  perfectly  conducting. 
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CHAPTER  1 

EQUATIONS  EXTENDING  OVER  THE  ENTIRE  WALLS  OF  THE  LINE 
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1.  INTRODUCTION 


In  this  chapter,  we  fonnulate  the  problem  of  radiation  by  a  monopole  over  a  ground  plane  in 
tenns  of  BIEs.  The  ground  plane  is  infinite  and  the  monopole  is  fed  by  an  infinite  coaxial  line. 
This  formulation  is  the  foundation  for  the  rest  of  the  report  and  contains  general  comments  that 
will  not  be  found  elsewhere  in  this  report. 

We  derive  BIEs  for  the  geometry  of  Figure  1.1.  All  surfaces  are  perfectly  conducting.  The 
inner  and  outer  radii  of  the  coaxial  line  are  a  and  b,  respectively.  The  line  supports  only  a  TEM 
wave  which  means  that  (references  1  and  2) 


,  a  +  b 

k - 

2 


<1 


(1.1) 


where  k  is  the  wavenumber  of  the  time -harmonic 


( Q+wot )  electromagnetic  wave  in  the  line. 


Z 


t 


(7 

1 

s 

(7 

h 

' 

Figure  1.1.  A  cylindrical  monopole  over  a  ground  plane,  fed  by  a  coaxial  line. 
The  excitation  of  the  line  occurs  at  z  =  -oo  and  results  in  a  TEM  wave  with  fields 
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E's'(r)  = 


Ve 


In 


\a) 


-ikz  YV  p~ikz 

p,  He  (r)  =  7--^  (p  ■ 


P 


In 


Uy 


P 


(1.2) 


where  V  is  the  voltage  of  the  inner  conductor  with  respect  to  the  outer,  and  Y  is  the  free-space 
impedance.  We  have  also  employed  cylindrical  coordinates  ( p,(p,z ).  Since  the  line  is  not 

infinite,  we  also  have  an  induced  wave  in  the  line  with  fields  { E! ,  H'  } .  The  total  fields  { E' ,  H' } 
in  the  coaxial  line  are  the  sum  of  the  generator  and  induced  fields 

E'  =  Eg  +  E'  ,  H'=Hg+H'.  (1.3) 


In  the  upper-half  space,  we  have  radiated  fields  {Er,H' }.  We  proceed  to  determine  these  fields 
by  first  deriving  integral  equations  and  then  solving  them  numerically. 


4 


NAWCADPAX/TR-2013/1 15 


2.  INTEGRAL  REPRESENTATIONS  IN  UPPER  HALF-SPACE 

In  this  section,  we  derive  integral  representations  of  the  electromagnetic  fields  in  the  upper- 
half  space.  From  these  representations,  we  will  eventually  obtain  integral  equations. 

We  apply  Green’s  second  identity  in  the  upper  half-space  of  Figure  1.1,  excluding  the  interior 
of  the  monopole.  Thus,  the  surface  integrals  will  extend  over  the  surface  S  and  the  infinite  plane 
z  =  0 .  There  is  also  a  surface  integral  over  a  hemisphere  with  center  at  the  origin  and  whose 
radius  tends  to  infinity.  We  omit  this  integral  since  the  radiating  fields  satisfy  the  Silver-Muller 
condition  and,  hence,  the  integral  will  tend  to  zero  as  the  radius  of  the  hemisphere  tends  to 
infinity.  From  Green’s  second  identity  (reference  3,  p.  509),  we  have 

J  {£'' (r) •  [ V x V x r, (r, r ')] -[VxVxEr (r)] •  r, (r, r ')} dV 

D+ 

=  |  |[n  x  Er  (r)]  •  V  x  Tj  (r,  r ')  -  V  x  Er  (r)  •  [n  x  fy  (r,  r ')]|  dS  (2.1) 

(z=0)uS 


where  D+  is  the  whole  upper-half  space  except  for  the  region  occupied  by  the  monopole,  and  n 
is  the  interior  unit  nonnal.  The  nonnal  is  equal  to  z  at  z  =  0  and  points  to  the  exterior  of  the 
monopole  on  S.  The  dyadic  T,  has  the  form 

T1(r,r')  =  -ikV  x  [g(r,  r ')  I  +  g(r,  r/)  I,  ]  (2.2) 


where  g  is  the  free-space  scalar  Green’s  function 


g(r,r')  = 


An  r-r' 


(2.3) 


and 


I  =  xx  +  yy  +  zz  ,  I;.  =  xx  +  yy  -  z.z  (2.4) 

are  the  identity  dyadic  and  its  image  about  the  xv-plane.  Moreover,  for  any  vector  a,  we  define  its 
image  a,  about  the  xy-plane  by 

a  =  (ax,ay,a:)oai  =(ax,ay,-az).  (2.5) 

The  dyadic  T  i  satisfies  the  differential  equation 

V xV xT,  (r,r')-A:2ri  (r,r')  =ikV x[<5(r,r')I  +  <5(r,r/)I(.]  (2.6) 


5 


NAWCADPAX/TR-2013/1 15 


and  the  boundary  condition 

zxr!(r,r')  =  0  at  z  =  0  .  (2.7) 

The  electric  field  satisfies  the  reduced  wave  equation 

V  x  V  x  E'  (r)  -  k2Er  (r)  =  0  (2.8) 

and  the  boundary  conditions 

«xE'(r)  =  0,  r  (E  S  ;  zxE'(r)  =  0,  r  e  {(x,y,z)  :  z  =  0,  (x,y)  <£  a)  (2.9) 


Combining  the  last  two  results,  we  find  that 


J [z  x  Er  (r)]  •  V  x  Tl  (r, r ')dS  -  ikz\[n  x  H  (r)]  •  T ,  (r, r ')dS  =  k2Z 

a  S 


fH'  (r')  ,  z'>0 
|h'<  (r/)  ,  z'<0 


(2.12) 


This  is  the  integral  representation  of  the  magnetic  field  in  the  upper-half  space.  If  we  use  the 
definition  of  the  dyadic  from  (2.2),  we  can  write  it  as 


2 J  [z  x  E'-  (r)]  •  [Wg(r ,  r ')  +  k2g(r,r')i,  ]  dS 
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3.  INTEGRAL  REPRESENTATIONS  INSIDE  THE  COAXIAL  LINE 

We  proceed  to  find  a  representation  of  the  magnetic  field  inside  the  coaxial  line.  We  apply 
Green’s  second  identity  to  the  region  D  of  Figure  3.1.  In  the  z-direction,  this  region  extends 
between  z  =  —d  and  z  =  0 .  In  the  lateral  direction,  it  is  bounded  by  the  walls  of  the  coaxial  line. 
Thus,  the  four  surfaces  that  bound  D  are  r,  cr,  Sa  and  Sb .  The  nonnal  is  directed  into  this  region. 

z 


t 


b 

Figure  3.1.  The  region!)  of  application  of  Green’s  second  identity. 


Let 

m  =  cruru  Savj  Sb.  (3.1) 

We  can  use  the  same  identity  as  for  (2.1)  to  obtain  an  integral  representation  for  the  total 
magnetic  field  in  D;  thus, 

J  {E'  (r)  •  [V  x  V  x  T(r,  r')]  -  [v  x  V  x  E'  (r)]  •  T(r,  r ')}  dV 

D 
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=  J  {[«  x  E'  (r)]  •  V  x  T(r,  r ')  +  [«xVxE'  (r)]  •  T(r,  r ')}  dS  (3 .2) 

ZD 

where  T  is  the  free-space  dyadic  Green’s  function,  defined  by 

r(r,r')  =  -ikV  x[g(r,r')l]  (3.3) 

and  satisfying  the  differential  equation 

V  x  V  x  T(r,r')  -  £2r(r,r')  =ikV  x  [<7(r,r')l] .  (3.4) 

The  total  electric  field  in  D  satisfies  the  reduced  wave  equation  (2.8)  and  the  boundary  condition 
«xE'(r)  =  0,  reSavSb.  (3.5) 


As  with  (2.1),  the  volume  integral  I D  in  (3.2)  yields 

pH'(r'),  r'.fl 
1  0,  r'£D 

while  the  surface  integral  Im  becomes 
Im  =  -ikZ  J  j[n  x  H'  (r)]  •  T(r,  r')}  dS 

+  |  |[«  x  E'  (r)]  •  V  x  T(r, r ')  -  ikZ  [«  x  H(  (r)]  •  T(r, r ')}  dS .  (3.7) 

<rur 


Combining  the  last  two  statements,  we  have  that 


J  [fixE'(r)]  -V xf(r,rV5 


-ikz\ {[/}  x  H'(r)]  •  r(r,r')}</£ 

ZD 


j£2ZH'(r'),  r  reD 
{  0,  r’eD’ 


With  the  definition  of  the  dyadic  in  (3.3),  this  becomes 


|[«xE' (r)]  •  [Wg(r,  r ')  +  £2g(r,  r ')l] dS 


-ikzj jj ~h  x  H'(r)J  x  Vg(r,r')j dS  =  |  ^ 


r'sD 

r'g£>‘ 


This  is  the  integral  representation  of  the  magnetic  field  in  D. 


(3.8) 


(3.9) 
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4.  ELIMINATION  OF  THE  ELECTRIC  FIELD 

In  these  section,  we  eliminate  one  of  the  unknowns,  the  electric  field,  from  the  magnetic  field 
integral  representations. 

When  z> 0 ,  we  get  from  (3.9)  that 

J[zxE'(r)]{vVg(r,r')  +  k2g(r,r')l]^  =  J[zxE'(r)]{Wg(r,r')  +  k2g(r,r')l]^ 

(J  T 

-/£Zj{[«xH'(r)]xVg(r,r')W,  z'X).  (4.1) 

m 

Since  the  electric  field  is  continuous  in  crossing  cr,  we  can  substitute  this  in  the  top  line  of  (2.13) 
to  get 

2 J [z  x  E'  (r)]  •  [  VVg(r,  r ')  +  k2g( r, r ')l]  dS  -  z2kZ  J  j[nxH'  (r)]  x  Vg(r,  r')}  dS 

T  ZD 

-ikz\ |[/7.  x  Hr(r)]  x  Vg(r,r')  +  ([«.  x  H'(r)]  x  Vg(r,r/))j  =  ikZW  (r') ,  r ’  e  D+ .  (4.2) 

s  v 

In  a  similar  way,  we  can  substitute  the  bottom  line  of  (2.13)  in  the  top  line  of  (3.9)  to  get 
J  [z  X  E'  (r)]  •  [ VVg(r,  r ')  +  k2g( r,  r')l]  dS 

T 

-,kf\ {[« x Hr(r)] x Vg(r,r')  +  ([« x H  (r)] x Vg(r,r/))  }^ 

-ikzj |[/7  X  H'(r)]  X  Vg(r,r')j  dS  =  ikZ  |h'  (r')  +  ^[Hr  (r/)].  | ,  r'  e  D .  (4.3) 

These  two  expressions  give  us  an  integral  representation  of  the  magnetic  field  everywhere.  The 
integrals  over  crdo  not  involve  the  electric  field. 

We  can  go  a  step  further  and  eliminate  the  electric  field  altogether.  On  the  annular  disk  r,  we 
have  that  z  =  —d  .  We  can  show  that,  as  d  —>  oo ,  the  integral  over  z  tends  to  zero  as  d  .  Taking 
this  limit,  we  can  write 

-  J  {[»  x  H'  (r)]  x  Vg(r ,  r ')  +  ([«x  IT  (r)]  x  Vg(r,  r/)^  }  dS 
sy 

-2j  {[w  x  H'  (r)]  x  Vg(r,  r')}  dS  =  W  (r')  ,  r'  e  D+  (4.4) 

m 

and 
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-  j  J {[/i x  H'' (r)]  x  Vg(r, r ')  +  ([«  x  H  (r)]  x  Vg(r, r/))  }  dS 

-fj[«xH'(r)]x  Vg(r, r ')}  dS  =  H' (r ')  +  i [Hr  (r/)];  ,  r '  e  D .  (4.5) 

m 

In  the  last  two  expressions,  the  region  D  extends  over  the  entire  space  between  the  walls  of  the 
coaxial  line. 

For  convenience,  we  define  electric  current  densities  on  the  various  surfaces 


Ja(r)  =  «xH'(r),re  Sa  (4.6) 

Jt(r)  =  «xH'(r),reSj  (4.7) 

Ja(r)  =  -zxH'(r),reff  (4.8) 

Js(r)  =  nxH'  (r)  ,  r  e  S  (4.9) 


and  substitute  in  (4.4)  and  (4.5) 


These  are  the  two  integral  representations  for  the  magnetic  field  in  terms  of  the  unknown  current 
densities.  Equation  (4.10)  clearly  displays  the  dependence  of  the  radiated  fields  on  the  fields 
inside  the  coaxial  line.  Conversely,  Equation  (4.11)  displays  the  dependence  of  the  total  fields 
inside  the  coaxial  line  on  the  geometry  of  the  monopole  and  the  resulting  currents  on  it.  Both 
expressions  indicate  the  existence  of  a  ground  plane  through  the  presence  of  the  image  of  the 
gradient  of  the  scalar  Green’s  function.  We  will  use  these  integral  representations  next  to  obtain 
integral  equations  for  the  unknown  current  densities. 
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5.  INTEGRAL  EQUATIONS  ON  AND  S 

The  derivation  of  integral  equations  for  the  four  unknown  current  densities  is  based  on  a 
theorem  of  Muller  (reference  4),  Theorem  46,  p.  205).  If  we  have  a  region  D  bounded  by  a 
closed  surface  S,  with  the  unit  normal  n  pointing  in  the  exterior  (away  from  D),  and  if  the 
surface  current  density  J  on  S  is  continuous,  then  at  any  point  r'of  S  (with  normal  n  )  we  have 
that 


n'x  J[j(r)xVg(r,r')]^  =  +-J(r')  +  n'xj[J(r)xVg(r,r')]^  (5.1) 

^  s' 

s ; 

where  the  expression  on  the  left  is  to  be  interpreted  in  the  sense  of  a  limit,  i.e.,  in  the  limit  as  we 
approach  the  surface  from  its  exterior  ( Se )  or  interior  ( Si ). 

If  we  apply  this  theorem  to  (4.10),  with  the  point  r'  on  S  and  the  approach  from  the  exterior, 
we  get 

-2/7' ’  X  J  Ja  (r)  X  Vg(r, r' ’)dS  +  J  3h  (r)  x  Vg(r , r' ')dS  +  J  JCT  (r)  x  Vg(r , r' ')dS 

[s„  sb 

-h'  x  J |j5(r)  x  Vg(r.r')  +  [Js(r)  x  Vg(r,r,')],}rfS'  =  '  J.s  (r')  ,  r'  e  S  .  (5.2) 

Again  using  (4.10),  we  approach  a  point  on  <j 
z  x  J  Jfl  (r)  x  Vg(r ,  r ' ')dS  +  J  Jb  (r)  x  Vg(r,r')dS 

si> 

+  ^ J { Js (r)  X  Vg(r,r')  +  [ J, (r)  x  Vg(r, r/)]  } dS 

+z'xj[JCT(r)xVg(r,r')]rfS'  =  J(T(r'),  r'ecr.  (5.3) 


We  note  that  we  can  also  use  (4. 1 1)  to  obtain  this  last  result.  We  also  observe  that 
Vg(r,r')  =  Vg(r,r/),  z' =  0 
and  that 

^  x  [  K  (r)  x  Vg(r,  r ')]  =  JCT  (r)  -  [f  •  K  (r)]  Vg(r,  r ') . 


(5.4) 


(5.5) 
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The  last  term  here  is  zero  because  the  current  density  does  not  have  a  z-component.  Also,  the 
derivative  of  the  scalar  Green’s  function  with  respect  to  z  is  equal  to  zero  when  z  =  z  =  0 .  Thus, 
in  place  of  (5.3),  we  have 


(5.6) 


We  can  also  use  (4.11)  to  obtain  integral  equations  on  the  walls  of  the  coaxial  line.  If  the 
geometry  is  exactly  as  shown  in  Figure  2,  then  there  is  no  problem.  If,  however,  the  geometry  is 
more  general,  then  we  cannot  use  (4.11).  As  an  example,  consider  the  case  where  the  radius  of 
the  monopole  is  greater  than  b.  Then,  for  a  point  of  Sa  or  Sh  with  z-coordinate  between  —h  and  - 

(h  +  /),  the  point  r/  is  inside  the  monopole,  where  the  radiated  magnetic  field  is  not  defined.  For 
this  reason,  we  develop  a  new  representation  for  the  magnetic  field  inside  the  coaxial  line. 
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6.  A  SECOND  REPRESENTATION  OF  THE  MAGNETIC  FIELD  INSIDE 

THE  COAXIAL  LINE 

For  reasons  presented  at  the  end  of  the  last  section,  we  proceed  to  obtain  another 
representation  of  the  magnetic  field  inside  the  coaxial  line.  In  place  of  (3.2),  we  write 


J  Je'  (r)  •  [V  X  V  X  r2  (r,r')]  -  [v  x  V  x  E'  (r)]  •  f\  (r,r')j  dV 

D 

=  J  j-E'  (r)  •  [«  x  V  x  r2  (r,  r ')]  +  [«  x  V  x  E'  (r)]  •  r2  (r,  r')  j  ^  ,  r  'eh  (6.1) 

m 

where  the  Green’s  dyadic  of  the  second  kind  is  defined  by 

r2  (r,  r ')  =  —ikV  x  [g(r,  r ')  i  -  g(  r,  r/)  i .  ]  (6.2) 

and  satisfies  the  differential  equation 

V  x  V  x  r2  (r,r')-  k2T2  (r,r')  =ikV  x  [£(r,r')I  -  ^(r, r/) I,  ]  (6.3) 

and  the  boundary  condition 

zxVxF(r,r')  =  0  atz  =  0.  (6.4) 

In  place  of  (3.6),  we  then  have 

ID=k2ZH'(r'),  r'efl  (6.5) 

and,  in  place  of  (3.7), 

4  =J{-E'  (r)  •  [/)  x  V  x  r2  (r,  r ')]  +  [n  x  V  x  E'  (r)]  •  T2  (r,  r ')}  dS .  (6.6) 

m 


As  in  Section  3,  the  surface  integral  over  z  vanishes  as  d  — >  oo.  Because  of  (6.4),  the  first  term 
vanishes  over  cr.  It  also  vanishes  over  the  walls  of  the  coaxial  line  since  the  tangential  component 
of  the  total  electric  field  is  zero  there.  We  are  thus  left  with 

4  =  -ikZ  J Ja (r) •  r2 (r , r ')dS  +  J Jb (r)  •  IT (r.r^  +  J Ja (r) •  r2 (r.r^  .  (6.7) 

[S*  Sb  a  J 

Combining  (6.5)  and  (6.7),  we  find 

J  Jfl (r)  •  r2 (r, r')dS  +  J  Jb (r)  •  T2 (r, r')dS  +  j  Ja (r)  •  T2 (r.r'JdS  =  ikW (r') ,  r'  e  D  (6.8) 

Sa  Sb  a 
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and,  if  we  take  definition  (6.2)  of  the  dyadic  into  consideration, 

- 1  { (r)  x  Vg(r,r')  -  [ JB  (r)  x  Vg(r,  r/)]  }  dS  -  J  { J,  (r)  x  Vg(r,  r ')  -  [ J, (r)  x  Vg(r, r/)]. }  dS 

sa  sb 

-J  { K  (r)  x  Vg(r,  r ')  -  [  JCT  (r)  x  Vg(r,  r,')],. }  dS  =  H'  (r') ,  r'  e  D .  (6.9) 

a 

But  on  <j 

Ja(r)x  vS'(r’r)  _[JCT(r)x  Vg(r,r/)],  =2Ja(r)xVg(r,r'),  z  =  0.  (6.10) 

Thus, 

-  J  {• 3  a  (r) x  Vg(r,  r ')  -  [  J8  (r)  x  Vg(r,  r/)]  }  dS  -  J  { Jh  (r)  x  Vg(r,  r ')  -  [  Jb  (r)  x  Vg(r,  r/)]  }  dS 

Sa  S„ 

-2 J  JCT  (r)  x  Vg(r,r’)dS  =  H'(r'),  r' e  D  .  (6.11) 

<7 

This  is  the  second  representation  of  the  magnetic  field  inside  the  coaxial  line. 
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7.  INTEGRAL  EQUATIONS  ON  THE  WALLS  OF  THE  COAXIAL  LINE 

In  this  section,  we  use  the  integral  representation  of  the  magnetic  field  that  we  obtained  in  the 
last  section  to  get  integral  equations  on  the  semi-infinite  walls  of  the  coaxial  line.  To  this  end,  we 
employ  (5.1)  and  (6. 1 1).  In  the  approach  to  the  inner  wall  we  get 

~n  x  |  { Ja  (r)  x  Vg(r,  r')  -  [ja  (r)  x  Vg(r,  r/)]  }  dS 

sa 

~n  x  |  { J6(r)  x  Vg(r,r')  -  [j,(r)  x  Vg(r,r/)].}  dS 

sb 

-2/i'x  f  JCT(r)x  Vg(r,r')dS  =  V(r') ,  feS,  (7.1) 

J  2 

a 

while  in  the  approach  to  the  outer 

~n  x  { { Ja  (r)  x  Vg(r,  r')  -  [  Ja  (r)  x  Vg(r,  r/)]  }  dS 

Sa 

~  n  x  |  ( Jb  (r)  x  Vg(r ,  r ')  -  [  Jb  (r)  x  Vg(r,  r/)] . }  dS 

sb 

-2n'x  f  Jff(r)xVg(r,rV5  =  ^JA(r') ,  r'e^.  (7.2) 

J  2 

a 

Following  this  approach,  we  can  ask  what  happens  to  (6.11)  as  the  observation  point  r' 
approaches  u.  If  we  can  get  an  integral  equation  on  cr,  then,  together  with  (7. 1)  and  (7.2),  we  will 
have  three  equations  in  three  unknowns  and,  thus,  we  will  be  able  to  solve  for  the  three  current 
densities.  We  note,  however,  that  none  of  these  three  equations  contains  information  about  the 
upper-half  space.  This  says  that  the  geometry  of  the  upper-half  space  does  not  influence  the 
behavior  of  these  currents,  which  does  not  seem  to  be  correct.  We  proceed  to  perfonn  this 
calculation.  First,  we  cross  (6. 1 1)  with  the  unit  normal  on  a 

J  z  x  { Ja  (r)  x  Vg(r,  r ')  -  J ,  (r)  x  Vg(r,  r/)}  dS 

Sa 

+  Jzx{j»x  Vg(r,  r')  -  Jb  (r)  x  Vg(r,  r/)}  dS 

s„ 

+  2jz x [JCT(r) x  Vg(r,r')]^  =  JCT (r')  ,  r'eD.  (7.3) 

CJ 

In  the  first  integral  above,  only  the  components  transverse  to  the  z-axis  are  involved;  hence,  we 
can  write  instead 

J  z  x  { Ja  (r)  x  Vg(r,  r')  -  J„  (r)  x  Vg(r,  r/)[  dS 

Sa 
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=  j  {■ Ja  (r)z  •  V  [g(r,  r ')  -  g(r,  r/)]  -  [z  •  Ja  (r)]  V  [g(r,  r ’)  -  g( r,  r/)][  dS  (7.4) 

which  is  equal  to  zero  because  of  (5.4).  The  same  is  true  about  the  second  integral  in  (7.3).  As 
the  observation  point  approaches  the  surface,  we  apply  (5.1)  to  the  remaining  integral  to  get  that 

|zx[j(j(r)xVg(r,r')]t/5,  =  0,  r'  e  <r  (7.5) 

<7 


with  the  last  statement  also  following  from  (5.5).  Thus,  we  have  the  identity  0  =  0  or  that  we 
cannot  obtain  an  integral  equation  from  (6.1 1)  when  the  observation  point  is  on  the  surface  cr. 
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8.  THE  SYSTEM  OF  INTEGRAL  EQUATIONS 

We  have  four  unknowns,  namely,  the  electric  current  densities  defined  in  (4.6)  through  (4.9). 
With  these,  we  associate  the  integral  equations  (5.2),  (5.6),  (7.1),  and  (7.2).  A  drawback  of  this 
system  is  that  the  surfaces  Sa  and  Sb  are  semi-infinite  in  the  z-direction  and,  hence,  they  will 

have  to  be  terminated  in  any  typical  boundary-element  scheme.  For  this  reason,  in  Part  II  of  this 
study,  we  return  to  the  scheme  where  the  parameter  d  is  finite  and  see  whether  we  can  expand 
the  system  of  integral  equations  to  accommodate  a  magnetic  current  density. 

One  remark  we  wish  to  make  is  that  (5.6)  can  also  be  obtained  using  the  approach  suggested 
by  Hansen  and  Yaghjian  (reference  5).  Here,  the  region  of  application  of  Green’s  second  identity 
is  the  interior  of  the  coaxial  line  and  the  entire  upper-half  space  except  for  the  region  occupied 
by  the  monopole  (Figure  3.1).  If,  with  it,  we  use  the  free-space  Green’s  function  in  (3.3)  we  get 

-J J(r) x Vg(r,r')dS  =  H(r')  ,  r'eflufl)  m  =  Sa  u Sh  u ac  u S  (8.1) 


where  crc  stands  for  the  metallic  part  of  the  .tv-plane.  Writing  this  compact  statement  out,  we 
have 


- 1  (r)  x  Vg(r,  r ')dS  -  J  J ,  (r)  x  Vg(r,  r ')dS  -  J  ib  (r)  x  Vg(r,  r ')dS 


-jjs(r)xVg(r,ryS 

s 


jH'Xr'),  r  'eD+ 
[H'(r') ,  r'eD 


(8.2) 


where 


J^(r)  =  zxHr(r),  recrc.  (8.3) 

In  passing,  we  mention  that,  in  (8.1),  there  should  also  be  present  an  integral  over  a  hemisphere 
with  center  the  origin  and  a  radius  that  tends  to  infinity.  In  the  limit,  however,  this  integral 
vanishes  because  the  radiated  fields  obey  the  Silver-Muller  radiation  conditions.  If  this  were  a 
problem  in  scattering  by  a  plane  wave,  as  in  reference  3,  then  we  would  have  to  use  a  stationery 
phase  method  to  show  that  the  integral  vanishes.  This  is  because  a  plane  wave  does  not  satisfy 
the  Silver-Muller  radiation  conditions.  Besides  the  scattered  fields,  the  sum  of  the  fields  of  the 
incident  and  reflected  plane  waves  would  appear  in  this  integral. 

If  r  e  a  and  we  pre-cross  (8.2)  with  the  unit  normal  along  the  z-axis,  then  the  surface  integral 
over  uc  becomes  zero  for  the  same  reason  as  (5.5),  and  we  get 


(8.4) 


17 


NAWCADPAX/TR-2013/1 15 


which  is  identical  with  (5.6).  We  note  that  this  is  the  only  occasion  in  which  the  surface  integral 
over  the  metallic  part  of  the  xy-plane  is  zero;  otherwise,  (8.2)  could  be  used  to  obtain  integral 
equations  over  the  whole  structure. 

Another  remark  is  that,  in  the  formulation  resulting  in  (8.1),  instead  of  the  free-space  Green’s 
dyadic  we  can  use  the  one  in  (2.2).  Because  of  the  property  (2.7),  the  integral  representation  will 
not  involve  a  surface  integral  over  ac .  Proceeding  as  in  Section  2,  we  get  in  place  of  (2.13) 

-  J  { ■ Js  (X) x  Vg(r,  r')  +  [J5(r)x  Vg(r,  r/)]  }  dS 

s 

-  { { Ja  («*) x  Vg(r,  r ')  +  [ Ja  (r)  x  Vg(r,  r/)]  }  dS 

s„ 

-  \  { J,  (r)  x  Vg(r, r ')  +  [Jb  (r)  x  Vg(r,  r/)],. }  dS  =  H  (r 1 ')  +  H,  (r/) ,  r'eD+uD  (8.5) 

sb 


where  H  stands  for  the  total  magnetic  field  in  either  of  the  two  regions.  This  statement  is  not 
entirely  correct.  For  example,  if  r'  =  ( p '  >  b,(p' ,z'  >  0) ,  then  r/  does  not  belong  to  the  region  of 

integration  D  and  the  last  tenn  in  (8.5)  should  be  equal  to  zero.  We  can  correct  for  this  by 
introducing  characteristic  functions  and  writing  a  separate  statement  for  each  of  the  two  regions. 
For  example,  in  the  upper  region 

-  j  {- Js  (X) x  Vg(r,  r')  +  [J5(r)x  Vg(r,  r/)]  }  dS 

s 

-  { {• 3 a  (r)  x  Vg(r,  r')  +  [  Ja  (r)  x  Vg(r,  r/)]  }  dS 

-  J{j4(r)xVg(r,r')  +  [Ji(r)xVg(r,r,')]f}dS  =  H'(r')  +  ^(r;)[H'(r;)]i  ,  r'sfl*  (8.6) 

sb 


where  %D  ■>  the  characteristic  function  of  the  region  D,  is  defined  by 


Xd{*) 


fl,  ref 
[0 ,  r  £D 


(8.7) 


Still,  we  have  to  worry  about  the  image  point  being  on  the  boundary  of  D  (the  walls  of  the 
coaxial  line).  We  can  write  a  statement  similar  to  (8.6)  when  the  observation  point  is  in  D,  with 
the  same  concerns. 

The  problem  with  an  expression  like  (8.6)  becomes  evident  when  we  try  to  obtain  integral 
equations.  The  image  of  a  point  of  S  may  correspond  to  a  point  somewhere  in  D.  Thus,  we  may 
generate  unknowns  not  only  on  the  metallic  boundaries  but,  also,  in  the  regions  D  and  D  .  This 
is  quite  unacceptable  and  we  abandon  this  approach. 
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CHAPTER  2 

EQUATIONS  EXTENDING  OVER  PART  OF  THE  WALLS  OF  THE  LINE 


20 


NAWCADPAX/TR-2013/1 15 


1.  INTRODUCTION 


In  Chapter  1,  we  derived  BIEs  for  the  geometry  of  Figure  1.1.  All  surfaces  are  perfectly 
conducting.  The  ground  plane  is  infinite  and  the  coaxial  line  extends  to  infinity  in  the  lower-half 
space.  The  inner  and  outer  radii  of  the  transmission  line  are  a  and  b,  respectively.  The  line 
supports  only  a  TEM  wave  which  means  that  (references  1  and  2 


(1.1) 


where  k  is  the  wavenumber  of  the  time -harmonic  ( c+IOJl )  electromagnetic  wave  in  the  line. 

Z 


t 


b 


Figure  1.1.  A  cylindrical  monopole  over  a  ground  plane,  fed  by  a  coaxial  line. 
The  excitation  of  the  line  occurs  at  z  =  -oo  and  results  in  a  TEM  wave  with  fields 


V  e  YV  e 

E?(r)  =  ^bV/U  H*(r)  =  — 

In  —  p  In  —  p 

K  o  )  v  a 


where  V  is  the  voltage  of  the  inner  conductor  with  respect  to  the  outer,  and  Y  is  the  free-space 
admittance.  We  have  also  employed  cylindrical  coordinates  ( p,(p,z ).  Since  the  line  is  not 
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infinite,  we  also  have  an  induced  wave  in  the  line  with  fields  { E! ,  H'  } .  The  total  fields  { E' ,  H' } 
in  the  coaxial  line  are  the  sum  of  the  generator  and  induced  fields 

E'=Eg+E\  H=Hg+H.  (1.3) 

In  the  upper-half  space  we  have  radiated  fields  { Er ,  H'  } . 

The  integral  equations  we  derived  in  Chapter  1  on  the  walls  of  the  coaxial  line,  extend  over  the 
entire  (infinite)  length  of  the  line.  Here  we  will  develop  a  system  of  integral  equations  that  will 
extend  over  a  portion  of  the  walls  of  the  line.  The  region  under  consideration  is  shown  in  Figure 
1.2.  Inside  the  line,  we  will  develop  integral  equations  extending  no  deeper  than  z  = -d  .  In  the 
process,  we  will  use  a  number  of  results  of  Chapter  1.  Equation  numbers  from  there  will  be 
preceded  by  the  Roman  numeral  I. 

z 


t 


b 


Figure  1.2.  The  region!)  of  application  of  Green’s  second  identity. 
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2.  INTEGRAL  REPRESENTATIONS  IN  THE  EXTERIOR  AND  INTERIOR 
For  convenience,  we  define  electric  and  magnetic  current  densities  on  rby 
Jr(r)  =  zxH'(r)  ,  Mr(r)  =  -zxE'(r)  ,  rer  (2.1) 

From  (1.4.2) 

2 J  [z  x  E'  (r)]  •  [  VVg(r,  r ')  +  k2g( r,  r ')l]  dS  -  z2kZ  J  j[nxH(  (r)]  x  Vg(r,  r')}  dS 

t  m 

-ikZ^ { J5 (r) x  Vg(r,  r ')  +  [ Js (r) x  Vg(r, r/)]  } dS  =  ikZH1  (r') ,  r'  e  D+ .  (2.2) 

s 

The  surface  w  is  defined  in  (1.3.1).  Taking  into  consideration  the  current  density  definitions 
(I.4.6)-(I.4.9)  and  the  ones  above,  we  write 

-2j  Mr  (r)  •  [Wg(r,  r ')  +  k2g(r,  r ')l]  dS 

T 

-ilkZ  J  3a  (r)  X  Vg(r,r')dS  +  J  Jb  (r)  x  Vg(r,r')dS  +  J  JCT  (r)  x  Vg(r,r')dS  +  J  JT  (r)  x  Vg(r,r')dS 

Sb  a  r 

—ikzj  { Js  (r)  x  Vg(r,  r ')  +  [  Js  (r)  x  Vg(r,  r/)]  }  dS  =  ikZH1  (r') ,  r'  e  D+ .  (2.3) 

s 

We  next  complete  (1.6.7)  by  supplying  the  integral  over  rfrom  (1.6.6) 
4=-jMr(r)-Vxr2(r,rVS 

T 

-ikz\  J  Ja  (r)  •  T2  (r,  r')dS  +  |  J4(r) -r2(r,rys 

[sa  s„ 

+J  JCT(r)  •  r2(r,rV5  +  J  Jr(r)  •  r2(r,r'K>  (2.4) 

(7  T  J 

Using  (1.6.5)  together  with  this,  we  get 

-ikZ  Jj8(r)-r2(r,ry5  +  Jj6(r)-r2(r,ry5  +  Jj<7(r)-r2(r,ry5  +  Jjr(r)-r2(r,ry5 

Uo  sb  a  r 

Mr(r)  -V  x  r2(r,r')t/1Sl  =  k2ZH'(r')  ,  r'efl.  (2.5) 

T 
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By  the  definition  (1.6.2)  of  the  dyadic,  the  last  integral  in  this  becomes 

j Mr (r) •  V x r2 (r, r' ')dS  =  -ik jMr (r) •  { VVg(r, r')  +  k2g(r, r')I - V [ Vg(r, tf].  - k2g(r, r,%}dS 

T  T 

=  -ik J Mr (r)  •  V { Vg(r, r ')  -  [Vg(r, r/)]  }  dS  -  ik}  J  jg(r,  r ')Mr (r)  -  g( r,  r/) [Mr (r)].  }  dS 

T  T 

(2.6) 

For  the  first  integral  on  the  right,  we  use  the  divergence  theorem  to  write 
jM,(r)-V{vg(r,r')-[Vg(r,r')],}« 

T 

=  j  { V  ‘  [Mr  (r)  ( Vg(r ,  r ')  -  [Vg(r,  r/)]. )]  -  V  •  Mr  (r)  ( Vg(r ,  r ')  -  [Vg(r,  r/)];.  )|  dS 

T 

=  \v- Mr (r) ( Vg(r, r ')  - [ Vg(r, r/)]/ ) ds  - { V  •  Mr (r) ( Vg(r, r ')  - [ Vg(r, r/)]  ) dS  .  (2.7) 

Y  * 

The  contour  /is  composed  of  the  two  contours  that  constitute  the  boundary  of  z.  The  unit  normal 
v  on  /  is  equal  to  the  negative  of  the  normal  on  each  of  the  two  walls  of  the  coaxial  line.  The 
component  of  the  magnetic  current  density  along  it  is,  by  (2.1),  the  tangential  component  of  the 
total  electric  field  on  each  wall.  Its  value  there  is  zero.  Thus,  the  contour  integral  is  equal  to  zero, 
and  we  have  that 

j  Mr  (r)  •  V  x  r2  (r,r')dS 

T 

=  I'qv.M,(r)(Vg(r,r')-[Vg(r,r,')],)rfS-tt3Jh(r,r')M,(r)-^(r,0[M,(r)]i.}rfS.(2.8) 

T  T 

Substitution  of  this  and  (1.6. 1 1)  in  (2.5)  gives 

-j  { (r)  x  Vg(r,  r ')  -  [  J8  (r)  x  Vg(r,  r/)]  }  dS  -  J  { Jh  (r)  x  Vg(r,  r ')  -  [  J,  (r)  x  Vg(r,  r/)]  }  dS 

s„  s„ 

~  2  j  Ja  (r)  x  Vg(r,  r  ')dS  -  J  { Jr  (r)  x  Vg(r,  r ')  -  [  J r  (r)  x  Vg(r,  r')]  }  dS 

(7  T 

+  y  J  V  •  Mr  (r)  ( Vg(r, r ')  -  [ Vg(r,  r/)]  )  dS  -  ikl\  {g(r,  r ')Mr  (r)  -  g(r, r/)  [Mr  (r)], }  dS 

T  T 

=  H'(r') ,  r'eD.  (2.9) 

This  is  the  total  magnetic  field  representation  in  D.  We  proceed  to  get  one  for  the  total  electric 
field.  We  begin  with  Green’s  second  identity 

J  {H'  (r)  •  [V  x  V  x  Tl  (r, r ')]  -  [v  x  V  x  H'  (r)]  •  T,  (r,  r')}  dV 

D 
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=J{["X  H' (r)] •  V x r, (r, r ')  +  ikY [n x E' (r)] •  T, (r, r')| dS  ,  r'eD.  (2.10) 

zu 

The  Green’s  dyadic  of  the  first  kind  is  defined  in  (1.2.2).  For  the  volume  integral,  we  get 

ID  =-k2YE'(r'),  r'eD.  (2.11) 

From  (1.2.2) 

J[nx  H'  (r)]  •  V  x  r,  (r,r')dS 

m 

=  -ik\[n  x  H'(r)]  •  v[Vg(r,r')  +  (Vg(r,r/))J^ 

m 

-ik^\n  x  H'(r)]  •  [g(r,r')I  +  g(r,iv')II.]c?5 

m 

=  ik\  V  •  [n  X  H'(r)][Vg(r,r')  +  (Vg(r,r/'))/]^  . 

m 

-ie\[n  X  H'(r)]  •  [g(r,r')I  +  g(r,rfc]dS  .  (2. 12) 


The  surface  divergence  theorem  has  been  applied  above  to  the  closed  surface  m.  The  result  is,  of 
course,  zero.  For  the  remaining  term  in  (2.10) 

J  [n  x  E'  (r)]  •  r,  (r,  r r)dS  =  -{ Mt  (r)  •  T,  (r ,  r ')dS  =  ik\ Mr  (r)  •  V  x  [g( r,  r')I  +  g(r,r%  ]  dS 

m  t  t 

=  ik\ {Mr  (r)  x  Vg(r, r ')  +  [Mr  (r)  x  Vg(r, r,!)]. }  dS .  (2.13) 

T 

In  arriving  at  this  result,  we  have  taken  into  consideration  the  fact  that  the  total  electric  field  is 
zero  on  the  walls  of  the  coaxial  line  and  that  the  tangential  component  of  the  dyadic  is  equal  to 
zero  on  the  z-  0  plane.  Substituting  the  above  results  in  (2.10),  we  get 

-^jv{nxH'(r)][Vg(r,rO  +  (Vg(r,r/));.]^  +  /kzj[nxH'(r)]-[g(r,r')I  +  g(r,r/)I,.]^ 

m  m 

+  j {Mr  (r)  x  Vg(r,  r ')  +  [Mr  (r)  x  Vg(r,  r/)]; }  dS  =  E'  (r 1 ') ,  r 1 '  e  D .  (2.14) 

T 

Moreover,  breaking  the  first  and  second  integrals  into  their  constituent  parts 

- ^  f v  •  [J. (r)][ V«<r, r'>  +  (’ Vgfr.r/)), ] dS  +  f v  ■  [ Js (r)][ Vg(r, r')  +  (' Vg(r, r/))f ] dS 

k  { Sa  Sb 
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+  J  V  •  [  Jr  (r)]  [ Vg(r,  r ')  +  ( Vg(r,  r/))f  ]  dS  +  2 J  V  •  [  JCT  (r)]  Vg(r,  r ' ')dS 

T  <7  ) 

+  ikZ  { {g(r,  r ') Jfl  (r)  +  g(r,  r/)  [ Ja  (r)], }  +  J {g(r,  r ') J,  (r)  +  g(r,  r/)  [ J,  (r)], }  dS 

U  5 b 

+  J  {g(  r,r')JT  (r)  +  g(r ,  r/)  [JT  (r)],  }dS  +  l\  g(  r,  r ')  JCT  (r)<tf 

r  <j  J 

+  j  {Mr  (r)  x  Vg(r,  r ')  +  [Mr  (r)  x  Vg(r,  r/)].  =  E'(r') ,  r'eD.  (2.15) 

T 

This  is  the  integral  representation  for  the  electric  field. 
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3.  INTEGRAL  EQUATIONS  ON 
From  (1.5.1)  and  (2.10),  we  get  that 

-z  x  |  { Jfl  (r)  x  Vg(r,  r ')  -  [Ja  (r)  x  Vg(r,  r/)] |. }  dS 

sa 

-z  x  { { 3b  (r)  X  Vg(r,  r')  -  [ Jb  (r)  x  Vg(r,  r/)]  }  dS  -  2z  x  J  Ja  (r)  x  Vg(r,  r')dS 

Sb  a 

-z  X  |  { Jr  (r)  X  Vg(r,  r ')  -  [ Jr  (r)  x  Vg(r,  r/)]  }  dS 

T 

+  ^ z  x  J  V  •  Mr  (r)  ( Vg(r,  r ')  -  [ Vg(r,  r/)], )  dS 

T 

- ikYz  x  J |g(r,r')Mr(r)  -  g(r, r/) [Mr (r)]; }  dS  =  Jr (r')  ,  r '  e  r  . 


We  can  simplify  this  expression  considerably.  For  the  third  tenn,  we  write 
-2z  x  |  JCT (r) x  Vg(r,r')dS  =  -2 j  Ja (r)  dS 


<7 


<7 


while  for  the  fourth 


-z  x  { Jr (r)  x  Vg(r,r ')  -  [ J r (r)  x  Vg(r,  r')]  }  =  -z  x  [ Jr (r)  x  Vg(r, r')] 


#g(r,r')  ,  „  r/)' 


+  z  X  {[z  X  H'  (r)]  X  Vg(r,  r/)j  =  - Jr  (r)  Cg(X,r  ’  +  z  x  J  H  (r) 


=  -Jr(r) 


Sg(r,r')  dg(  r,r/) 
dz  dz 


=  Jr(r) 


dz 

og(r,r') 

dz 


dz 


Thus,  in  place  of  (3.1),  we  can  write 

-z  x  |  { Jfl  (r)  x  Vg(r,  r ')  -  [ Ja  (r)  x  Vg(r,  r/)]  }  dS 

Sa 

-  £  x  J  {  JA  (r)  x  Vg(r ,  r ')  -  [  JA  (r)  x  Vg(r , r/)]? }dS-l\ (r ) dS 

S„  CT  dZ 

+  J  Jr  (r)  r/)  ^  +  y  J  V  •  Mr  (r)z  x  ( Vg(r,  r ')  -  Vg(r ,  r/))  dS 

T  T 

-ikYzx  jMr(r)[g(r,r')-g(r,r/)]^  =  jjr(r'),r'er. 


(3.1) 


(3.2) 


(3.3) 


(3.4) 


27 


NAWCADPAX/TR-2013/1 15 


This  is  the  integral  equation  for  the  electric  current  density  on  r.  For  the  magnetic  current 
density,  we  have  from  (2.15) 


^  x  j  V  •  [ J„ (r)][Vg(r, r ')  +  (Vg(r, r/)), ] dS  +  {  V  •  [ J, (r)] [  Vg(r,  r ')  +  ( Vg(r, r/))f  ] dS 


+  J  V  •  [  Jr  (r)]  [ Vg(r,  r ')  +  ( Vg( r ,  r/)),  ]  dS  +  2 J  V  •  [  JCT  (r)]  Vg(r,r')dS 


+  ikZz x <  | {g(r,r')Ja(r)  +  g(r,r;)[jfl(r)].}r/5  +  |{g(r,r')Ji(r)  +  g(r,r/)[J6(r)].}r/5' 


+  J  {g(r ,  r ')  Jr  (r)  +  g(r ,  r/)  [  Jr  (r)],  }dS  +  l\  g(r,  r ')  JCT  (r)<tf 

r  <t  J 

+  z  X  J  {Mr  (r)  X  Vg(r,  r ')  +  [Mr  (r)  x  Vg(r,r/)]. }  dS  =  -  ^Mr  (r 1 ')  ,  r'er. 


We  can  re-write  this  as 


Jv.[J„(r)]lx[Vg(r,r’)  +  Vs(r,r/)]rfS+Jv.[J,(r)]lx[Vg(r,r')  +  Vg(r,r;)]rfS 


+ 1 V  ■  [J,(r)]zx  [Vg(r.r')  +  Vg(r,r/)]<iS  +  2j  V  [J„(r)]  zxVg(r,r')dS 


+  ikZ  |  J  z  x  Ja  (r)  [g( r,  r ')  +  g(r,  r/)]  +  J  z  x  J,  (r)  [g(r,  r ')  +  g(r,  r/)]  dS 


+  J  z  x  JT  (r)  [g(r,  r ')  +  g(r,  r/)]  dS*  +  2j  g(r,r')z  x  Ja  (r )dS 


+  f  M»  ^r’r/)  JS  =  -^(r') ,  r'er. 

J  <7Z  7 
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4.  SYSTEM  OF  INTEGRAL  EQUATIONS 

Besides  the  integral  equations  (3.4)  and  (3.6),  the  system  of  equations  will  include  equations 
from  the  rest  of  the  structures  that  make  up  the  problem.  Using  the  procedures  of  Chapter  1,  we 
get  from  (2.9)  that 

-  n’  x  |  {, Ja (r)  x  Vg(r ,  r ')  -  [ Ja (r)  x  Vg(r, r/)]  }  dS 

s„ 

-  n'  x  { { J,  (r)  x  Vg(r,r')  -  [jb  (r)  x  Vg(r, r/)]( }  dS  -  2h'  x  {  Ja (r)  x  Vg(r,  r')dS 

St  o’ 

-  n  x  J{ Jr (r)  x  Vg(r.r')  -  [ J,  (r)  x  Vg(r,  r/)], }  dS 

T 

+  ^n'xJv-Mr(r)(Vg(r,r')-[Vg(r,r/)]()^ 

k  T 

-ikYh'x  J{g(r,r')Mr(r)-g(r,r,')[Mr(r)]/}rf5  =  ^Ja(r')  ,  r '  e  Sa  (4.1) 

T 

and 

-  /)'  x  |  { Jfl (r)  x  Vg(r,r')  -  [ Ja (r)  x  Vg(r, r/)]. }  dS 

s„ 

-«'x{{j,(r)x  Vg(r,r')  - [Jb (r)  x  Vg(r , r')]( } dS  -  2h' x  {  Ja (r) x  Vg(r, r')dS 

sb  <J 

-  n  x  J{ Jr (r)  x  Vg(r.r')  -  [ J,  (r)  x  Vg(r,  r/)], }  dS 

T 

+  ^77'xJv-Mr(r)(Vg(r,r')-[Vg(r,r/)]()^ 

k  T 

-  ikYti  x  J {g(r,r')Mr  (r)  -  g(r,r/)  [Mr(r)],}  dS  =  Ub(r')  ,  r ’<=Sb.  (4.2) 


From  (1.4.2) 

ikYzx  [Mr(r)g(r,r')<i1S'  +  — zx 
{  k 


J  Mr(r)  •  Wg(r,r')dS 


-fxJjs(r)xVg(r,rV*S,  =  -Jff(r')  ,  r'  e  cr 

s 


(4.3) 
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and 


2 ikYh'  x  J Mr  (r)g(r,r +  2  —  h'  x  |Mr  (r)  •  VVg(r,r')^ 

k 

T  T 

-M'xfj,(r)x  Vg(r,  r ' ')dS  +  J  J,  (r)  x  Vg(r,r')dS 

U  sb 

+  J  K  (r)  x  Vg(r,r')dS  +  J  Jr  (r)  x  Vg(r,r')dS 

<J  T  J 

~  n'  x  J  { Js  (r)  x  Vg(r,  r ')  +  [ J s  (r)  x  Vg( r,  r/)]  }  dS  =  ]-  Js  (r ' ') ,  r'  e  5 


(4.4) 
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5.  REMARKS 

We  have  some  concern  here  as  to  whether  the  equations  we  just  derived  are  the  right  ones  for 
determining  the  unknown  currents.  We  will  not  know  the  answer  to  this  until  we  compare 
numerical  results  between  this  system  and  the  one  of  Chapter  1.  If  this  system  provides  the 
correct  answers,  then  it  is  the  preferred  one  because  it  also  provides  substantial  computational 
savings  if  we  keep  d  small.  We  note  that,  in  this  fonnulation,  current  densities  are  assumed 
piecewise  differentiable  in  the  interior  of  the  line.  In  Chapter  1,  it  is  sufficient  that  these  currents 
are  continuous. 

Equations  (2.5)  and  (2.14)  give  the  fields  in  the  interior  of  the  coaxial  line.  They  are  the  result 
of  applying  Green’s  second  identity  in  the  region  bounded  by  the  closed  surface  m,  defined  in 
(1.3.1) 

m  =  (jUTyjSauSb.  (5.1) 

We  make  two  remarks.  The  first  is  that  the  integral  over  the  surface  t(z  -  -d )  can  be  converted 
to  an  integral  over  the  two  cylindrical  surfaces  that  join  this  surface  to  the  identical  disk  at 
z  =  -oo .  Indeed,  the  walls  and  the  two  disks  make  up  a  closed  surface.  Application  of  the 
appropriate  Green’s  second  identity  (in  each  of  the  two  cases),  with  the  observation  point  outside 
this  surface  ( -d  <  z  <  0  ),  leads  to  the  result  that  the  integral  over  r  is  equal  to  the  integral  over 
the  other  three  surfaces.  As  remarked  in  I,  the  integral  over  the  disk  at  infinity  is  equal  to  zero; 
hence,  the  conclusion  that  the  integral  over  z  is  equal  to  the  same  integral  over  the  cylindrical 
walls  extending  from  -  d  to  -oo .  In  the  process,  the  terms  involving  the  circumferential 
components  of  the  electric  field  drop  out,  because  the  walls  are  perfectly  conducting,  and  we 
end-up  with  the  same  representations  as  in  Chapter  I.  We  see,  then,  that  we  can  start  with  the 
present  case  and,  as  a  corollary,  derive  the  equations  in  Chapter  I. 

The  second  remark  is,  that  when  we  use  the  present  approach  to  compute  the  current  densities 
on  G7,  then  we  can  use  the  integral  representations  to  compute  the  fields  in  the  coaxial  line  only 
in  the  region  enclosed  by  m.  Differently  stated,  we  are  not  able  to  compute  the  fields  below  depth 
-d.  In  I,  we  can  compute  the  fields  at  any  depth  inside  the  line.  Perhaps,  this  is  the  trade-off 
between  the  two  methods.  We  are  confident  that  the  present  approach  is  correct  and  that  it  is  the 
preferred  one  computationally  since  what  we  are  striving  for  is  to  get  the  fields  correctly  at  the 
opening,  cr,  of  the  coaxial  line. 
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6.  FAR  FIELD 

We  will  present  two  ways  to  compute  the  far  field.  First,  we  point  out  that  we  can  express  the 
far  field  in  terms  of  an  integral  over  cr  (and  an  integral  over  S).  The  integral  over  cr,  however, 
contains  the  magnetic  current  density  on  cr,  a  quantity  we  have  not  computed  above.  We  are  thus 
led  to  consider  other  representations  of  the  far  field.  These  representations  will  be  for  the 
magnetic  field.  To  obtain  representations  for  the  electric  field,  we  recall  the  following. 

In  the  far  field,  the  magnetic  field  can  be  expressed  in  the  fonn 

H' (r )~- — f(r)a(r),  r  — »  oo.  (6.1) 

r 


From  Maxwell’s  equations 


E'(r)  =  — VxH'(r) 
ikY 


Q~ikr  ! 

x  a(r )  +  — — /(r)V  x  a(r)  > , 


r  — >  co  . 


(6.2) 


But 


V 


xa(r) 


-ik— — f(r)rxa(r)  +  C)[r  2), 


r  — >  oo 


(6.3) 


and 

Vxff(r)  =  0[r~1^  ,  r  — >  co  .  (6.4) 

Hence 


c 

E'  (r) — Z - f(r)fxa(r),  r  — >  oo 

r 


(6.5) 


or 


Er  (r)  =  -Zr  x  Hr  (r)  ,  r  — >  oo . 


(6.6) 


We  turn  to  the  derivation  of  the  far-field  representations  of  the  magnetic  field.  The  first  comes 
from  (2.3).  We  observe  that 


g(  r,r')~- 


—ikr'  ^ikr-r 

e  e 
4  Tir 


■  00 


(6.7) 
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Vg(r,r')  — ikg(r,r')r'  ,  VVg( r,r') - k2g(r,r')r'r' ,  r'  — >  oo  (6.8) 

and  substitute  in  (2.3) 

ikZl r  (r')  ~  -2j Mr  (r)  •  [-&2g(r,  r')r'r'  +  &2g(r,  r')I  J  dS 


+  2k2 Zr' ’  x  J  Ja  (r)g(r,r')dS  +  J  J,  (r)g(r,  r')dS  +  J  JCT  (r)g(r,  r ')</S  +  J  JT(r)g(r,r')dS 


-k2Z\  Js(r)xr'g(r,r')  + 


Js(r)x^g(r,r;') 


dS  ,  r'  — >  oo  . 


(6.9) 


We  observe  that 

Mr  (r)  •  [-r 'r'  + 1]  =  -Mr  (r)  •  [r '  x  (r'  x  I)]  =  [r'  x  Mr  (r)]  •  (?'  x  I)  =  -r'  x  [r'  x  Mr  (r)]  (6.10) 


and  that 


Js(r)x^. 


=  x 


( y  ■  Js  (r)) [-Z  •  r  I  -  (z  •  Js  (r)) (y  ■  r 


+  y 


(z-J5(r)) (x-r  )-(x-Js(r))  -z 


-  z 


=  X 


y-r  )(-z-J5(r))-(z-r  )(j>- Js(r)) 


+y 


(x-Js(r))^v-r'j-(v-Js(r 

J,(r))-(x-r  )(-z-Js(r)) 


x  •  r 


z-r  1 1  x 


+  z 


x  •  r 


J5(r))-(v-^'  Js(r))  =r  x[js(r)].. 


(6.11) 


Substitution  of  the  last  two  results  in  (6.9)  gives 


H'  (r') - ilkYr' : 


5'xjMr(r)g(r,rV5 


ilk?  x  J  Ja  (r)g(r,r')dS  +  J  Jb  (r)g(r,r')dS  +  J  JCT  (r)g(r , r ' ')dS  +  J  Jr  (r)g(r,  r ' ')dS 

U  Sb  a  v  J 

ikr'  x|{j5(r)g(r,r')-[J5(r)].g(r,r,')}^,  r'->co  (6.12) 


where  g  is  given  by  (6.7). 

A  second  representation  of  the  far  field  can  be  obtained  from  (1.6.1)  with  T2  replaced  by  T, 
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J  {E'  (r)  •  [V  x  V  x  r,  (r,  r')]  -  [v  x  V  x  E'  (r)]  •  T,  (r,r')}  dV 

D 

=  J  j-E'  (r)  •  [n  xVxT,  (r,  r')]  +  [n  xVxE'  (r)]  •  r,  (r,  r')}  dS ,  r'  e  Df .  (6.13) 

zn 

The  region  Df  is  the  upper-half  space,  excluding  the  monopole  structure  and  the  reflection  of 

the  region  D  about  the  xv-plane.  We  can  think  of  it  as  representing  the  far- Held  region.  Thus,  the 
contribution  of  the  volume  integral  is  zero  while,  for  the  surface  integral,  we  have 

|  j-E'  (r)  •  [77  x  V  x  T,  (r,  r ')]  +  [77  x  V  x  E'  (r)]  •  r,  (r,  r ')}  dS 

m 

=  -J  {[z  X  E'  (r)]  •  V  X  r,  (r,  r')i  dS 

<J 

+ !{[,-,  x  E'  (r)]  ■Vxri(r,r')- ikz[n  x  H'  (r)]  •  T,  (r,  r ')}  dS 

T 

-  ikz\  [77  x  H' (r)]  •  T,  (r, r')dS  - ikZ^[h  x  H' (r)]  •  r,  (r, r')dS  (6. 14) 

sa  s„ 

or 

-J{[ixr  (r)]  •  V  x  Yl  (r,  r ')}  dS  -  J  {M,  (r)  •  V  x  T,  (r,  r ')  +  ikZJT(r)  •  T,  (r,  r')}  dS 

(7  T 

-  ikZ  J  J5  (r)  •  r J  (r,  r  ')dS  -  ikZ  J  (r)  •  (r,  r')dS  =  0  ,  r '  e  Df .  (6.15) 

sa  Sb 

We  combine  this  with  (1.2.12)  to  get 

-1 kzj  Js  (r)  •  (r,r')dS  -  J  {M,  (r)  •  V  x  T,  (r,  r')  +  ikZ3T  (r)  •  (r,  r ')}  dS 

S  T 

-  7lZ  J  J5  (r)  •  r,  (r,  r')dS  -  ikZ  j  JSh  (r)  •  (r,r ')dS  =  k2ZH'  (r') ,  r  ’eDf  (6.16) 

s»  s„ 

or,  in  terms  of  the  scalar  Green’s  function  (see  (1.2.2) 

IT  (r')  =  -J{  Js  (r)  x  Vg(r,  r')  +  [js(r)  x  Vg(r,  <)],  j  dS 

s 

-IK  (r)  x  Vg(r,  r ')  +  [  JSa  (r)  x  Vg(r,  r/)]. }  dS 

-|K (r)  X  V^^r’  r'}  +  [Js>  (**)  x  V^(r’  r')]( }  dS 

sb 

- 1{ Jr (r)  x  Vg(r, r')  +  [ Jr (r)  x  Vg(r, r/)]  } dS 

T 
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1 

ikZ 


J  {Mr  (r)  •  V  [Vg(r,  r ')  +  (  Vg(r,  r/)),  ]  +  k2  [Mr  (r)g(r ,  r ')  +  (Mr  (r));  g(r ,  r/)]}  dS  , 


r' e  D f .  (6.17) 


We  proceed  to  evaluate  this  expression  in  the  far  field  using  (6.8) 

H'  (r')  ~  -ikr'  x|{js(r)g(r,r')  +  [Js(r)].  g(r,r/)}</£ 

s 

-  ikr  x  { {jSa  (r)g(r.r')  +  [jSo  (r)].  g(r,r/)}</S  -ikr  x  J  {jSj  (r)g(r,r')  +  [j5j  (r)].  g(r,r/)J  dS 

sa  sb 

-  ikr  x  |  Jr  (r)  [g(r,  r ')  +  g(r,  r/)]  -  ikYr  x  jr  x  J  Mr  (r)  [g(r,  r')  +  g(r,  r/)]  dS^  ,  r'  -»  oo 

(6.18) 

where  g  is  given  by  (6.7).  To  obtain  this  result,  we  used  (6.10)  and  (6.1 1).  We  also  used  the  fact 
that,  on  t,  the  electric  and  magnetic  current  densities  are  transverse  to  the  z-axis.  As  a 
consequence 


Finally, 


(6.19) 


V[Vg(r,r/)].  ~  V 


-ikg(  r,r,V, 


=  -/lV[g(r,r/)f']  ~  -k2g(r,r;)r^i 


oo 


(6.20) 


The  far  field  expression  (6.18)  differs  substantially  from  the  one  in  (6.12).  It  may  even  be 
argued  that  (6.12)  is  simpler  to  compute.  Equation  (6.18)  demonstrates,  however,  that,  when  it 
comes  to  the  far  field,  we  have  an  image  theory  in  place.  We  proceed  to  show  this.  Define  a 
current  density  Js  on  the  image  5)  of  S  (about  the  xy-plane)  as  follows 


zxjs.  (r)  =  -zxjv(r),  r  g  Si 
z-Js_(r,)  =  z-Js(r),  re5  )' 

This  statement  is  in  agreement  with  how  electric  currents  are  imaged  in  the  presnce  of  a  perfectly 
conducting  half-space  with  normal  along  the  z-axis.  We  make  similar  definitions  on  the  images 
of  the  rest  of  the  surfaces  in  (6.18).  For  the  magnetic  current  density,  we  define  its  image  on  r; 
as  follows 


zxMr.(r,.)  =  zxMr(r)  ,  rer. 


(6.22) 
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We  recall  from  (6.19)  that  the  magnetic  current  density  does  not  have  a  component  along  the  z- 
axis. 


We  examine  the  second  term  in  the  first  integral  in  (6. 18).  By  (6.7) 


j[Js(.r)lg(.r  >*!)<£  = 


dS. 


(6.23) 


Using  (6.21)  and  the  fact  that 
f  •  r  =  r  •  r. 


(6.24) 


we  get  that 

~—ikr' 

J[JS  (r)].  g(r,r l)dS  =  —  j  J,  (r, )  e"'"  dS, . 

s  s 


(6.25) 


This  is  an  integral  defined  over  the  image  of  the  surface  S.  It  is  worth  writing  the  whole 
expression  (6.18)  out 


H'  (r')  ~  -ik 


y  4 nr'  j 


-  ik 


-  ik 


y  4  nr'  j 
v  Anr'  j 


r  x  { 


r  x  < 


'x  (js(r)  Qm"r  dS  -  J  Js.  (r. )  Qikf'r>  dS \ 
[s  s, 

Jjs<(r)e»r<K-  J 


|jSi(r dS-  |  Jw(r,)elf'^ 

Sb  (sb\ 


-ik 


y  Anr'  j 


x  J  Jr  (r)  Qikf'  r  dS  -  J  Jr  (r, )  eikf'  r‘  dS, 


ikY 

6 

^  Anr'  j 

r~  t 

r  x 

r  x 

| Mr  (r)  Qikf"  r  dS  +  \ Mr_  (r, )  e^"'  dS, 

T  T; 

> ,  r  — >  oo  , 


(6.26) 


Thus,  we  have  an  expression  for  the  far  field  in  tenns  of  integrals  over  the  surfaces  S,  Sa,  Sh ,  z 
and  their  images  about  the  xv-plane.  We  note  that  there  is  no  integral  over  cr. 
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CHAPTER  3 

MONOPOLE  AS  AN  EXTENSION  OF  THE  CENTER  CONDUCTOR 
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1.  INTRODUCTION 

In  this  chapter,  we  use  results  from  Chapters  1  and  2  to  develop  a  system  of  integral  equations 
for  geometries  as  in  Figures  1.1  and  1.2.  In  these  figures,  the  monopole  is  an  extension  of  the 
center  conductor  of  the  coaxial  transmission  line.  The  diameter  of  the  monopole  in  Figure  1.2  is 
not  restricted  in  any  way.  It  can  be  smaller  than  that  of  the  center  conductor  or  even  greater  than 
that  of  the  outer  conductor.  Also,  the  height  h  of  the  monopole  above  ground  is  variable  and  it  is 
allowed  to  become  zero.  Figure  1.1  is  a  special  case  of  Figure  1.2.  We  proceed  to  derive  integral 
equations  for  the  geometry  of  Figure  1.2.  Equation  numbers  from  Chapters  I  and  II  are  preceded 
by  Roman  numerals  I  and  II,  respectively. 


z 


b 


Figure  1.1.  Monopole  as  extension  of  center  conductor  of  coaxial  line 
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z 


t 


- > 

b 


Figure  1.2.  A  more  general  monopole  geometry 

For  h  >  0 ,  we  denote  the  surface  of  the  structure  in  the  upper-half  space  ( z  >  0  )  by  S.  If  h  -  0 
,  then  we  consider  the  part  that  lies  on  the  xy-plane  as  part  of  the  inner  conductor  Sa.  We  first 
consider  the  case  h  >  0 . 
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2.  CENTER  CONDUCTOR  EXTENDING  ABOVE  THE  PLANE  (h>  0) 

The  unknowns  are  given  in  (I.4.6)-(I.4.9).  The  currents  J  s  and  J  s.  must  be  equal  on  the  xy- 
plane.  The  integral  equations  are  given  by  (1.5.2),  (1.5.6),  (1.7.1),  and  (1.7.2). 

If  the  coaxial  line  is  to  be  terminated,  as  in  Figure  2,  then  we  have  two  additional  unknowns 
and  they  are  given  by  (II.2.1).  The  integral  equations  are  given  by  (II. 3.1),  (II. 3.2),  and  (II.4.1) 
through  (II.4.4). 
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3.  CENTER  CONDUCTOR  EXTENDING  TO  THE  PLANE  (h  =  0) 
The  unknowns  and  integral  equations  are  as  in  the  previous  section. 
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4.  REMARKS 

We  must  make  sure  that  at  the  point  r'  =  (a,(p', 0),  integral  equations  (1.5.2)  and  (1.7.1)  give 
the  same  answer.  These  equations  are 

-2  n  x  (jjr)x  Vg(r,r')dS  +  J  3h  (r)  x  Vg(r,r')dS  +  J  JCT  (r)  x  Vg(r,r')dS 

U 

-«'x j{js(r)xVg(r,r')  +  [Js(r)xVg(r,r/)].}^  =  ^Js(r')  ,  r 'eS.  (4.1) 

s  2 

and 

-n '  x  |  { Ja  (r)  x  Vg(r,  r ')  -  [Ja  (r)  x  Vg(r,  r/)] . }  dS 

sa 

-  n'  x  |  { J,  (r)  x  Vg(r ,  r ')  -  [Jb  (r)  x  Vg(r,  r/)]  }  dS 

s„ 

-2 n'x  f  Ja (r) x Vg(r,r ')^  =  (O  ,  rTS,  (4.2) 

J  2 

a 

If  r'  =  {a,(p’,  0) ,  these  equations  become 

-2  p'  x  J  Ja  (r)  x  Vg(r,r')dS  +  J  Jh  (r)  x  Vg(r,r')dS  +  {jcr(r)x  Vg(r,  r')dS 

U  sb 

-  /3'x  J{js(r)x  Vg(r,r')  +  [J5(r)x  Vg(r,r')].}^  =  '  J  v  (r') ,  r'  =  (a,<p',0)  (4.3) 

5  2 

and 

-P'x|{j.(r)x  Vg(r,  r’)  -  [j„  (r)  x  Vg(r,  r’)]_ }  dS 

Sa 

-  p'  x  |  { Jb  (r)  x  Vg(r,  r')  -  [  Jb  (r)  x  Vg(r,  r ')]. }  dS 

sb 

~2p' x  [  JCT  (r)  x  Vg(r,r')dS  =  ^  Ja(r') ,  r'  =  (a,p',0).  (4.4) 

J  2 

<T 

The  p'  -component  of  both  equations  is  zero.  From  (4.3) 

V  •  Js  (r 1 ')  =  2z  •  J  Ja  (r)  X  Vg(r,  r')dS  +  J  Jh  (r)  x  Vg(r,r')dS  +  J  Jff  (r)  x  Vg(r,r')dS  (4.5) 

1  U  s>>  J 
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while  from  (4.4) 

V  •  J„  (r ')  =  2 z  •  J  JCT  (r)  x  Vg(r,r')dS  +  2z  •  {  Ja  (r)  x  Vg(r,  r')dS  +  J  J,  (r)  x  Vg(r,r')dS  . 

1  o  U  Sb  J 

(4.6) 

From  these  two 

=  r'  =  (a,«>',0).  (4.7) 

For  the  z-component,  we  get  from  (4.3)  that 

^  •  J s  (r')  =  -2 <p'  ■  J  Js(r)  x  Vg(r,r')dS 

2  s 

-2*'-  Jja  (r)  x  Vg(r,  r ')</S  +  J  Jfe  (r)  x  Vg(r,  r ')</£  +  J  Jff  (r)  x  Vg(r,r')dS  (4.8) 

U  Sb  a  J 


while  from  (4.4) 


a  (r')  =  -2^'  ■  J  JCT  (r)  X  Vg(r,  r')dS  =  -2j  Ja  (r)  •  [<p'  xV'g(r,r')]^ 

(7 

<?g(r,r')  ,  ~,^g(r,r') 


=  -2/J,(r). 


— z- 


dp 


—  +  P 


dz 


dS  =  0. 


(4.9) 


With  this  last  result,  we  can  write  for  (4.8) 


z  •  J5  (r ')  =  -4 <p’  •  {  Ja  (r)  x  Vg(r,  F ')dS  +  J  Jb  (r)  x  Vg(r,r')dS  +  J  Js  (r)  x  Vg(r,r')dS  . 

U  s>  $  J 

(4.10) 

From  (1.8.4) 

p'  •  JCT  (r ')  =  -$  •  J  Ja  (r)  x  Vg(r,r')dS  +  J  Jh  (r)  x  Vg(r,r')dS  +  J  J,  (r)  x  Vg(r,  r')dS  , 

U  sb  s  J 

r'=(a,<p',0).  (4.11) 


From  the  last  two 

z- Js(r')  =  4p'-JCT(r'),  r'  =  (a,<p',0) .  (4.12) 
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But  p'  ■  JCT(r')  is  the  normal  component  of  an  electric  current  density  on  a  perfect  conductor 
and,  hence,  it  must  be  zero.  We  conclude  then  that 

z-Js(r')  =  0,  r'  =  (a,^',0).  (4.13) 
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PART  2 

INTEGRAL  EQUATIONS  FOR  THE  FIELDS  OF  A  COAXIAL  LINE  RADIATING 

INTO  A  HALF-SPACE 


ABSTRACT 

This  is  the  second  Part  in  a  report  on  the  formulation  of  the  problem  of  radiation  of  a  coaxial 
line  into  a  half-space  in  terms  of  BIEs.  In  it,  we  use  the  results  of  Part  1  to  derive  BIEs  for  the 
fields  of  a  coaxial  line  radiating  into  an  otherwise  empty  half-space.  The  relationship  of 
wavelength  to  the  radii  of  the  line  is  such  that  the  input  wave  is  a  TEM  wave.  We  take  advantage 
of  the  circular  symmetry  of  the  line  to  reduce  the  vector  integral  equations  to  three  scalar 
equations.  The  unknown  electric  current  densities  on  the  walls  of  the  line  are  expressed  as 
infinite  series  in  the  natural  modes  of  the  line,  the  coefficients  of  the  modes  being  the  unknowns. 
We  point  out  that  any  of  the  three  integral  equations  can  be  solved  numerically  using  the  method 
of  weighted  residuals.  This  method  results  in  an  infinite  system  of  linear  equations  with  an 
infinite  number  of  unknowns.  The  actual  solution  of  a  truncated  version  of  such  a  system  will  be 
given  in  Part  3  of  this  report.  We  conclude  by  deriving  expressions  for  the  far  fields  in  the  half¬ 
space  in  which  the  coaxial  line  radiates. 
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1.  INTRODUCTION 


We  use  here  the  integral  equations  developed  in  Part  1  to  determine  the  fields  of  the  structure 
in  Figure  1.1:  a  coaxial  transmission  line,  tenninating  in  a  ground  plane  and  radiating  over  the 
upper-half  space.  All  surfaces  are  perfectly  conducting.  The  ground  plane  is  infinite  and  the 
coaxial  line  extends  to  infinity  in  the  lower-half  space.  The  inner  and  outer  radii  of  the 
transmission  line  are  a  and  b,  respectively.  We  assume  that  the  line  supports  only  a  TEM  wave. 
For  this,  it  is  sufficient  that  (references  1  and  2) 


a  +  b 
2 


<1 


(1.1) 


where  k  is  the  wavenumber  of  the  time -harmonic  ( e+,(0‘ )  electromagnetic  wave  in  the  line. 

4  ^ 


(7 


a 


Figure  1.1.  The  semi-infinite  coaxial  line  of  inner  radius  a  and  outer  radius  b  is  fed  at 
z  =  — oo  .  At  z  =  0 ,  it  opens  up  into  a  half-space,  its  outer  conductor  becoming  a 
plane  that  extends  to  infinity.  All  surfaces  are  perfectly  conducting. 


The  excitation  of  the  line  occurs  at  z  =  -oo  and  results  in  a  TEM  wave  with  fields 


Eg(r)  = 


V  e 


In 


\a) 


-ikz  YV  p-ikz 

P,  H*(r)  =  — —  <p. 


P 


In 


Uv 


p 


(1.2) 


where  V  is  the  voltage  of  the  inner  conductor  with  respect  to  the  outer,  and  Y  is  the  free-space 
admittance.  We  use  cylindrical  coordinates  (p,(p,z).  Since  the  line  is  not  infinite,  we  have  also 

an  induced  wave  in  the  line  with  fields  { E'  ,H'  }.  The  total  fields  { E'  ,H' }  in  the  coaxial  line  are 
the  sum  of  the  generator  and  induced  fields 

E'  =  Eg  +  E'  ,  H=Hg+H.  (1.3) 
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In  the  upper-half  space  we  have  radiated  fields  { E'  ,Hr }.  We  proceed  to  determine  both  sets  of 
fields  by  first  deriving  integral  equations  and  then  solving  them  numerically. 
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2.  FIELDS  INDUCED  BY  COAXIAL  LINE 


Since  the  generator  fields  are  independent  of  the  angular  coordinate,  and  since  the  geometry 
has  circular  symmetry,  we  expect  the  induced  fields  to  be  independent  of  (p.  Thus,  according  to 
Appendix  A,  we  can  have  a  TEM  wave  traveling  down  the  line,  given  by 


E'(r)  = 


V  e 


+ikz 


In 


f 


P,  H'(r)  = 


P 


YVe 


+ikz 


<P- 


In 


vU 


P 


(2.1) 


In  addition,  we  can  have  the  zeroth-order  mode  of  a  TM  wave,  given  by  (A3. 23),  (A3. 24),  and 
(A3. 27)  in  Appendix  A 


m  =  p„n(p) 

(2.2) 

m  EUp,<P,z)  =  ihJ’---  pJAp) 

(2.3) 

™HUp,v,z)  =  ikYeA-Pj(p) 

(2.4) 

with  v0n  a  root  of  the  equation 

■/o(h„fl)7oM)-'/o(’/o,i)7oK,fl)  =  0  »  n  =  1,2,... 


and 


jk2 -v0„2  , 

k 1  >  Von" 

r^V0n2  ~k2  3 

•  V0n  >  ^ 

while 


(2.5) 


(2.6) 


P^P)^J^V^nP)- 


•h (V) 

E0(v0;,a) 


Yo(vonP) . 


(2.7) 


Similarly,  we  can  have  the  zeroth-order  mode  of  the  TE  wave  (see  Appendix  A) 


HI ne  (A  A  *)  =  Pon  ^  R0 n  (P) 

(2.8) 

HL(p’<p>z)  =  i/c0n  qIk°":  rJ(p) 

(2.9) 

EL  (A  <P,  z)  =  ikZ  eilc°nZ  Rj  (p) 

(2.10) 
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with  ju0n  a  root  of  the  equation 


J<i (Pona)Yo  (Ponb) ~J')  (Ponb)Yo (Po„a)  =  0  ,  n  =  1,2,... 


and 


while 


(2.11) 


(2.12) 


KXp)  =  J^Po„P) 


■A,  (Ao igfO 

YUPdna) 


Yo(Po„P)- 


(2.13) 


The  first  few  roots  of  (2.5)  and  (2.1 1)  are  given  in  (reference  3)  for  selected  values  of  %  =  b/ a  . 
Also,  Mathematica®  (reference  4)  has  special  subroutines  that  compute  the  roots  of  these 
equations. 

From  the  above,  we  can  write  that  the  magnetic  field  induced  in  the  coaxial  line  is 


H'(r)  =  -A 


YVe 


+ikz 


(p 


In 


P 


oo  / 

+  V^lBJkY^  Pj (p)<p  +  C„  eiK°"z 


<nKn  (p)P  +  Pon  R0n(p)2 


(2.14) 


where  A,  Bn ,  and  Cn  are  unknown  constants.  The  constant  A  is  dimensionless,  Bn  has  dimension 
of  length,  and  Cn  has  dimension  of  length-mhos. 
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3.  INTEGRAL  EQUATIONS  AND  ELECTRIC  CURRENT  DENSITIES 

From  (5.6),  (7.1),  and  (7.2)  of  Chapter  1  of  Part  1,  the  integral  equations  for  the  present 
problem  are 


z  x  f>  a  (r)  x  Vg(r,r')dS  +  J  J,  (r)  x  Vg(r,  r' ')dS  =  JCT  (r')  ,  r'  e  <r 


i  x  J  f Ja  (r)  x  Vg(r,  r')  -  [ J0  (r)  x  Vg(r,  r/)]  }  dS 

sa 

~n' x  { { J,  (r)  x  Vg(r,r')  -  [Jfi  (r)  x  Vg(r,r/)]. }  dS 


-2/7'xjjff(r)xVg(r,r')^  = 


2  J*(r'),  r'eSb 


The  current  densities  in  these  equations  are  defined  in  (4.6),  (4.7),  and  (4.8)  of  Chapter  1  of  Part 
1. 

Jfl(r)  =  HxH'(r),re5a,  J4(r)  =  nxH'(r)  ,  r  e  Sb  ,  JCT(r)  =  -zx  H'(r)  ,  r  e  a  (3.3) 

where  H'  is  defined  in  (1.3). 

In  terms  of  the  fields  in  (1.2)  and  (2.14),  we  have 


JCT  (r)  =  (i  -  A)  p  +  iV^lBnkYPj  (p)p  -  Cnrc0nR0n'  (p)q> 

In  I  -\p  n='K 
\a 


while 


Ja < 00  =  Td\  ( -  A e+*  )z  +  VU  B*ikY ^  Po,!  ia)z  -  Cn  e""s  p0n2R0n  (a)cp 

In  L  n=1 


J,  (r)  = - ( e-ife  -  A  e+'fc )  z  -  u£  { BJkY  ^  Pj  (b)z  -  Cn  e^z  p0„2R0ll  {b)q> }  (3.6) 

In  —  b  "=1 

l  a 


52 


NAWCADPAX/TR-2013/1 15 


But  from  (2.7) 


P»n  0)  =  V, 


On 


j o'  ( v0n fl)  -  y0(l/o"a)  Yd  (Yo na) 

Yo(Yo  na) 


Vn 


Yo(VOna) 


(y*na)Yo(yQna)-Jo(y0na)YV  (y0na) 


xaY0(v0na) 


(3.7) 


where,  above,  we  made  use  of  the  Wronskian  relationship  for  Bessel  functions  (reference  3). 
Similarly,  from  (2.13) 


,»  = 


1 

^>'(Ao„«)L 


J0  (Porfi)  Y<)  (Ao  na)  J 0  (Ao/A)^0  (Ao/A) 


^MonaY0  (Ao /A) 


(3.8) 


With  these  two  results,  we  can  write  for  (3.5) 


K (r )  =  -j^r  ( e*  -  A  )  f  ■ - —  £ 

aln(j)v 


ikYBn  _  Ao»Q  e/%,z  ^ 

^(WA)  ^(AoA) 


(3.9) 


where  ^  =  b  /  a  . 

From  (3.15)  and  (4.10)  of  (reference  1),  we  can  write  (2.7)  and  (2.13)  in  the  form 


pQ,XP)  =  J^nP)-  Yo(yonP)’  *o„(AWo (  Ao„ A)  “  ^  ^0 ( Ao„ A) 

7«(Vo^)  roW) 


(3.10) 


Proceeding  as  above,  we  find  that 


n>>  =  -  ,.,2  ...  *„„(0  =  - 


xbY0(v0nb) 


xp0nbYXP0nb) 


(3.11) 


and 


J  b<?)  =  -  777~XT  ( ^  ^  e+<fe )  z  +  27  J  |  e^  z  -  e-*  p  J. . 


In  (z)b 


7Tb  tt  \Y0(y0nb) 


YXPo„b) 


(3.12) 


In  the  next  section,  we  substitute  the  current  densities  in  the  integral  equations  and  resolve  the 
latter  into  their  components. 
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4.  SCALAR  FORM  OF  INTEGRAL  EQUATIONS 

We  substitute  next  expressions  (3.4),  (3.9),  and  (3.12)  for  the  electric  current  densities  in  the 
integral  equations  (3.1)  and  (3.2).  First,  we  write 

J<T(r)  =  ua(r)p  +  va(r)v ,  J0(r)  =  va(r)tp  +  ha(r)z  ,  JA(r)  =  vb(r)<p  +  hb(r)z  (4.1) 

where 


YV 


Ua  (r)  =  (1  -  A)  +  ikYV Z  Bn  P0n  (P)  ,  Va  (r)  =  ~^Z  CnKOnRo!  (P) 


n- 1 


n- 1 


v  (r)  =  Mo"C "  e'*--', 

na^Yx^0na) 


a^{x)  Ka 


L>„z 


(4.2) 

(4.3) 


2V  tA  iuC  z  ,  FL  /  ,fe  ,  +/faN  likYV^  B, 

\  ^o»  »  e'*b„z  A  (r)  = - —  e  -Ae+,i-  + - >  - " — e 

"  1  ’  nb  hUv0„b) 


u(r)  =  -^,Zw  i\ 
xb  „=1  F,(/vA) 


Aln(j) 


'A,2 


(4.4) 


Compute 

f  x  [  (r)  x  Vg(r,  r’)]  =  J„  (r)  ALL  -  Vg(r,  r 1 ')z  ■  J„  (r) 

OZ 

=  [va  (r)^  +  ha  (r )z]  ~  V§(r’  r ')z  ■  [va  (r)p  +  ha  (r )z] 

=  vfl (r)  V~K (r)V,g(r, r ')  (4.5) 

oz 


where  the  subscript  t  stands  for  “transverse  to  the  z-axis”.  A  similar  expression  exists  for  the 
outer  conductor.  Substitute  in  (3.1) 


dg(r,  r') 

dz 


<P 


^(r)V(g(r>r') 


</S+ J 


vb  (r)  r  V  -  K  (r)Vfg(r,rf) 
oz 


dS 


=  u*(r')p'  +  v{r(r')<p'  ,r'  e  cr  . 


(4.6) 


This  is  the  first  integral  equation  in  terms  of  scalar  functions. 

Compute 

J  8  (r)  x  Vg(r,  r ')  -  [  JB  (r)  x  Vg(r,  r/)] 

=  [va  (r  )q>  +  ha  (r  )z]  x  Vg(r,  r')  -  [(va  (r)£  +  ha  (r  )z)  x  Vg(r,r/)]. 

=  va  (r)  {(p  x  Vg(r,  r ')  -  [^  x  Vg(r,  r/)] . }  +  ha  (r)z  x  V  (g(r,  r ')  -  g(r,  r/)) 
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=  v8(r) 


=  va(r) 


^dg(rj')  „dg(r,r')  ^dg{ r,r.  )  ^dg( r,r.  ) 


d  z 


dp 


dgD(r,r’)  „  dgN (r , r') 
dz  dp 


dz  dp 

+  ha(r)zxVgD(r,r') 


+  ha  (r)z  x  V  (g(r,  r ')  -  g(r,  r/)) 

(4.7) 


where 


&>(*■>  O  =  g(r,r')  -g(r,r/)  ,  gw(r,r')  =  g(r,r')  +  g(r,r/) 


(4.8) 


the  Dirichlet  and  Neumann  Green’s  functions  for  the  z  -  0  plane.  Equation  (4.7)  is  also  valid 
when  the  subscript  a  is  replaced  by  b.  We  also  compute 


JCT  (r)  x  Vg(r,  r ')  =  [ua  (r  )p  +  (r)0]  x  Vg(  r,  r ') 

=  [pva  (r)  -  (pua  (r)]^r’r)  +  z 

dz 


o  v  /  cr  v  / 

p  dcp  dp 


(4.9) 


Substitute  these  expressions  in  the  first  of  (3.2)  (r '  e  Sa  ,  n  =  p  ) 


~P'X  hW 


-P'x  fiv4(r) 


dgD(r,r')  ^dgN{ r,r’) 
dz  dp 

dgD(r,r')  „dgN(  r,r') 
P - z 


dz 


dp 


+  ha(r)zxVgD(r,r')}dS 


+  hb(r)zxVgD(r,r')\dS 


-2 p'  x  | |[pv(T  (r)  -  (pua  (r)]  V  }  +  z 

=  ]pJr'‘IV  -  hjr')t\,  r'ej,. 


u  ( r)l^g(r,r')_ 

p  d(p  dp 


>dS 


(4.10) 


Since 


p  =  cos cpx  +  sin (p y  ,  (p  =  -sirupx  +  cos(py 


(4.11) 


we  can  perform  the  vector  operations  in  (4.10)  to  get,  for  the  cp'  -component, 


dgN  (GO 


l  dp 


dS-  Jv6(r) 


+  2\ 


dgN(r,  r') 

dp 

dg(  r,r')l 


dS 


1  <7g(r,r ) 

ua(r) - ' - vCT(r) 

p  d(p  dp 


1 


dS  =  -va(r’),  r'eSa 


(4.12) 
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and  for  the  z-component 


- J  (r) sin(^  -cp')  dgD^V  )  +  K (r )p'  ■  VgD (r, r ')| dS 
- J  jv4 (r) sin(^  -<P')  -  +  hb(r)p'-  v§d (r,i r ')| dS 


2j  [sin(p  -  <p')va  (r)  -  cos(^  -  (p')ua  (r)] 


^(^)-dS  =  Ua(  r'),  r'eSfl 


Similar  expression  is  obtained  from  the  second  of  (3.2)  ( r'  e  Sb  ,  n  =  -p' ) 


i  dp  i  dp 
-  2|L(r)IAMt  _v,(r) ^Pl]dS  .  '(O  ,  r-  6  S„ 

J I  p  oq>  op  J  2 


dgN(  r,r') 


va  (r)  sin((p  -  <p')  )  +  K  ir)P'  •  V gD  (r,  r ')  j dS 


+ J  jv4(r)sin(p  -  <p')  -  +  hb  (r)p'  •  VgD(r,r')|^ 


+  2j [sin(^>  -  (p')va  (r)  -  cos(^  -  (p')ua  (r)] 


^pPdsP-ht(r'),  r'  e  St 


We  return  to  (4.6)  and  we  write  it  in  terms  of  its  components 


-j  vfl  (r) dg{y’ r  ^  sin(ff>  -  cp')  +  ha  (r  )p'  •  V,g(r,  r ')  dS 

dz  J 

-f  vb (r) d8^' r  ^ sinO ~ <P')  +  K ir)P' •  vr<g(r> r ')  dS  =  ua(r'),r'  ecr 
i  oz 


(4.13) 


(4-14) 


(4-15) 


(4-16) 
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S„  L 


Va  (r) dg(y,X  1  cos i<p  -  <p')  -  ha  (r )<p'  •  V,g(r,  r ') 

oz 


dS 


S,<- 


Vb  (r)  ^r,r  -  cos(^  -  <p')  -  hb  (r ■  V,g( r,  r ') 

oz 


dS  =  va(r') ,  r'  g  cr 


The  system  of  integral  equations  is  given  by  (4. 12)  through  (4.17). 


(4.17) 
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5.  SIMPLIFICATION  OF  SCALAR  INTEGRAL  EQUATIONS 
We  recall  the  definition  of  the  scalar  Green’s  function 

-ikR  i 

g(r,r')  =  -— — ,  R  =  |r - r'|  =  [p2  -  2  pp’ cos((p  -  (p')  +  p’2  +(z-z')2]2 . 

4  nR  L  J 

We  observe  that 

f.Wg(r,r  ",  =  _  pWg(  r.r1)  =  „ 

•U  dcp'  V  d(p 


(5.1) 


(5.2) 


Since  the  various  current  density  functions  are  independent  of  the  second  term  in  the  two 
integrals  in  (4.17)  is  zero;  thus,  we  can  write 

f  Va (r)  dg(r,r  1  cos(^  - +  f  vb (r)  1  Cos((p  -<p')dS  =  va( r ')  ,  r '  e  a  .  (5.3) 

•  dz  t  dz 

This  statement  makes  good  physical  sense  because  it  says  that  the  angular  linear  current  density 
on  cr  depends  only  on  the  angular  linear  current  densities  on  the  walls  of  the  coaxial  line  and  not 
on  the  z-directed  ones. 


We  encounter  the  same  kind  of  term  in  (4.12),  and  we  rewrite  that  expression  as 


J  v-(r> 


dgN(r,r') 

dp 


dS-\vh(r) 


dgN(r,r') 

dp 


dS-  2jv„(r) 


=  U(r'),  r'ES,. 
dp  2 


(5.4) 


A  similar  expression  results  from  (4.14) 


f  (r’ r  )  jo  ,  f  r^dgN(r,  r)  r  <?g(r,r)  1 

I  va(r) - - - +  I  v6(r) - - - ^  +  2  va(r) - - - dS  =  -vb(  r),r  g  . 

s  dP  S,  dP  i  dp  2 


(5.5) 


From  (5.1),  we  see  that  the  scalar  Green’s  function  in  cylindrical  coordinates  is  an  even 
function  of  cp-  cp'  about  zero.  Thus,  its  Fourier  series  contains  only  cosine  terms;  in  fact 
(reference  3,  p.  487) 


g(r,r')  =  --1-XGcos/(C-^')fn  J, (ap)J,(apr) e  ad°  (5.6) 

4 Jo  4a2 -k2 

where  s,  is  the  Neumann  symbol 
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£0=l;  e,= 2 ,  /  =  1,2,... 


(5.7) 


Because  of  this,  the  terms  in  (4.16)  containing  the  sine  function  dropout  and  we  get 


dg(r,r') 

dp' 


dS+\hb{  r) 


sb 


dg(r,r') 

dp' 


dS  =  ua(r')  ,  r'  e  a  . 


(5.8) 


Again,  this  statement  makes  good  physical  sense  because  it  says  that  the  radial  component  of  the 
linear  current  density  on  cr  depends  only  on  the  z-component  of  the  linear  current  density  on  the 
walls  of  the  coaxial  line.  We  can  visualize  the  current  coming  up  in  the  inner  conductor, 
displacing  itself  radially  at  the  termination  of  the  line,  and  returning  via  the  outer  conductor. 
Similarly,  from  (4.13)  we  write 


j  K  M  ^°(l/r')  dS  +  J  hb (r)  ,'s°vV  '  dS 

{  dd  l  dd 

+  2j  ua  (r)cos(^  -  cp') 


dgD<jy) 


dz  2 


and  from  (4.15) 


(5.9) 


\K(*)^PpP<lS-\hb(  r) 


pgo^y) 


dp’  l  dp' 

~  2 J iia  (r)  cos(^>  -<p') 


dS 


,yg(±y)ds=  iA>(r,)jr.eS 

dz  2 


(5.10) 


From  definitions  (4.2),  (4.3),  and  (4.4),  we  see  that  the  last  three  equations  involve  both 
unknown  and  known  quantities.  By  contrast,  (5.3),  (5.4),  and  (5.5)  involve  only  unknown 
quantities,  i.e.,  there  is  no  generator  that  excites  these  circumferential  currents.  Indeed,  a  quick 
examination  of  the  geometry  of  our  problem  and  the  way  we  excite  the  coaxial  line  reveal  that 
there  is  no  reason  for  these  currents  to  exist;  therefore,  in  (4.2).  (4.3),  and  (4.4),  we  set  all  Cn  =  0 
.  From  (2.14),  and  the  discussion  in  Section  2,  we  conclude  that  the  arrangement  we  have 
supports  a  TEM  and  a  TM  wave,  but  it  does  not  support  a  TE  wave. 
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6.  SOLUTION  OF  THE  SYSTEM  OF  INTEGRAL  EQUATIONS 

The  system  of  integral  equations  consists  of  (5.8),  (5.9),  and  (5.10).  The  three  scalar  functions 
that  appear  in  them  are  defined  in  (4.2),  (4.3),  and  (4.4)  and  involve  an  infinite  number  of 
unknowns.  We  can  use  the  method  of  weighted  residuals  (reference  6)  to  convert  each  equation 
into  an  infinite  system  of  linear  algebraic  equations.  Thus,  we  will  have  three  separate  systems  of 
equations,  each  involving  the  same  unknowns.  If  we  truncate  the  three  systems  into  finite 
systems  of  N  + 1  equations  in  N  + 1  unknowns,  then,  as  N  tends  to  infinity,  the  solution  for  each 
system  must  converge  to  the  same  value;  thus,  we  only  need  to  numerically  solve  one  of  the 
three  integral  equations.  We  chose  to  deal  with  (5.8).  We  will  do  this  in  Part  3  of  this  report.  In 
the  next  section,  we  derive  far-field  expressions. 
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7.  FAR-FIELD  REPRESENTATIONS 

We  derive  a  far- Held  expression  for  the  geometry  of  this  section.  From  Green’s  second 
identity  (reference  7,  p.  509),  we  have  that  in  the  upper-half  space 

J  {ir  (r)  •  [ V  X  V  X  r2  (r,  r ')]  -  [ V  x  V  x  Hr  (r)]  •  T2  (r,  r ')}  dV 

D+ 

=  J  jlT (r) •  [/) x V x T2 (r, r ')] -[«xVxHr (r)] •  T2 (r, r ')} dS  (7.1) 

z=0 


where  H'  denotes  the  radiated  magnetic  field  (total  magnetic  field  in  upper-half  space)  and  T2 
is  the  second-kind  Green’s  dyadic.  The  latter  satisfies  the  differential  equation 

V  x  V  x  r2  (r,r')  -  k2T2  (r,r')  =ikV  x  [<7(r,  r')  I  -  S( r,r/)I,.  ]  (7.2) 

where  I  is  the  identity  dyadic  and  I,  its  image  about  the  z  =  0 ,  and  the  boundary  condition 


zxVxR  (r,r')  =  0,  z  =  0. 

Following  standard  procedures,  we  can  show  that 

J|H  '  (r)  •  [ V  x  V  x  r2  (r,  r ')]  -  [’ V  x  V  x  H  (r)]  •  T2  (r,  r ')}  dV 

D+ 

=  ikV'x  H'  (r')  =  -k2YEr  (r') ,  z'  >  0 . 

Also, 

J  \W (r)  •  [n  x  V  x  r2 (r, r ')]  -  ikY[n  x  W (r)]  •  T2 (r, r ')}  dS 

z= 0 

=  -ikY\ {[z  x  Er  (r)]  •  r2  (r,  r')}  dS  =  -ikY$ Ma  (r)  •  I\  (r,  r ')dS . 


From  the  last  two  equations,  we  have  that 
E'(r')  =  f  fMCT(r)-r2(r,rVS,  z' >  0  . 

K 


<7 


The  dyadic  in  this  expression  is  given  by 
r2(r,r')  =  -ikV  x  [g(r,  r ')  I  -  g(r ,  r')  T;  ] 


(7.3) 


(7.4) 


(7.5) 


(7.6) 


(7.7) 
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Where 


g(r,r')  =  - 


-/A'lr-rl 


An  r-rl 


(7.8) 


With  this,  we  have  that 


z=0 


(r)  x  Vg( r,  r ')  -  [Ma  (r)  x  Vg( r,  r/)]  •  I, } 

=  Ma  (r)  x  V,g(r,r')  -  [Mct  (r)  x  V,g(r,r')]  •  I,.  +  Ma  (r)  x  z  )  + 

/pa(r  r h 

=  2M„ (r) x V,g(r,r')  +  2M„ (r) x 2  =  2M„ (r) x Vg(r.r')  „„ 

OZ 


MCT(r)xz 


^g(fA') 

d  z 


I, 

(7.9) 


which  we  substitute  in  (7.6)  to  get 

Er  (r ')  =  2j  Mct  (r)  x  Vg(r,  r ')dS  ,  z'  >  0  . 


(7.10) 


Similarly,  from  (2.13)  of  Chapter  1  of  Part  1,  we  have 


H''  (r ')  =  W7  |X  (r)  •  Wg(r, r')dS  -  i2kY  f  MCT  (r)g(r, r ')dS  ,  z'  >  0  . 
ikZ J  J 


(7.11) 


This  equation  can  also  be  obtained  by  applying  one  of  Maxwell’s  equations  (Ampere’s  law)  to 
(7.10). 


In  the  far  field  (  r'  — >  oo  ) 

-/A(r'-r-r') 


-/A(r'-r-r') 


g(r,r') - —a — - —  ,  Vg(r,r')  ~  ik  — — ; — r' ,  VVg(r,r')  ~  k2  - 


4  nr' 


Anr' 


-ik(r'-rr') 

Anr’ 


'V'V 

-r  r  . 


(7.12) 


Using  this  in  the  previous  two  equations,  we  get 

ik  Qikr' 


E'(r')  =  fMa(r) eikr'f'  dS  ,  r' 

2  nr  J 


— >  oo 


(7.13) 


and 


H'(0  = 


ikY  e 


-ikr 


2nr' 


-r  x 


x  J  MCT(r) 


Qikr  'V  dS 


r  — >  oo . 


(7.14) 


These  are  the  far-field  representations. 
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We  will  attempt  to  compute  the  integral  in  these  expressions.  From  (1.2)  and  (2.14) 


Mff  (r)  =  Mff  (p, cp,0)  =  zxEr (p, <p, 0)  =  z  x E' (p, cp, 0) 
1  +  A 


=  V 


- 1 

HZ)P  n= 1 


(p> 


<P- 


(7.15) 


We  substitute  in  the  integral 


I  =  J  Mct  (r)  ekr'y  dS  =  f£  pd p 


1  +  A  ^  n  „  „  ,,  N 

,  7  ,  ip) 

In  (x)P  »=i 


J 


(7.16) 


For  the  second  integral,  we  write 


p2;r  ^  »,  p2;r 

I,  =  d<p<pe,ir'r  =  J  (-xsin^>  +  ycoscp)e 

p  1.71— (p'  -  ^ 

=  J  [-xsin(^  +  ^')+ vcos(^  +  ^')Je 

=  -p'fWsin^eftsin'9>cosr 

J —(o'  J—cp' 


ik  sin  &'pcos{(p-(p') 


dcp 


ik  sin  «9'/?  cos  ip 


dy / 


eos^e'isin,9>cosr  dy/  =  p'i2nJx(ksm9'p)  (7.17) 


and  substitute  in  (7.16) 


I  =  cp'ilnV  [  J,  (A:  sin  $'/?) 

J  a 


1  +  A 

- 1 

HZ)P  n= 1 


GO 

YjBn^nPo!(p) 


pd  p 


1  “f"  ^  b  ^  ^ 

=  cp'ilnV  \  — —  f  A,  (A:  sin S' p)d  p  -  iJ^BJ 0„  f  J,  (A:  sin  9'p)Pj  (p)pdp  [ . 

[Hz)Ja  t!  Ja  J 


(7.18) 


For  the  first  of  the  last  two  integrals,  we  write 


rb  \  r  kb  sm  S'  f 

7i  =  |  J ! (k sin 9' p)d p  =  —  ,  I  .  J0(u)du  = 

Ja  kSmd  J^sm5 


J0  {ha  sin  9')  -  J0  (kb  sin  9') 
k  sin  9' 


(7.19) 


while  for  the  second 


I 2  =  fV)(A:sin 9'p)Pj (p)pdp  = -v0n  \b  J^ksimS'p) 

J  a  Ja 


A  (YonP)~ 


(VQna) 

Yo(vona) 


Yl  (Vo„P) 


pd  p 


V0nP 


v0n2  -  k 2  sin2  9 


-  lv0nJ2  Oo nPVx  ik  sin  S' p)  -  k  sin  9'Jl  (v0np)j2  ( k  sin  S' p)\ 


VpnP 

Yo(vo„a)  v0n2  -k2  sin2  S 


,W,Yi  (y0nP)J i(k  sin  S' p)  -  kYx  (i v0np)J2  (k  sin  S' p)\ 
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R>„  P 


v02  -  k 2  sin2  9' 
v0nksin9'p 
v0n2-k2  sin2  9' 


J 2  (Yon  P)  -  Y2  (V0n  P) 

Y0(v0„a) 


Jx  (k  sin  9'p) 


A(y0  „p)~  Jy^nCl\  Yi  ( vo„  P) 

Yo(VOna) 


J2  ( k  sin  9'p) 


2v. 


On 


Vq2  -  k1  sin2  9' 
v0nk  sinS'p 


•AOo  nP)~ 


■/p  (V()na) 

Y0(v0na) 


YMo  nP) 


Jx  ( k  sin  9'p ) 


b 


v0  2  -  k 2  sin  9 


1 2 


A  (vo„P)- 


J0  (R) na) 


v0nk  sin  9'p 
2  -k1  sin2  S' 
v0nk  sin  9'p 


WO  nH>  v/  .YMonP)  \J2(k$m9'p) 

Yo  (vona) 

Jx  (Yo  nP)  -Jv0(.V°"alYx(^nP) 

Y0  (V0na) 


2J\  (k  sin  9  p) 

.  n, - J2(ksvn9p ) 

k  sin  9  p 


vQ2  -  k 2  sin2  9' 
v0nkb  sin  9' 


J 1  Oo n  P)  ~  Y1  (Vfl„  p) 


Yo(^,P) 


v0  2  -k2  sin2  9' 
v(]Hka  sin  9' 
v0lI2  -k2  sin2  9' 
v0nkb  sin  9' 
v02-k2  sin2  9' 
v()nka  sin  9' 


Ji<y0nb)~ 


Jq  (Aq«A) 

Yo(vo„a) 


Yl(V0nh) 


Joi^Oip) 


v02  -  k 1  sin2  9' 


Jx(v0na)-^-^Uv0na) 

Yo(VOna) 

Yo(y0nb) 

Jx(y o„a)  -  Yx(y0na) 


J0  (k  sin  9'p) 
J0(kb  sin  9') 
J0  (ka  sin  9') 
J0(kb  sin  9') 
J0  (ka  sin  $') 


Yo(y0,p) 

2k  sin  9'  J0(kb  sin  9')  2k  sin  9'  J0(ka  sin  9') 


-  +  - 


n(v0;  ~  k2  sin2  *9')  Y0  ( v0„b )  7i(v()n  -  k 2  sin2  9')  Y0  ( v0na ) 


2k  sin  9' 


7t(v,2  -  k2  sin2  9’) 


J0(ka  sin  $')  J0(kb  sin  9') 
Y0(v0na)  Y0(v0nb) 


We  substitute  the  last  two  equations  in  (7.18)  and  use  (2.6) 


T  tA\  +  A  J0(ka  sin  9')  -  J0 (kb sin 9') 

1  =  cp  i27rV  ( - 

in(x)  ksin  9' 


2A:sin  9' 

Vn~  2 
On 


«  -  kz  sin2  9’) 


J0(ka  sin  9') 

J0  (kb  sin  9')  } 

Y0(v0na) 

-  J0  (kb  sin  9') 

Y0(v0„b)  Jj 

Un(^) 


£sini9' 


(7.20) 


64 


NAWCADPAX/TR-2013/1 15 


2k  sin  3'  ~k 

+ - LBr  - - - 

«= 1  yO, 


n  „=i  y2  -k1  sin7  3’ 


J0  (ka  sin  3') 

J0(kb  sin  3') 

Y0(v0na) 

Y0(v0nb)  \ 

From  this  and  (7.13) 


-ikr 


Er(r)  =  -5- _ Vk  \l  +  AJ^ka  sin  ,9)  ~  jq  (kh  sin  ^ 

r  k  sini9 


Vv2-^2 


2k  sin  3  ^ 

+ - >  B , 

n  B=1  v0n  -k  sin  3 


J0  ( ka  sin  3)  J0  ( kb  sin  3) 

Yo  (y0na)  Y0(v0nb) 


> ,  r  — >  oo 


while  from  (7.14) 


n  /r  )  _  >  e~f _  Yyk  l  +  A  J0(ka sin  3)  -  J0(kb sin  3) 
r  M/f)  ksini9 


z  J  2. 

y0„  ~k 


2k  sin  3  -A  n 

+ - 2X  2  ,  2  •  2  Q 

v0n  —  a:  sin  ,9 


«=1 


J0  (ka  sin  3) 

Jq  (kb  sin  3) 

Y0(v0na) 

Y0(v0nb)  \ 

> ,  r  — >  oo 


(7.21) 


(7.22) 


(7.23) 


where,  in  these  two  expressions,  we  changed  from  primed  to  unprimed  coordinates.  These  are 
the  far- Held  representations  of  the  electric  and  the  magnetic  field  in  terms  of  the  unknown 
coefficients  of  the  current  density  expansions. 
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8.  CONCLUSIONS 

For  a  coaxial  line,  as  shown  in  Figure  1.1,  we  derived  three  BIEs.  We  also  wrote  the  unknown 
current  densities  as  expansions  in  the  modes  of  the  coaxial  line.  We  concluded  in  Section  6  that 
we  need  deal  with  only  one  of  the  integral  equations  to  numerically  determine  the  coefficients  of 
these  expansions.  This  will  be  explicitly  done  in  the  parts  that  follow.  In  Section  7,  we  present 
the  expansions  of  the  far  fields. 
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APPENDIX:  THE  ELECTROMAGNETIC  FIELDS  OF  A  COAXIAL  TRANSMISSION  LINE 
Al.  INTRODUCTION 

In  this  note,  we  develop  the  eigenvalues  and  eigenfunctions  for  an  infinite  coaxial  line.  The 
inner  and  outer  walls  of  the  line  are  perfectly  conducting,  with  radii  a  and  b,  respectively  (see 
Figure  Al).  The  axis  of  the  line  is  the  z-axis.  Derivations  can  be  found  in  many  books.  We  have 
consulted  Jones  (reference  8),  Stratton  (reference  9),  and  Tai  (reference  10).  The  time 
dependence  of  the  electromagnetic  fields  is  +icot.  Thus,  Maxwell’s  equations  are 

V  x E  =  -ikZR  ,  V  x  H  =  ikYE  ,  k  =  ,  Z  =  Y  1  =  ^/u/ s  (Al .  1) 

where  //  and  e  are  the  constitutive  parameters  of  the  space  between  the  two  conductors.  On  Ca 
and  Q,,  the  tangential  component  of  the  electric  field  is  zero 

nxE  =  0.  (A  1.2) 


y 


Figure  Al.  Cross  section  of  coaxial  transmission  line.  Ca  and  Q,  represent  the  inner  and  outer 
circles.  The  region  between  the  two  conductors  is  denoted  by  D. 

A2.  GENERATION  OF  WAVES  IN  THE  COAXIAL  LINE 

It  is  well  known  (reference  9)  that,  in  a  source-free  region,  electromagnetic  fields  can  be 
expressed  in  terms  of  two  scalar  functions.  These  functions  are  components  of  Hertz  vectors  of 
arbitrary  direction.  In  this  case,  we  take  the  direction  to  be  that  of  the  z-axis  and  write 

II(r)  =  77(r)z  ,  M(r)  =  M(r)z  (A2.1) 
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where 

V277(r)  +  k2n(r)  =  0,  V2M(r)  +  k2M(r)  =  0  .  (A2.2) 

The  fields  generated  by  these  Hertzian  potentials  are  (reference  9) 

E(r)  =  VxVx[77(r)z]-;HVx[M(r)z],  H(r)=VxVx[M(r)z]  +  ;17Vx[i7(r)z].  (A2.3) 


A3.  TRANSVERSE  MAGNETIC  WAVES 

If  in  (A2.1)  we  let  M(r )  =  0 ,  we  find  that  the  magnetic  field  does  not  have  a  component 
along  the  z-axis.  The  z-component  of  the  electric  field  is  given  by 

™  Ez (r)  =  z  ■  V  x  V  x  [77(r)z]  =  +  k2ll(r)  .  (A3. 1) 

dz 

From  (A  1.2),  we  have  that  ™  Ez  must  be  equal  to  zero  when  p  =  a  or  p  -  b  .  For  this,  it  is 
sufficient  that 


n  (a,(p,z)  =  n(b,(p,z)  =  0 .  (A3. 2) 

For,  if  this  is  the  case,  it  follows  from  the  Hugoniot-Hadamard  theorem  (reference  11)  that  the 
second  derivative  with  respect  to  z  is  also  zero  there.  Thus,  from  (A2.2)  and  (A3. 2),  we  have  the 
boundary-value  problem 

V277(r)  +  k2  IT(r)  =  0 ,  n{a,cp,z)  =  n{b,cp,z)  =  0 .  (A3.3) 

We  use  cylindrical  coordinates  ( p,(p,z )  and  write 


n{p,(p,z)  =  f(p,(p)h(z). 

Substitution  in  (A3. 3)  gives 

v;f{p,(p)  _h\z)  [  fc2_r2 
f(p,<P )  h{z) 


(A3. 4) 


(A3. 5) 


where  V,  is  the  transverse  component  of  the  Laplacian  and  v2  is  a  constant  that  will  be 
determined  from  the  boundary  condition.  Expressing  the  Laplacian  in  cylindrical  coordinates,  we 
have 


P 


d 

dp 


P 


df_ 

dp 


+  ^-  +  v/V/(AC)  =  0. 
d(p~ 


(A3. 6) 
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Writing, 

f{p,cp)  =  P{p)0{(p)  (A3. 7) 

we  obtain 
d 

P  dp 

where  the  constant  m  is  an  integer  since  the  solution  in  the  angular  direction  has  period  2  n.  The 
function  Q  is  then  of  the  form 


.  dP(p) 


P(p) 


2  2 

+  V p  =  • 


&”(<P)  2 

- =  m 

0{(p) 


(A3. 8) 


f  cos  mm 

<*>»  =  .  ,  m  =  0,1,2,.... 

m  \s\nm(p 

The  remaining  equation  in  (A3. 8)  is  Bessel’s  equation 

pd  pdP(p)  +^vip2_m2  J/>(p)  =  0 

d  p  y  d  p  J 

and  its  solution  is  the  linear  combination  of  a  Bessel  and  a  Neumann  function 
p,»  (P)  =  amJm  (vmp)  +  bmYm  (vmp) . 

These  are  the  functions  that  also  must  satisfy  the  boundary  conditions 

PnM)  =  pm(b)  =  o. 

This  results  in  the  homogeneous  system  of  equations 

amJm{vma)  +  bmYm{vma)  =  Q 
amJm(ymb)  +  bmYm(ymb)  =  0 

which  has  non-trivial  solutions  only  if  its  detenninant  is  zero 
J  (v  a)Y  (v  b)-J  (v  b)Y  (v  o)  =  0. 

m  V  m  /  m\  m  J  m  V  m  '  mV  m  ' 


(A3. 9) 


(A3. 10) 


(A3. 11) 

(A3. 12) 


(A3. 13) 


(A3. 14) 


This  is  a  transcendental  equation  with  roots  vmn,  n  =  1,2,...  .  For  these  roots,  the  solutions  of  the 
system  are 
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Jm(vmnb)  n 

Y  (v  a )  "m  Y  (v  b) 

m  v  mn  /  m  \  raw  / 


(A3. 15) 


For  (A3. 1 1)  we  can  then  write 


p„  ip)  =  (v„/>)  -  ^4^7  F" 


(A3. 16) 


and  for  (A3. 7) 


/  e(P><P)  =  Fmn(P)0  e  (?)  ■ 


mn 

o 


m 

o 


(A3. 17) 


These  functions  are  orthogonal,  in  the  sense  that 

e(P’<P)f  e(P’<P)dS  =  AmmSmm'  »  { /  e  (A  <P)f  o  ( A  =  0  .  (A3. 18) 

"  mn  m'n  J  mn  m'n 

D  o  o  Doe 


They  can  also  be  orthonormalized. 
Returning  to  (A3. 5),  we  write 

h"(z)  +  (k2  -  vmn2 )h(z)  =  0 
which  has  solutions 

*1(2)  =  =*“-' 
where 


~VJ  ’ 

k2  >  vm2 

-*\Km2  ~k2  > 

•  Vmn  >  k2 

(A3. 19) 


(A3. 20) 


(A3. 21) 


For  the  upper  root,  the  solution  with  the  plus  sign  represents  a  wave  traveling  along  the  negative 
z-axis  while  the  one  with  the  minus  sign  represents  one  that  travels  along  the  positive  z-axis.  For 
the  lower,  we  have  exponential  decay  along  these  directions.  From  (A3 .4),  (A3. 17),  and  (A3 .20), 
we  can  write, 

I7+  e  (A  cp ,  z)  =  e±iAm,:  Pmn  (p)0  .  (cp) .  (A3 .22) 

mn  m 

o  o 
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These  are  the  eigenfunctions  that  comprise  the  z-component  of  the  electric  Hertz  potential. 
From  this  we  can  construct  the  electromagnetic  field  according  to  (A2.3).  Thus,  the 
eigenfunctions  for  the  z-component  of  the  electric  field  are 


' E*+-e{p,<p,z)  =  vmn2e±i ^ 

mn 

o 


Pmn(P)0  eiV) 


(A3. 23) 


while  those  for  the  p-  and  cp  -components 


TM 


EP±e(p,(p,z)  = 


d2n*  e(p,cp,z ) 

mn 

o 

dpdz 


=  e±'"""'"  Pmn  (P)®  J(P) 


and 


d2^  e{p,(p,z ) 


TM  tt’  (p~ 


E<p  e{p,(p,z)  = 


iX 


pdcpdz 


=  +- 


P 


Pmn  (P)0  e  (?)  ■ 


(A3. 24) 


(A3. 25) 


We  note  that  the  ^-component  automatically  satisfies  the  boundary  conditions.  For  the 
magnetic  field 


pn±  e{p,(p,z) 

ikY  mn  ikY  , 

™ HPt, (AM  =  - - ^ - =  —  e‘“-:  P„ (P)0  ,  M 


P 


dcp 


P 


(A3. 26) 


and 


sn±  eC o,<p,z) 

mn 

™  H’\ (p,<p,z)  =  ikY - - =  ikY Pj (p)0  ,(<P) ■ 

mn  UQ  m 

o'  o 


(A3. 27) 


Again,  we  see  that  (A3. 26)  satisfies  the  boundary  condition  that  the  component  of  the  magnetic 
field  normal  to  the  perfectly  conducting  surfaces  must  be  zero  there.  This  condition  follows  from 
(A  1.2). 

A4.  TRANSVERSE  ELECTRIC  WAVES 

If  in  (A2.1)  we  let  77 (r)  =  0,  we  find  that  the  electric  field  does  not  have  a  component  along 
the  z-axis.  The  ^-component  of  the  electric  field  is  given  by 
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TE  J71’  —  Vfc'7’ 


E<P  (r)  =  ikZ 


dp 


(A4.1) 


For  (A  1.2)  to  be  satisfied,  it  is  sufficient  that  the  normal  derivative  of  the  magnetic  Hertz 
potential  on  the  two  conductors  is  zero;  thus,  we  have  the  boundary-value  problem 


V2M(r)  +  k2M(r)  =  0,  =  =  0. 

da  db 


(A4.2) 


Proceeding  as  in  the  TM  case,  we  write 
M(p,<p,  z)  =  g(p,<p)h(z) 


(A4.3) 


and  obtain 


V,2g(Aff)  ti'(z)  2  2 

g(P,<P)  Kz) 


where  g  satisfies 


P 


d 

dp 


A 

v  dP  J 


+  ^4  +  A2Ag(A^)  =  0- 

d(p~ 


Letting 

g(p,<p)  =  R(p)0(<p) 


(A4.4) 


(A4.5) 


(A4.6) 


where  (Pis  defined  in  (A3. 9),  we  get  that  R  satisfies  Bessel’s  equation  (A3. 10)  with  v replaced 
by  //;  thus, 

R,n  (P)  =  CmJm  (, PmP )  +  d„J,n  (PmP)  •  (A4-7) 

From  the  boundary  condition  (A4.2), 


CmJn!(P,na)  +  d,X(P,na)  =  0 
c  J  (ub)  +  dY  (ub)  =  0 

and  from  this  we  get  the  condition 


(A4.8) 


Jm  (PnP)Y,n  (PX)  -  J ( Pmb)Y, (//„fl)  =  0  . 


(A4.9) 
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This  is  a  transcendental  equation  with  roots  jumn ,  n  =  1,2,...  .  For  these  roots,  the  solutions  of  the 
system  are 


mn  f  L' mn  f  L' mn 

Ym  (Pm„a)  Ym  (jJmnb) 


(A4.10) 


and  for  (A4.7)  we  write 


Y  ( u  a ) 

m  \r*mn  / 


(A4.ll) 


For  (A4.6)  then  we  write 


g  e(p,cp)  =  Rnn(p)0  e{(p). 


mn 

o 


m 

o 


(A4.12) 


These  functions  are  orthogonal  in  the  sense  of  (A3. 18)  and  can  be  orthonormalized. 

As  with  (A3. 19),  the  solutions  of  h  in  (A4.4)  are  given  by 

KnW  =  *±iK-Z  (A4.13) 


where 


Knm  =  < 


\[/c2  —  u  ,  k  >  u 

V  r*mn  ’  r*mn 

■  I  2  7T  2  j  2 

V  ’  r*mn 


(A4.14) 


For  the  eigenfunctions  of  the  magnetic  potential  in  (A4.3),  we  can  write 


M±  e  (A  <p,  z)  =  Q±iKm"z  Rmn  (p)0  e  (cp) . 


(A4.15) 


From  (A3 .23),  the  corresponding  functions  for  the  electromagnetic  fields  are 


d1M± 


'  Hz±e  (p,  ^,z)  =  z-  VxVx  \M±  ez]  = 


Pz2 


+  k2M±e=Jumn2e±i^Rmn(p)0  e(<p) 


(A4.16) 
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d1M±  e 

mn 

TE  Hp±e  (A  A  z)  =  —— -  =  ±iKmn  Q±iK"mZ  rJ  (p)0  e  (<p) 

mn  CDUZ  m 

o  '  o 


(A4.17) 


d2M+~e 

mn  jjs- 

TE  ^  =  ±  —  q±>k""':  Rmn(p)0  e  (<p) 

mn  pOipOT.  p  m 


(A4.18) 


k7  mn  ik7  , 

--r^  =  -— *±iK-z  RM®  M) 

p  dcp  p 


(A4.19) 


mn 

TE  (p,<p,z)  =  ikZ  — -JL  =  ikZ  Q±iK-nZ  Rj  (p)0  e  (<p) . 

mn  OO  m 

o'  o 


(A4.20) 


We  note  that  all  boundary  conditions  are  satisfied. 


A5.  TRANSVERSE  ELECTROMAGNETIC  WAVES 


In  this  case,  the  z-component  of  the  electric  and  magnetic  fields  is  zero.  From  (A3.1)  and 
(A4.16),  this  implies  that  the  fields  have  a  c±:kz  dependence.  From  (A2.3),  the  part  of  the  electric 
field  due  to  the  electric  Hertz  potential  is 


'  E(r)  =  V  x  V  x  [77(r)z]  =  Vx  [V77(r)  x  z]  =  -zV277(r)  +  —  [V77(r)] 

dz 


=  z  k2n(r)  + 


d2IJ(r)  f  dIT(r) 


dz 1  '  dz 


(A5.1) 


Since  the  z-component  is  zero,  we  have  that 


rE(r)  =  V, 


dll  (r) 


(A5.2) 


Write 


n“E±(r)  =  eabe±(p),  p  =  (x,y). 


(A5.3) 


From  this 


0  =  V  •  [  TEM  E*  (r)]  =  V  (e±ife )  •  e2  (p)  +  e±ife  V  •  e1  (p)  =  e±ikz  V  •  e1  (p) 


(A5.4) 
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so  that 


V  •e±(r)  =  0 .  (A5.5) 

Since  in  (A5.2)  the  electric  field  is  expressible  in  terms  of  a  scalar  function,  we  can  let  that 
function  be  e±,fc  F±(p) ,  so  that  we  can  write  for  (A5.2) 

e±(p)  =  -VF±(p).  (A5.6) 

From  (A5.5) 

V2F±(p)  =  0.  (A5.7) 


The  general  solution  of  this  is 

K(P )  =  Um(p)0  e(<P)  (A5.8) 

m  m 

o  o 


where  the  angular  function  is  defined  in  (A3. 9)  while  the  radial  one  satisfies  the  differential 
equation 


p2Um  (p)  +  pUm(p)-rn2Um(p)  =  0. 


(A5.9) 


This  is  an  Euler-type  equation  with  solutions 

U0(p)  =  A0\np  +  B0-,  Um(p)  =  Amp""  +  Bmp  ,  m  =  l,2, 


(A5.10) 


If  we  reconstruct  the  angular  component  of  the  electric  field  using  these  solutions,  we  find 
that,  for  the  boundary  conditions  to  be  satisfied,  we  must  have 


Amam+Bma-m=0 
Ah' "  +  BJfm  -  0 


,  m  =  1,2,...  . 


(A5.ll) 


This  system  has  only  the  trivial  solution  unless  a  -  b  ,  a  case  of  no  interest;  thus, 

A,n=B,n=  0.  m=  1,2,... 


(A5.12) 


and,  for  the  potential,  we  get  the  simple  expression 
V±(p)  =  A0lnp  +  B0. 


(A5.13) 


The  unknown  constants  are  determined  by  the  potential  difference  between  the  two  conductors. 
For  example,  if  the  inner  conductor  is  at  potential  V  and  the  outer  at  zero,  then  the  boundary 
conditions  give 
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An  In  a +  B„  =  V 

0  0  (A5.14) 

A0  In  b  +  B0  =0 

from  which  we  determine  non- trivial  values  for  the  constants. 

If  we  follow  the  same  line  of  reasoning  using  the  magnetic  Hertz  vector,  we  will  find  that  we 
must  have  the  radial  derivative  of  the  potential  equal  to  zero  on  the  two  conductors.  This  will 
result  in  (A5.12)  being  true,  but  also  that  Ag  =  0  ;  thus,  the  potential  will  be  equal  to  a  constant  ( 

B0 )  and  the  fields  will  be  equal  to  zero.  We  see  then  that,  in  the  TEM  case,  the  fields  are 
determined  from  a  single  scalar  function,  the  electric  potential. 

A7.  CONCLUSION 

We  have  derived  the  eigenfunctions  that  may  be  present  in  the  space  between  the  two  perfect 
conductors  of  a  coaxial  line.  The  TEM  mode  will  always  be  present,  irrespective  of  frequency. 
The  presence  of  other  modes  depends  first  and  foremost  on  the  feeding  arrangement.  If  the  fields 
of  the  source  are  independent  of  the  angular  direction,  so  will  the  total  fields.  If  not,  then  a 
number  of  modes  is  possible,  depending  on  the  dimensions  of  the  inner  and  outer  conductor.  We 
note  that  both  in  the  TM  and  the  TE  case  we  have  a  term  (in  =  0)  independent  of  (p. 
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PART  3 

NUMERICAL  SOLUTION  OF  INTEGRAL  EQUATION  FOR  COAXIAL  LINE 

ABSTRACT 

This  is  the  third  part  of  the  report  on  the  formulation  of  the  problem  of  radiation  of  a  coaxial 
line  into  a  half-space  in  terms  of  BIEs.  In  Part  2,  we  showed  that  the  problem  can  be  reduced  to 
solving  a  single,  scalar  integral  equation.  Here,  we  convert  the  scalar  integral  equation  into  an 
infinite  system  of  linear  algebraic  equations.  We  also  express  the  coefficients  of  the  system  in 
tenns  of  single  integrals,  and  proceed  to  show  how  to  compute  them. 
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1.  INTRODUCTION 

We  proceed  with  the  numerical  solution  of  the  integral  equations  obtained  in  Part  2.  The 
geometry  is  as  in  Figure  1.1.  In  Part  2,  we  found  that  the  coefficients  of  the  current  density 
expansion  can  be  detennined  by  solving  a  single,  scalar  integral  equation.  In  Section  2,  we  begin 
with  the  scalar  integral  equation,  substitute  in  it  the  expansions  for  the  unknown  linear  current 
densities  and  perfonn  the  two-dimensional  integrations  to  end  up  with  one  equation  in  an  infinite 
number  of  unknowns.  The  coefficients  of  the  equation  are  Sommerfeld-type  integrals.  In  Section 
3,  we  use  the  orthogonality  properties  of  the  expansion  functions  to  generate  an  infinite  system 
of  linear  algebraic  equations.  The  coefficients  of  this  system  are  integrals  of  four  distinct  forms. 
In  Sections  4  through  7,  we  discuss  how  the  integrals  are  computed  using  Mathematica® 
(reference  1).  In  Appendix  A,  we  evaluate  analytically  the  integrals  used  in  Section  3. 

4  ^ 


<j 


G 


Figure  1.1.  The  semi-infinite  coaxial  line  of  inner  radius  a  and  outer  radius  b  is  fed  at 
z  =  -oo  .At  z  =  0,  it  opens  up  into  a  half-space,  its  outer  conductor  becoming  a  plane 
that  extends  to  infinity.  All  surfaces  are  perfectly  conducting. 
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2.  INTEGRAL  EQUATION 

The  integral  equation  is  given  by  (5.8)  of  Part  2  and  we  repeat  it  here 


1 K  (■ ■)gSff)dS  +  |  h„  (r  f^’dS  =  u„(r’),r’ea 

S„  Su  ^P 


dg(  r,r') 


where 


YV 


K  (O  =  (1  -  A)  -f  —  +  ikYV^B.Fj  (p) 

^yx)  P  n= 1 

a^{x)  na  n=lY0(v0na) 

b ln(/)v  xb  „=i  Y0(v0nb) 


and 


g(r,  r ')  =  -  —  X  cos  K<P~  <P')F,  (p,p',  z,  z') 
4^r  /To 


with 


(2.1) 


(2.2) 

(2.3) 

(2.4) 


(2.5) 


poo 

£0=1;  G=2>  /  =  1,2,...;  Fl(p,p',z,z')  =  J/(ap)Jl(ap')e 

J  0 


z-z'\yja2  -k2 


ad  a 


yja2  -k2 


(2.6) 


We  want  to  accelerate  the  convergence  of  the  integral  in  (2.6).  To  this  end,  we  split  the 
integral  in  two  parts,  a  static  part  and  the  rest.  From  (2.6),  the  integral  corresponding  to  the  static 
Green’s  function  is 


<• 00  |  t\ 

Hj  (p,  p' ,  z,  z')  =\  Jt  ( ap)J,  ( ap ')  e  '  da  (2.7) 

J  0 

and  what  remains  in  (2.6)  is 


G,  (j o,  p',  z,z')  =  F,  (p,  p\ z,z')~  H,  (p,  p',z, z') 

poo 

=  J0  j ,{ccp)J ,{ap') 


-\z-z'\yla2-k2 

a  C  -\z-z\a 

-e  1  1 


\ja2  -  k2 


da . 


For  the  integral  equation  in  (2.1),  we  can  then  write 


(2.8) 
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l  oo 

—  £' Sl  J  ha  Cr)  C0s  l^P  —  V') 

^7T  /=0  s 


„  dG,(a,p',  z,  0) 
dp' 

dH,(a,p',z,0 ) 


-g£,|;,„(r)Cos%^')  ^ 

■  j-  j  */,(■•)  cos/(p  -«>')— (y,  ’ z’ o> 

4^"  /=0  ^  4? 

1  ^  r ,  .  s  „  ndH,(b,p',z,0) 


dS 

dS 


,  Z^/|^(r)cos^-^')' 

4;r  /=o  s„ 


dp' 


dS 

dS  =  wCT  (r') ,  r'  e  a . 


(2.9) 


From  (2.2),  (2.3),  and  (2.4),  the  current  densities  are  independent  of  the  angular  direction.  We 
can  then  replace  (2.9)  by 


£f“  h(rdG°(a’P’’Z’0)dz-“C  h(vdH°(a’P'’Z-0) 

2  J-*  a  pp'  2  J-00  a  dp' 

2J-»  4,'  2J-  7  4?' 


c/z 


<iz  =  w  (r')  ,  r'  e  cr . 


(2.10) 


We  deal  with  each  of  these  integrals  separately.  For  the  second  one,  we  use  (2.3)  and  (2.7)  to 
write 


a  r° 
2' 
ikW 


a  r  o 

/2=-J  VO 

V  J-oo 

B. 


dH0(a,p',z,0) 


YV  fo 


dH0(a,p',z,  0) 


dp' 


dz=  ■’  J  (e^-de^T"0 - 

21n(2')^“°0^  dp' 


dz 


X  n= 1  ^(V,,,,**) 

YV 

Jo 


f 


e«..*  dH0(a,p',z,Q)dz 


dp' 


2  (/ j  SI  *4  -  ^  )  |r  1>-C  («« U.  (<*>’)  e- 


^  »=i  ° 


dp 


r  oo 

yJo  daJ0(aa)J0(ap')e~za  . 


(2.11) 


We  require  that  z  >  £  >  0  so  as  to  make  the  integrals  in  cr  converge  absolutely  and  uniformly. 
This  allows  for  interchanging  operations;  thus, 


h  = 


YV 


lim-^-f  daJ Aad)J Aap')\  &(e'fe-de  lfc)e 

£■->0  f)n'  JO  js  V  / 


2  In  ( 


ikYV  x-— .Z?  ^  r°°  r°°  -3 

- £v,"  T lin)  T7  L  ^«4(««)-/o(«/?')f  dzQ'^':  d 

n  f^Y0(v0na)  dp  V 


l0Ky0,i 

YV  dr ® 

„  7  \  lm?T7Jn  daJ0(aa)J0(ap’) 
21n(^J  *->°  <^9  Jo 


(ik-a)s  -{ik+a)e 

— - A- - 

ik  -a  ik  +  a 
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ikYV 


Z- 


^  r00  ,  .  .  ...  e  1  } 


(7  r™  e 

lin^T^Jo  daJo(aa)Jo (ap')~ 


x  ttY0(v0na)  s^o^p'Jo 
YV  r® 

- - — -  lim  d  aJ0  (aa)Jl  (ap')a 

2\n{%)  '?_>oJo 


//l0„  +a 

(ik-a)s  ~{ik+a)e 

6  ,  C 

- b  ^4 - 

ik  -a  ik  +  a 


ikYV  ^ 

+ — Z 


£)  — (//1q 

B  p<x>  e  v  ' 

lim  daJAaa)Jx{ap')a 

£->0  Jo  7 


(2.12) 


We  re-write  this  as 


h  = 


YV 


2  In  (*) 

ikYV 

2  In  (*) 


p  oo 

lim  d  aJ0  (aa)J \  ( ap ') 

£->0J0 


-e 


+  Ae 


-((&+«)£■ 


p  00 

lim  d  aJ0  ( aa)Jx  {ap') 

£•-»  0  J  0 


(ik-a)s  ~(ik+a)e 

- b  ^4 - 

a  -  ik  a  +  ik 


+ ,J7 1  i^b)  S  I" ' ^  e"(u*'  "* 

^  „=1  Yo(yo„a)‘^°  lK  n+a 


(2.13) 


In  the  second  and  fourth  terms,  the  limits  can  be  moved  inside  the  integrals.  The  first  and  third 
are,  by  definition,  Weber-Schafheitlin  discontinuous  integrals  (reference  2,  pp.  398-410). 
Dropping  the  limit  notation,  and  thus  having  only  one  such  integral,  we  have  that  (reference  2,  p. 
404  and  reference  3,  p.  100) 


poo 

J0  Jfl(ua)Jfl+l(va)da  = 


</,+1),  0  <  u  <  v 


u  v 


1 


U  =  V 


2  u 

0,  u  >  v 


(2.14) 


We  then  have 


YV(A-l)  ikYV  r® 

I7  = - - — - 7 — 7  [  daJ n(aa)J,  ( ap ') 

2  21 n(Z)p'  21n(j)b  P  > 


1 


A 


ikYV  ^ 

+  ^rX 


B 


a-ik  a  +  ik 
da 


np  „=iT0(v0„a)  n  „=1  T0(v0„a) Jo 


M On  +  a 


(2.15) 


The  fourth  integral  in  (2.10)  is  the  same  as  the  second  with  a  replaced  by  b  and  a  reversal  in 
sign  because  of  (2.3)  and  (2.4);  thus,  we  have 
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b  f°  ,  ,  , dH0(b,p',z,0 )  /11T 


o  ru 

ia=-\  hb(r) 

/  J  —  oo 


2  J-°°  ° ' '  dp' 

ikYV  iK„B„ 


dz  = 


21n(^)j0°  daJ0(ab)Jx(ap’) 


1  A 

-  +  - 


a  -  ik  a  +  ik 


IKI V  vp  lAqnDn  r»  r  .  ,  . 

Zj  v  /  7  \  in  d0{C(b)Jl(C(p  ) 

*  tiY0(v0nb)J° 


da 


i\n  +  a 


We  turn  to  the  first  integral  in  (2.10) 
a  r°  ,  ,  .dG0(a,p',z,0) 


n  /»u  , 

/i=^|  Mr)- 

/  J  —00 


4?' 


<iz 


=  f  II  t/z/z« (r)  ^7  j”  Jo  («a  Vo  (M) 


ae 


z'ja  -k 


\la2  —  k2 


-e' 


da 


IT  ro 


2 1„  ( ^  £ ' *  (e_lt  • -  ■ A' )jp'  i  MaaWap') 

■  — L,*'  J“  <fc /,('  J,ina)J„Utp') 

X  ttY0(v0na)J-^  dp  Jo 


ae 


4a2-k2 


ae 


\la2  —  k2 


-  e~ 


-e' 


da 


da  . 


We  re-write  this  as 


£  = 
ikYV 


2i^y  ^  i°  daj°(aa)j°  <^>j:  *  ( e":  -  ^  e_“ ) 


ae 


zda2-k~ 


B„  d 


'■4»z 


\Ja2  -  k2 


-e 


ae 


\ja2  -k2 


YV  dr 

2h^)^Jo  daJ0(aa)J0(ap’) 


-a 


+  - 


-e 


A 


\Ja2  -  k2  v  ik  -  yja2  -k2  ik  +  \Ja2  -  k2 


1 


ik  -  a 


—  Z  f  ,£  daJ0(aa)J0(ap') 

n  t(Y0(v0na)  dp  Jo 


a 


\la2  -k2  {i/ l0n  +  sja2  -k2  j  M« 


+  a 


We  compute 


a  1 

Va2  -  k 2 1  ik  -  \Ja 1  -  k2  j  ^  _  a 


\la2  -k2  [ik  -  yja2  -k2  j  -  a(ik-a ) 
\/a2  -  k 2  1 -  Va2  -k2  j (ik-a) 


(2.16) 


(2.16) 


A 

ik  +  a 

(2.18) 
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Similarly,  for  the  third  integral  in  (2.10),  we  have 
hr0,.  ^dG0(b,p',z, 0) 


n  < 

=  A  k(t)- 

1  J  —  00 


2  J-°°  dp' 

YVk2  d  r°° 


-dz 


I  V  K  O  poo 

- — —  daJJab)JJap') 

21n(j)4?'J°  oV  '  oV  ^ 


« 


+  ik  +  \J  a~  —  k 


V«2  -  k 2  (ik-a)[^a2  -k2  +aj 


a 


+  .4 


(2 


—  ik  +  yj  a~  —  k~ 


yja2  -k2  {ik  +  a)[yja2  -  k 2  +aj 

poo 

da 
'  Jo 


ik3YV^ 

^  tiY0(v0nb)  dp 


J0  (ab)J(;  (ap')  |//t0n  +a  +  \la2  -k2  j 
yja2  -k2  [i/ lon  +\Ja2  -k2  j(i/ l0)i  +a)|«  +  V«2  -k2 


In  the  next  section,  we  convert  (2.10)  to  an  infinite  system  of  linear  algebraic  equations. 


(2.23) 
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3.  SYSTEM  OF  EQUATIONS 

We  write  (2. 10)  in  the  fonn 

-(A  +  h  +  h  + 1  a  )  =  K  (O  ,r' e  a 


(3-1) 


where  the  four  terms  on  the  left  are  given  by  (2.15),  (2.16),  (2.22),  and  (2.23),  while  the  term  on 
the  right  is  given  by  (2.2).  We  proceed  to  collect  coefficients  of  the  unknown  terms 


A 


YVk  d  r®  "o 

Jo 


^  poo 

da 
/  Jo 


a 


-  ik  +  \la2  -  k2 


+ 


21n (%)  dp  a^Ja2  - k 2  (ik  +  a)[^a2  - k 1  +aj 

YVk 2  d  r«  J0(ab)J0(ap')la  -  ik  +  Ja2  -k2  J 

- 7 - V -  da - ,  - .  - f-t - 7 - r — 

21n {%)  dp  0  a^ja2  - k 2  (ik  +  a)[^Ja2  - k 2  +aj  21n (z)p 

ikYV  r®  J0(aa)Jl(ap')  ikYV  r®  J0(ab)Jl(ap')) 


21n(^) 

ik2YV 


YV 


n 


YB. 


a  +  ik 

1  3 


21n(/f ) 


p  oo 

da 

Jo 


a  +  ik 


poo 

da 
'  Jo 


J0(aa)J0(ap')\ 

[iK 

7  +  a  +  yfa 

■2-lc) 

1 

yja 2  -  k 2  j 

iA0n  +  \la~  — 

\ 

¥) 

(i\n  +«) 

[a  +  ^ja2  -  k 2  j 

3 


poo 

da 
'  Jo 


J0(ab)J0(ap')  | 

\is l0„  +a  +  \Ja 

■ 2-k 2) 

1 

^Ja2  -  k2  [i/ l0n  +yja2  - 

kl) 

(//L0„  +  «) 

| a  +  yja 2  -k2  j 

i/ 1, 


'On 


®  J {]  (aa)J]  (apr)  i\ 


2 


poo 

da 

Jo 


'On 


i^on+a  k  Y0(v0nb) 


poo 

J0  J0(ab)Jl(ap1)- 


da  n  „  <  ,  . 

- ,  2  P0«  (P) 

U0n  +a  k  i 


rur  3 
2 In (z)  dp 


poo  ( 

da  — 
'  Jo 

a 


J0  ( aa)J0  (ap')^a  +  ik  +  \Ja2  -  k2 


-k2  ( ik-a)i\ja 2  -  k 2  +a 


2  ,  2 

a  -k 


YVk 2  3  JMb)JQ 

21n (%)  3p  0  a^Ja2  - k 2  (ik-a)lja2  - k 2  +a 


YV 


ikYV  r®  J0(aa)Jl(ap')  ikYV  r®  J0(ab)Jx(ap') 


poo 

da 

I  Jo 


2  In  (x)p'  21n(j)Jo  a  -  ik 

We  combine  terms  above  to  get 


21n(j) 


poo 

da 

Jo 


a  -  ik 


(3.2) 
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A 

2]n(z) 

poo 

+  ik\  da 

Jo 


k2^r  da 
dp'  ^ 


[ J0  (aa)  -  J0  (ab)]  J0  ( ap ')  -  ik  +  \]a2  -  k2  j  ^ 


ayja 2  -k2  [ik  +  a)[ \ja2  -  k2  +  a 

[/0  (aa)  -  J0  ( ab )]  Jx  (ap') ) 
a  +  ik 


P 


ik3  ® 

+—Yb, 

n  „=i 


_d_ 

dp 


poo 

da 
'  Jo 


J0  (aa) 

J0(ab) 

J0(ap')[iX0n  +  a  +  \Ja2  - k2  j 

J0(V0na) 

1  7o(n>^)J 

\ la 2  -  k 2  j 

iX0n  +\la2  -k2  J 

(a0„  +a) 

+  \ja2  -  k2  j 

+ 


iX, 


'On 


k  Yo(VOna)P'  k~ 


p  00 

da 

Jo 


J0  (a  a)  J0  (ab) 
Yo(vona)  Y0  (v0nb) 


L0\v0nL 

iX0n  +  a 


Jx(ap') 


-y^iP) 


21n(j) 


n  x  [j0(aa)- /0(ab)]/0(a//)  a  +  /£  +  V«2  ~k2 

- Jc 2  —  [  da - - \  - r 

^P  °  a\ja2  -k2  {ik  -a)Na2  -k2  +a\ 


1  f 00 

H - /«:  da 

p'  •’° 


[  J0  (aa)  -  J0  (ab)]  Jx  (ap') 


a  -  ik 


(3.3) 


To  convert  the  integral  equation  to  a  system  of  linear  algebraic  equations,  we  first  integrate 
both  sides  of  this  with  respect  to  p’  from  a  to  b.  In  effect,  we  perform  the  integration 


rb  1  rab  \  rab  f 

Jx  (ap)d p  =  Jx  (t)dt  = - J0  (t)dt  = 

Ja  J  aa  py  J  aa 


1  cab  T  ' 

a Jaa 


J 0(aa)  -  J 0(ab) 


a 


(3.4) 


and  also  use  the  results  in  Appendix  A.  The  end  result  is 

?  f 00  I 

-k  da- 

Jo 


A 


21n(/) 


[J0(aa)  -  JQ(ab)] 

r(- 

7  - 

ik  +  yja 2  -k 

7) 

a^ja2  -k2  ( a  +  ik 

)(v 

cr  -  k 2  +aj 

1 

In  (/)+/&]"  da 


[J0(aa)~  J0(ab)]~ 


a 


(«  +  ik) 


ik3  00 

+—Yb, 

n  B=i 


poo 

-  da- 

Jo 


\jQ(aa)  -  J0(ab)] 

J0  (aa)  J0  (ab)  [ 

iX0n  +  a  +  \J a~  - 

/oK«)  ^o(VOJ' 

\ja2  -  k 2  j 

iX0n+\ja  —  j 

(a0„  +  a 

)(a  +  ^a2  -  k 2  j 

1 
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AM.  +iArrfg. 

k  Y0(v0na)  k~  Jo 


[  J0(aa)  -  J0(ab )] 

J0  (aa) 

J0  (ab)  1 

_Yo(VOna) 

^o(VOn^)  J 

a(a0„  +  a) 

21n(j) 


2  r00 

-k  da 

Jo 


|/o(aa)  "A) 


2  i  2 

a  -k 


a^J a2  -  k 2  (t 


a-ik) I  V«2  -  k 2  +a 


Jo  a  (a-ik) 


(3.5) 


We  next  multiply  (3.2)  by  p'P0m  ( p')  and  integrate  from  a  to  b 


A 


x  in  Car) 


2  f00 

k  da 

Jo 


[  J0(aa)  -  J0(ab)\ 

J0  (a  a)  J0  (ab) 

|a  -ik  +  \Ja2  -k2  j 

J0(v0ma)  Y0(v0mb)_ 

yja 2  -  k2  ( 

r 

a  +  ik^(y]a2  - k2  + 

_  -i  ^ 

l(«2  -Mm2) 

a 


<•00 

-ik\  da- 

Jo 


[  J0(aa)  -  J(l(ab)] 


J0  (aa)  J0  ( ab ) 
Yo(y0ma)  Y0(v0mb) 


(a  +  ik)(a2  -vj) 


a 


ilk 3 


n 


.2  S*. 

K=1 


poo 

da 

Jo 


J0  (aa)  J0  (ab) 

J0  (aa)  J0  (ab) 

(z'2o»  +ar  + 

Va2  -  k2  j 

a2 

J0(v0„a)  Y0(v0nb)_ 

1 - 

3 

a)  r0(v0mfc)J 

\la2  —  k2  ^iA()n  +  Va2  -k2  j 

(iX0n  +a) 

|a  +  Va2  -£2  j 

(«2  -n0m2) 

iX, 


p  oo 

da 

Jo 


J0  (aa) 

J o(ab) 

J0  (aa) 

J0(ab) 

Yo(V0nb)_ 

Y0(v 0ma) 

Y0(yomb)_ 

a 


(iAv„+a)(a2  -vj  t 


k 2 


1  _  1 

IftWO]2  [^o(vb,«)]2 


s„ 


n 


In  (z) 


?  r00 

&  da 

Jo 


[J0(aa)  -  J0(ab)] 

J0  (aa)  J0  (ab) 

|a  +  ik  +  -\Ja2  -k2  j 

J0(v0ma)  Y0(v0mb)_ 

yja2  -  k 2  ( 

a -ik)[^ja2  -k2  +aj 

l(«2  -n0„,2) 

a 
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+ 


da 


[J0(au)  -  J0(ah)] 

J0  (a  a) 

J0(ab) 

iYo(VOma) 

Y0(y0J)_ 

(a  -  ik) 

(  2  2 

-v0m 

m  =  1,2,... . 


(3.6) 


The  last  two  expressions  constitute  an  infinite  system  of  linear  algebraic  equations  in  an 
infinite  number  of  unknowns  (A  and  the  Bn ).  If  in  (3.5)  we  set  all  the  Bn  equal  to  zero,  then  we 
get  for  A 


h  r 

a(a-ik)  Jo 


[J0(aa)-  J0{ab)\  [a  +  ik  +  ^Ja2  -k2  j 


-k  +  a 


A  = 


ln(/f)  +  z^J(  da 


[  J0(aa)  -  J0(ab)]~ 


a 


(a  +ik) 


k2  [da 

Jo 


a^Ja2  -k2  (a -ik){yja2  -k 
[  J0(aa)  -  J0(ab)]~  (a  -  ik  +  V«2  -k 


a\ja2  -  k 2  (<2  +  ik)  V 


a2  -  k 2  +a 


(3.7) 


We  can  combine  the  integrals  to  reduce  this  to 

Hx)-ikr  dai^-j^f 

^  _ _ || _ a\jan  -k1 

In (x)ukrdaW^Z^L  ' 
y  ’  J» 


(3.8) 


The  issue  we  must  address  now  is  whether  we  should  use  such  combinations  in  (3.5)  and  (3.6). 
What  we  have  accomplished  in  splitting  up  the  Green’s  function  is  to  increase  the  rate  of 
convergence  of  the  integrals  in  a  by  an  order  of  magnitude.  Had  we  not  split  the  Green’s 
function,  the  integrals  in  (3.8)  would  behave  as  a  2 ,  as  a  —>  .  In  (3.7)  and  (3.8),  the  integrands 

behave  as  a  3 .  We  have  thus  gained  an  order  of  magnitude  in  the  rate  of  convergence.  In  fact,  in 
the  integral  that  does  not  involve  the  static  term,  the  behavior  is  a4.  In  the  static  term,  we 
altered  the  original  behavior  of  a  2  by  integrating  part  of  the  expression,  with  the  rest  behaving 
as  a  3 .  This  is  true  of  all  four  terms.  Since  the  acceleration  of  convergence  has  been 
accomplished,  we  should  feel  free  to  recombine  the  remaining  terms.  The  only  reservation  here 
is  the  following.  Looking  at  the  numerator  of  (3.7),  we  see  that  the  first  integrals  behaves  as  a  3 , 
while  the  second  can  be  split  into  two  integrals,  one  behaving  as  a~A  and  the  other  as  a  3 . 
Whether  the  last  two  terms  are  much  smaller  than  the  first  cannot  be  concluded  from  these 
observations.  We  proceed  then  to  combine  terms  in  (3.5)  and  (3.6). 
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First,  we  perform  the  following  calculations 


^0  n 

iA0n  +  a  +  a 

Ja2  -  &2 

k2a(a0n+a )  v 

'a2  -k2  [i\n  +V 

(a0„  +  «) 

|«  +  V 

*a2  —k2  j 

Min  ^ 

a2  -  k2 1 

*'4>B  +>/ 

a2  -  k2  j 

+  V 

& 

|-* 

2a| 

+  a  +  m 

1 a2-k2J 

k2a\ 

/«2  -k2 

a2  -k 

J) 

(*'4>„  +«) 

+  V 

b2-£2) 

iA0nJa2  -k2 

+  Ja2  -k2  J 

|a  +  iA0n\la2  - 

k 2  +a2  - 

1  1 
(N 

~k2a(iA o„ 

+  a  +  J  a  —  k~  j 

k2aja^  -  k2  |//tDn  +\Ja 2  -  k 1  j 

(U0n  +«) 

+ Ja2  -  A'"  j 

1 

iA0nJa2  -k2 

+  Z>^QW  +  yj  (X  —  k  j 

a  +iA0n\[t 

z2  -£2  -k2 

-  k2a  [i/ l0)I  +  a  +  Ja2  -  k 2  j 

k2aja2  —  k2  \iA0n  + 

a Ja 2  —  k2  j 

|(a0;,  +«) 

+  Ja2  -  k 2  j 

1 

+  z'/L0;(  +  Ja2  ~k2  j 

\a 

iA0iiJa2  -  k 2  - k 2 

<N 

1 

<N 

Ljs 

SC 

+ 

A 

sc 

i _ H _ i 

1  1 

<N 

1 

CN 

1 

(N 

k2aja2  - 

~k2\ 

i\„  +  Ja2-k2  j 

(a0„  +  «) 

+  V 

a2  -k2  j 

1 

iA0nJa2  -  k 2  -  k 2 

yz  +  a0n  +  \ 

la2-k2]) 

a  +  iA0llJa 2  -  k 2 

k2aja 2  —  k2  [iA 10/ 

+  Ja2  —  k2  j 

(iA0n  +a) 

[a  +  Ja2  —  k2  j 

i/L0n  ■Ja2  -k2  -  k2 

[a  +  Ja2  -l 

c2  j  a  +  iA0n  |  a  +  Ja 2  -k2  j 

k2aja 2  —k2  [iA0n  +  Ja2  -  k 2  j 
iA0lIJa2  —  k2  —  k2 

(a0„  +  «) 

a  +  Ja 2  —  k2  j 

k2aja 2  —k2  [i\n  +  Ja 2  -  k2  j 


Also, 


k2a\a 

-ik  +  Ja2  -  k 2 

)  ika 

Ja2  -k2  (a  +  ik){ja2  -k2  +aj 

(a2 -v0m2) 

|  (a  +  ik) 

(«2  -yJ) 

k2a\ 

'  55 

1 

+ 

'  2  i  2 

a  -k 

j  -ikaja2  -k2 

^Ja2  —  k2 

+  aj 

Ja2  —  k2 

(a  +  ik)(J a2  - k2  +  aj| 

(  2  2  ^ 

a  -Vo  m  J 
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i/ L0n  +  a  +  \la2  -k2  \a2 


a0„a 


la~-k~  iA0ll+yja  -k~ 


+  a)[a  +  4cc2  -k2  j(a2  -  v0m2  j  k2  (iA l0n  +a)(a2  - v{)2  j 


i'/L„  +a  +  yja 2  —k2  )ak2  -iAn„\laz  -  kz  I z'/L  +  V«2  ~kz  II a  +  \laz  -k 


k24ct2  -  k 2  |  //L  +V«“  -kz 


+  a)  a  +  \Ja2  - k 1  Ma2  -v0m2 ) 


/!„  +  a  +  -  A'2  a&2  -iA~yJaz  -  kz  I  //L  | «  +  \jaz  -kz  |  +  -  A:2  +az  -k 


2  /  2  .  2  /  2 


&2Vtf2  -k2  y\n  +yja2  -k2  j(//t0n  +  a)|ar  +  ^Ja2  - k 2  j(a2  -v0w2) 

|//10b  +«  +  \Ja 2  —k2  j ak2  - iA0n \Ja2  -k2  a[a  +  iA0n  +  4  a 2  - k2  j  +  iA0ll4~a2~ 
k2  4  or  -  k2  |//L0n  +4 a2  -k2  j(/20„  +  a)|a  +  V«2  -&2  j(«2  -v0m2) 


2  ,2  i  2 

a  -k  -k 


a  +  z/L  +  VeT  -  £2  (A:2  -iAn„\laz  -  k'  \a  -iAn„\Ja  -kz  \  iAnAaz  -kz  -  k 


2  ,2  ,2 


+  «)(«  +  4a2  -k2  Vcr  -v0m2) 


k^a~  -k  z/L  +y]cr  -k~ 


a\iAn.,\la  -k  -k~\\  a  +  iAn„  +  Vcr  -k  \a  +  iAnJa~  -k 


k24oc2  ~k2  \  iA0n  +4 a2  ~k2  j [iA0n  +  a)ya  +  4cc2  -  k 2  j(a2  -v0m2j 
{iAGn4oc2~^-k2  -  k2  +  4a2  -k2  j(«  +  /20b  ) 
k2 4 oc2  —  k2  (iA0n  +4 a2  -k2  \ (iA0fl  +a)  a  +  4a 2  - k 2  \ (a2  -  v0(„2  j 


a  iA0ny]a2  -  k 2  -  k 2 


k2rrr?r.„+rrr?\(a2  ~vj-) 


(3.11) 
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k2  a  +  ik  +  ya2  - k2  )a 


a 2  -k~  (a  -  ik)\ya2  -  k 2  +  a)(a2  -v, 


)  («  "  ik)[al  ~  voJ )  V«2  ~k2  ( a 2  -  vj ) 


(3.12) 


We  substitute  these  results  in  (3.5)  and  (3.6).  From  (3.5) 


- 7 — A\n(x)  +  ik\  da 

2Hx)  Jo 


[  J0  (a a)  -  J0  (ab)\~ 
aJa 2  - k 2 


+*fB  L_hW_ 

"I  Y0(v0na) 


/•oo 

da- 

Jo 


[  J0(aa)  -  J0(ah)] 


J0  {a  a)  J0  ( ab )  / . 


Y0  (y0„a)  Y0(v0nb) 


iAn„da2  -  k2  -k2 


ay  a2  -k~  \i/ L„  +ya1  -  k 2 


= \ — -<  In {y)-ik[  da 

2  In  (z)  W  J» 


[J0(aa)~  J0(ab)f 
aJa2  -  k 2 


(3.13) 


while  from  (3.6) 


-A  \  £  da 


[  J0(aa)  -  J0(ah)] 


J0(aa )  J0(ocb) 
Y,(y0ma)  Y0(v0mb) 


la 2  ~  k 2  (a2  -v0m2) 


21nMy  B  — _ 1  s 

n  -  I  [Y0(v0„a)]~  j 


/•CO 

da 

Jo 


J0(aa)  J0(ab)  J0(aa )  J0(ab )  /  i2  _^2  2\ 

Y0iy0na)  F0(Ko„h)J[70(i/0ma)  Fo(n0,„h)Jv  °"  _ r_ 

a2  -k2  ia0n  +Ja2  - k 1  Ma2  -v0m2) 


/•CO 

=  da 

Jo 


[/0(aa)-  J0(ab)\ 


J0(aa )  J0(cxb) 
Y0{v0ma)  Y0(v0mb ) 


/a2  -k2  (a2  -v0„,2) 


(2  ,  /«  =  1,2,... . 


(3.14) 


93 


NAWCADPAX/TR-20 13/115 

The  last  two  expressions  constitute  the  final  system  of  equations. 
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4.  CALCULATION  OF  THE  FIRST  INTEGRAL  IN  (3.13) 
The  first  integral  in  (3. 13)  is 

j_  e<°da  [Maa)-Jo(.ab)f 


(4.1) 


The  square  root  in  this  (and  all  other  expressions)  must  be  chosen  carefully.  It  appeared  first  in 
the  Green’s  function  Fourier-series  expansion,  specifically  in  the  integral  in  (2.6).  For  this 
integral  to  converge,  the  square  root  must  be  real.  With  s  >  0  but  small,  we  replace  the  original 

expression  by  yja2  -  k2  +  is  .  This  root  must  have  a  positive  real  part  for  a  >  k 


\ la 2  —k2  +is  =  (a2 -k2^ 
«  \la2  -k2 


+  A 


.•  f 


1  + 


i  i 

~  i—  tan  | 

4  g  2  {a  -k 


£ 

a2  -  k‘  j 


a  >  k 


(4.2) 


which  indeed  has  a  positive  real  part.  For  k  >  a  ,  we  must  choose  the  branch  of  the  arctangent 
function  that  gives  a  positive  real  part  and  has  a  dominant  imaginary  part  as  s  — »  0+ .  The 
appropriate  branch  is  the  one  between  n  /  2  and  'in  12 


^ la 1  -k2 


+  l£  = 


\a2-k2) 
iyjk2  -  a2 


+  £ 


1- 


_i  ,i 
4  l2 

e 


£ 

k2  -a 2 


k2-a- 


2  4k^ 


+  i\lk2  -  a2  ,  k 


>  a . 


(4.3) 


a 


A  simple  calculation  shows  that  this  is  equivalent  to  setting 
yja2  -  k2  =  i\/k2  -  a2  -i£  . 


In  general,  we  write 


\la2  -  k2 


yja2  -  k?  , 
i\Jk2  -  a2  , 


a  >  k 
a  <  k 


(4.4) 


(4.5) 


Using  this  information,  we  write  for  (4.1) 


T  -  _/  [k  sir,  [Jo(aar)-  ^oC^)]2  +  p  r]r/  K  (a  a)  -  J0  (ah)]2 

J°  ~sjk2  —  ~2  k 


ic 


-  -il ,  + 12 . 


(4.6) 


a^Ja 
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The  first  integral  in  (4.6)  can  be  computed  using  a  polynomial  approximation  of  the  Bessel 
function  (reference  4,  pp.  369-370).  Since 


.  a  +  b 

k - 

2 


<1 


(4.7) 


we  see  that  the  argument  of  the  Bessel  function  is  smaller  than  two.  We  elaborate  further  on  this. 
Let 


a  -  pb  ,  0  <  p  <  1  (4.8) 

Then,  in  place  of  (4.7),  we  have 

kb  <  — -  <  2  .  (4.9) 

1 +  P 

We  also  have  that 

k<vm  (4.10) 

where  v0l  satisfies 


■/o(1/o„«)roM)-^o(1/o,1^o(,,o„a)  =  0.  n  =  1,2,...  (4.11) 

Using  the  values  of  v0l  from  reference  4,  Table  9.7,  p.  415,  we  can  provide  bounds  for  kb  and  k. 
The  results  are  shown  in  Table  4. 1 . 


Table  4. 1 .  Bounds  on  kb  and  k. 


p  =  a  b 

Ui  (m’  ) 

kb< 

k  <  (m'1) 

l 

II 

& 

0.8 

12.56 

1.111 

12.56 

1.25 

0.6 

4.697 

1.250 

4.697 

1.67 

0.4 

2.073 

1.429 

2.073 

2.5 

0.2 

0.7632 

1.667 

0.7632 

5 

0.1 

0.3314 

1.818 

0.3314 

10 

An  alternative  way  is  to  use  Mathematica®  (reference  1)  to  compute  this  and  all  other 
integrals.  This  is  what  we  did  in  practice.  We  first  re-write  the  first  integral  in  (4.6)  as  the  sum  of 
two  integrals 

1  \jQ{aa)  -J0(ab)]2  \Jk 2  -a2  1  p \J0(aa)  -J0(ab)Y  a  1 

/,  =  ,  - - - da  +  —  - ,  — da  =  — 

k2  Jo  a  k 2  Jo  je^a2  k‘ 

(4.12) 


(In  +4) 
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For  the  second  integral,  we  use  the  transformation 

t2  =  k2  -  a2 


(4.13) 


to  re-write  the  integral  in  the  form 
[  J0  {aa)  -  J0  (ah)]2  a 


7i2  I 


-da  = 


a" 


£  |yo  (Vh2  - 12  a)  -  J0  (Jk2  -t2b) 


dt . 


(4.14) 


For  the  second  integral  in  (4.6)  we  use  the  following  transformation  (reference  5,  p.  174) 
t2  =  a2  -k2 
to  get 


(4.15) 


w; 


J0(b\jt2  +  k2  )-J0(a\lt2  +  k' ) 


t2  +  k2 


dt . 


(4.16) 


We  split  the  interval  of  integration  into  two  intervals 


W< 


J0(b-jt2  +  h2  )-J0(a\Jt2  +  k2 )  J() (b\]t2  +  k2 )  -  J() (a^t1  +  k2 ) 

- - . — tv - —  c/t  +  J  - - ; — tv - —dt 


t2  +  k2 


t2+k2 


I 21  +  ^22 


(4.17) 


where  L  is  a  large  positive  integer.  In  Mathematica®  (reference  1),  we  set  it  equal  to  10,000.  The 
last  integral  in  (4.16)  is  computed  by  using  the  asymptotic  form  of  the  Bessel  functions 


(4.18) 


J0(z)  ~  J — cos(z-;r/4),  z  — »  oo  . 
nz 


Thus, 


[Mba)-J0(aa)f  =  +  (26a) 

na  [a  b  a  b 

^cos[a(h-a)]  +  sin[a(h  +  a)] 

yfab  U; 


(4.19) 


and 
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1 


1  + 


a 


\Ja2  —  k2 


or  j  1- 


(k\ 
V a) 


2 

(2 


\JXJ 

2 


We  substitute  in  I 


22 


w; 


J0(b\jt 2  +k2)-J0(a\Jt2  +k 2) 


t2  +  k2 

da  1 


dt  =  [ 


[J0(aa)-J0  («/>)] 


(2 


\Ja2  -k2 


da 


11  lb. 

7i\a  a3  27T^L2  +k2^j 


r  i  r 

— I — 

b  y 


We  substitute  (4.12),  (4.14),  (4.17)  and  (4.21)  in  (4.6)  to  get 

T _  rk  [J0(aa)-J0(ab)f  0(a  a)- J0{abj\ 

Jo  (Ti  i  ik  I  2  7T 

adk  -a  ad  a  -k 


k 


\^[Uaa)  JMb)]4k-  a~da+W 


dt 


l 


J0  (. byjt 2  +k2  )  -  J0  (< ayft 2  +k2) 


r  +k2 


-  dt  + 


2^(z2  +k2) 


fi  n 

— +- 

\a  b 


This  concludes  the  description  of  how  we  evaluate  the  first  integral  in  (3.13). 


(4.20) 


(4.21) 


(4.22) 
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/ (a; a,b;  v0n )  =  [J0(aa)  -  J0 ( ab )] 


J0(aa )  JQ{ab) 

Jo(v0na)  Y0(v0nb)_ 


(5.3) 


and  (4.5),  we  write 
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Won2'!' 

k 2 -a2  - 

■«2v 

/  2 

- 

-k2 

iay 

Jk 2  -  a2  | 

(^o„2 

-a2  j 

o'  -k2 


a 


(v'o „2-«2)  -Jk2  -a2  (v0„2  -a2) 


Then 


f k 

=  da 

Jo 


f(a;a,b;v0„)  ( 

A 

K2-k2-Jk2 

-a2  -k 

1 

/'av 

,/:2  -  a2 

'v0n2  -k2  +i^k2  -  a2 

) 

=  v, 


2  rk  f(a;a,b;v0n)da  .  /  2  77 Ck  f(a;a,b;v0n)ada 


On 


a 


/  -«2) 


+  ?'Vvo„2  -  k2  Jo 


-Ik2  -a2  (v0„2  -  a2' 


=  l/0„2/ll  +WV0„2  k2  I \2  ' 


Let 


=  \jk2  -a2  ,  a  =  \Jk2  -s2  ,  sds  =  -ada . 


Then 


T  _  [k  f(a;a,b;v0n)ada  _  ^ /(Vfc2  - s2  ;a,b;v0n)ds 

12  Jo  /Ti  IT/  2  ?\  Jo 


4k^aI{ 

and  (5.6)  becomes 


L)„ 


> 


Vo„2-^2+^2 


r  _  2f k  f{a-,a,b-,v0n)da  .  /  2  17 f * /(V k2  -s2;a, b; v0n )ds 

1_Vo"Jo  2  _ 2  \  +AV0„  *  J0- 


/  2  2 

a^0„  -« 


Von-k2+S2 


We  also  write  the  second  integral  in  (5.4)  as  the  sum  of  two  integrals 


poo 

I2  =  J  da- 


f(a; a,b;v0n)\  Jv0n2 -k2  4a1  -  k 2  -  k 2 


\]a2  -  k2  ^Vq2  -  k 2  +  yja2  -k2  j 


a 


=  da 

Jk 


/(a;- 

a>b;v0n)\ 

^o„2  -^2V 

^a2  -  k2  -  k2  j 

aV 

1  2  , 2 1 

a  -k  | 

^o,,2  ~k2  +7 

(5.5) 


(5.6) 


(5.7) 


(5.8) 


(5.9) 
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(5.10) 


where  A  is  a  large  number.  The  last  integral  is  evaluated  asymptotically,  while,  for  the  first,  we 
use  the  transformation 


tdt  =  ad  a 


(5.11) 


to  get 


We  choose  K  so  that  the  upper  limit  of  integration  is  a  large  integer  L;  thus, 

L  =  4k2 -k2  ,  or  K  =  yjL2  +k2  .  (5.13) 


In  the  second  integral  in  (5.10),  we  may  rationalize  the  denominator  to  get 


(5.14) 


where  the  lower  limit  of  integration  is  large  enough  so  that  the  integrand  does  not  exhibit  any 
singularities. 

We  replace  (5.9),  (5.12),  and  (5.14)  in  (5.4)  to  get 
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+ 


«,  f(a;a,b;v 0l,)| 

(a2V 

^o„2 

-  k2  - 

2  'Ja1  -  k 2  j 

A2  A  /  2 

(2V«  - 

-£2I 

(«2  -^o„2) 

doc  —  J\  +  iJ2  + 1/3  +  J 4 . 


(5.15) 


The  first  three  terms  are  computed  numerically  in  Mathematica®  (reference  1),  while  the  last 
one  is  evaluated  asymptotically. 

The  asymptotic  evaluation  proceeds  as  follows.  From  (5.3) 

....ill  1  sin(2  aa)  sin(2  ba) 

f(a;a,b;vn„)  = — \ - + - + - 1 - -  + - - - - 

net  \  aY()(v0na)  bY0(v0nb )  aY0(v0na)  bY0(v0nb) 


1  1 

cos  a(b-a)  +sin  a(h  +  a)  . 

L  v  L  v  1  >  ) 

Yo(y0nb)_ 

j  1  T/  1  M/  1 

V  7 

(5.16) 


while 


a2\lvon  -k2  ~  k 2 

a^a2  - k 2  (a2  -v02  j 


aJv~  -k  -vn 


1- 


a  Jl- 


"-T 


1- 


fy  V 

V0  n 


On 


On 


1- 


°-2^aJv0n2  -k2  -v, 

:a"3|a>/von2  -k2  ~Vi 
«Vvo„2  -k2  -  v0,2  o„2  ~k 


2{cc) 


1  + 


Uk^ 
2  [a) 


v  a  j 


1  + 


(v  V 

V0  n 

v  a  J 


f  /A2 


2 \a) 


1  (  k  \ 
1  +  -  - 
2 \a) 


2  A,  \2 


v, 


On 


v  a  j 


a 


or 


2  2 

Von 

a 3 


(5.17) 


Thus, 

f{oc\  a,b ;  v0„ )  (a2  Jv 02  -k2  -  v02  yja2  -  k 2  j 

a^ja2  -k2  (a2  -  v0„2  j 

1  J  1  1  sin(2aa)  sin(2ha) 

na  \aYo(vona)  bYo(vo»b)  aYo(vo„a)  bYo(vo„b) 
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1  1  cos  a(b-a )  +sin  or(fe  +  a)]j  ^ v0ii2  - k2  v0n2 

Y0(v0na)  Y0(v0nb)  \  sfab  J  a2  «3 

h;-k2\  1  ,  1 

nal  [ aY0(v0na )  bY0(v0nb) _  ' 


From  this  and  (5.15) 


(5.18) 


(5.19) 


This  concludes  the  description  of  how  we  evaluate  the  second  integral  in  (3.13).  From  (5.15)  and 
(5.19),  the  final  expression  is 


f(a;a,b;v0n)da 

-«2) 


+  i 


f  (a;a,b;v0n)ada 


(5.20) 
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6.  THE  FIRST  INTEGRAL  IN  (3.14) 

We  proceed  to  evaluate  the  integrals  in  (3.14).  The  first  integral  in  (3.14)  is 


=  r 

Jo 


[j0(aa)-  J0(ab)] 


J0(aa )  J0(ab) 


\Ja2  -k2 


Yo(v0ma)  Yo(vomb) 

) 


-ad  a 


a~  -v. 


(6.1) 


We  write  this  as 


=  r 

Jo 


[J0(«a)-  JQ(ab)\ 


J0(aa )  J0(ab) 
Yo(voma)  Yo  (Y0mb) 


4a2  -k 2  (a2  -v0;„2) 


ad  a  =  f 

Jo 


f(a;a,b;v0m)ada 
4  a2  -k2  [a2  -V0m2) 


=  -i 


.r*  f(a;a,b;v0m )ada  r°°  f(a;a,b;v0m)ada 


"  (cr  - vj )  +  •  4^e(a-  - rj ) 


-  -// ,  +  1 2  . 


(6.2) 


We  use  (5.7)  to  re-write  this 


=  r 

Jo 


[yo(«a)- J0(<2h)] 


J0(aa )  J0(ab) 


YM>ma)  yo(Yo  nP), 


\Ja2  -k2  (a2  -v0„,2) 


=  f 

Jo 


f(a;a,b;v0m)ada 

4a^e{a2-v J) 


=  .r*/(>/^2-.y2;a,^;Vom)<fe  ,  p  f(a;a,b\v0m)ada  =  -T  ±  T 
v0m2-k2+s2  h  r-^r,  ,  -  1  2 


4  a2  -k2  (a2  -v0m2) 


(6.3) 


The  first  integral  is  computed  in  Mathematica®  (reference  2).  For  the  second  integral,  we 
write 


W, 


k  f(a;a,b;v0m )ada  r»  f(a;a,b;v0m)ada 


4a2 -k2  (a2  -  v0m2 )  ^  4a 2  -k2  (a2  -  vj ) 


1 2 1  +  1 22  • 


(6.4) 


Using  (5. 1 1),  we  re-write  the  first  integral  in  the  form 


j  _  cl  f(4k2  +t2  ;a,b-,v0m)dt  |  p  f(a;a,b;v0m)ada 

2  Jo  f1  2  '-2i  Jv1 


I, 


f(4k2+t 2  ;a,b;v0m)dt 


/2  - 


(u,,,2  -£2' 


Vv0,„2-E+io-' 


f(4k2+t2;a,b;v0m)dt 

(u>m2  -^2) 


/2  - 
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+ 


f(a;a,b;v0m  )ada 
Va2  -  k 2  (a2  -v0„,2) 


(6.5) 


where  Z  is  a  sufficiently  large  integer  so  that  the  last  integral  is  free  of  singularities.  The  first  two 
of  these  integrals  are  computed  numerically  in  Mathematica®  (reference  2)  while  the  last  is 
evaluated  asymptotically.  We  note  that  the  singularity  at  v0m  is  removable  and  that  the  integrand 
has  a  finite  value  there. 


For  a  ->  oo 


1 


f(a;a,b-v0m)  =  — 


1 


1  sin(2aa)  sin(2  ba) 


na  \aY0(v0ma)  bY0(v0mb)  aY0(v0ma)  aY0(v0mb ) 


1  1 

-  +  ■ 


Yo(y0ma)  Y0  (V0mb) 


COS 


a(b-a )  +sin  a{b  +  a ) 


yfab 


O^a  ') 


and 


a 


4ar^e{ 


a'  -vn 


')  ,  L  (k' 


a  Ji¬ 


lt- 


k 

\aj 


1  + 


V  a  ) 


\a) 

1  + 


1- 


fy  V 

y0m 

V  a  ) 


1  (k42  f'’  ^ 


+ 


\ocJ 


Vn 


V  a  ) 


a 


cc 


so  that 


f{cr,a,b\vQm)a 
4a2  -k2  (a2  -vj) 


na 


1  1 

■  +  ■ 


aYo(v0ma)  bY0(v0mb ) 


+ 


O^a  4  j . 


We  can  then  write 


4k 


CO 

I 


f(a;a,b;v0m)ada 


+k< 


Y^¥(a--vJ)  2  n(C+k>) 


1  1 

+ 


aYo  (v0ma)  bY0(v0mb) 


From  (6.3),  (6.5),  and  (6.9),  the  final  expression  for  the  third  integral  is 


(6.6) 


(6.7) 


(6.8) 


(6.9) 
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*  fjylk2  -s2  ;a,b;v0m)ds  |  A^-a^-kt*  f(\lk2  +  t2  ;a,b-v0m)dt 
“  V0m  ~k2  +S2  +J°  {--(yj-k 2) 


r L  f(yjk2  +t2;a,b;v0m)dt 

1 

1  1 

J^0m2-i2+10^  j2  _ 

1  2  72) 

v-  ~k ) 

2  n 

L2+k2) 

_aY0(v0ma)  bY0(v0mb)_ 

(6.10) 


106 


NAWCADPAX/TR-20 13/115 


7.  THE  SECOND  INTEGRAL  IN  (3.14) 
From  (3.14),  the  fourth  integral  is 

J0  (a  a)  J0  ( ab ) 


1  = 


=  r 

Jo 


J0  (a  a)  J0  (ab) 

Yo(vona)  Yo  (vo„b) _ 

4a1  -k1  i^iXi)n  +4a2  -k2  -v0„,2) 


la(Ja2-k2  -k2) 

Yo(y^)  Y*<yJ>)  1  ’ 


a 


da . 


This  is  similar  to  the  previous  integral.  Letting 
J0(aa)  J0(ab ) 

g(a\(xM  v0l, )  =  — - -  -  — — —  • 

Yo  (vona)  IiW) 

and  following  Section  5,  we  write 


g(a;a,b;v0n)g(a;a,b;v0m)Uv0n2  -k2  4a2  -k2  -k2  )a 

I=  \  -  — - r - - — da 


=1, 


[ 


g(a\ a,b; v„n )g(a; a, 6; ) E/v0.2  -  k2  4k2  -a2  -k2) 
i'lk2  - a 2  Uv02  - k 2  +i^k2  -a2  ((ft"  j 
g(a-,a,b;va:l)g(a;a,b-,v0J^JvJ^-k2  sla2  -t!  -t! 


a 


4a2  -k2  y v0n 2  -/c2  +4a2  -  k 2 

Using  (5.5)  in  the  fonn 


a"  -  v, 


Om 


) 


da 


-  da  =  /,+/,. 


-  k2 

4k 2  -  a2  ^y]v0„2  -k2  +  i4k 2  -  a2  j 


^v0n2-k2a2  [f  v0„2 
4k2  -a2  (v0h2  -  a2)  (v0„2-«2) 


we  write 


r* 

W. 


r*  g(a;a,b;v0n)g(a;a 

V 

1  2 
u>„  ■ 

-k2  a2 

-  4-  i  - 

i  i 

u>„2 

Jo  /| 

(«2  -u>,„2) 

(N 

-a2 1 

(vo„2  -«2) 

(u,,2  -«2) 

ada 


=  _y  2  r*  g(a-,a,b\v0„)g(a\a,b;v0m) 


JO  /  2  2  \  /  2  2 

(u„  -a  jiUm  -a  ) 


(7.1) 


(7.2) 


(7.3) 


(7.4) 
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.  f~i  77 f*  g(a;a,b;v0n)g(a;a,b;v0m) 

^V°"  ~k  Jo  /,  2  2/,,  2  „2\/.  2  „ 


rrt«  =  ->A  W<V.’  -*W 


lk--a’-(vj  -er’)(v,.!  -a!) 


Using  (5.7),  we  re-write  this  as 


=  _  ir.g(a;a.<’;vD  )g(a;8.6;K  ) 

Jo  K!-«!)K2-«!) 


+  WV'o  n 


p*  g(V^r2  -t2 ; a, b; v0n )g(jk2  -t2 ; a, b; v0m )(k2  -t2^) 

°  (vo«2  _ + 12 )(v'0m2  -k2  +t~^j 


V On  A 1  WV0»  ^  A 2  * 


We  turn  to  the  second  integral  in  (7.3) 


As  with  (5.12),  for  the  first  of  these  integrals,  we  write 


(7.5) 


(7.6) 


(7.7) 
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g(W2  +  k2  ;a,b;v0n)g(^Jt2  +  k2  ;a,b;v0m)[Jv0n  - k2 1  - k 


Jv„„,  -k  +  i<r 


V0  n  - 


k2  +  ^)£2 -(vj  -k2) 


while  for  the  second 


722  -  £ 


g(a;a,b;v0n)g(a;a,b;v0m)ia2Jv0n2 -k2  -v0n2Ja2  - k 2  )a 


la2  -k2  (a2  -v02)(oc2  -v0m2) 


As  with  Sections  5  and  6,  we  make  the  upper  limit  of  integration  in  the  second  integral  in  (7.8) 
be  an  integer.  We  can  then  write 


Wt 


+k  +  k- ;a,b-,v0m)Uv0n  -k'l-k‘ 

g(sl?T¥ ; a, b; k0„ )g(V<;  +k = ; a, b; )(/>V  t - k2 J 

•JvnJ'-k1  +KT4  I  I  2  77  ,  .\r.2  /  2  ,  2  \"| 

[VvoH  +tj\t  -(v0m  -k  J 


lV0n  -< 


g(a;a,b;v0l,)g(a;a,b;v0m)[a2Jv02  -  k 2  -v0„V«2  -  k 2  )a 


I  2  77/  2  2  \  /  2  2 

/a  -k  - v0n  ){a  -v0m 


(7.10) 


The  last  integral  is  computed  asymptotically.  We  can  approximate  the  integrand  by 

g(a;a,b;  v0n)g(a;a,b; v0m ) ( a 2  Jv0„2  -  k 2  -  v02  V a 2  -  k 2  J a 
yja2  -  k 2  (a2  -V0  2  )(a2  -v0m2) 

A:-k2f  1  ,  1  1  ,  ^..-4, 


D  n  -k  1  ^  1  ^n,  -  4, 

— i - - 1 - +  0{oc  ) . 

na  aYQ(vQma)YQ(v0na)  bY0(v0mb)Y0(v0nb) 


(7.11) 


Then 


g(a;a,b;v0n)g(a;a,b;v0m)ia2Jv02  -  k 2  -  v0n2yja 2  - k 2  \a 


lex2 -k2  (a2  ~v0n2)(a2  -v0m2) 


Von 2  ~k2 


2 k (l2  +  k2 )  [aY0 ( v0ma)Y0 (v0na)  bY0 ( v0mb)Y0 (v0nb)  J  ' 


(7.12) 
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From  (7.6),  (7.8),  (7.9),  and  (7.12)  we  have  for  the  fourth,  and  last,  integral 

2  tk  g(a;a,b;v0n)g(a;a,b;v0m) 


>4 


-ada 


+  l\V, 


(^o„2  -«2)(^o,„2  -«2) 

k  g(yJk2  -t2 ;  a,  b;  v0n  )g(-Jk2  -t2 ;  a,  b;  v0m )  [k2  - 12 ) 


dt 


I, 


(vo«2  -k2  + 12  )(v0m2  -  k2  +  V  j 
+k 2 ;  a,  b;  v0„  )g(J?  +  k 2  ;a,b;  v0m  )^v0,2  -  k2 t-k2^ 


V0  n  ~k~  +t 


t2  - 


[VJ  -k1) 


dt 


JVl'o».2-*2  +io- 


2  ,  2 

R)„  ~k 


2/r(L2  +k 2) 


g(Vr  +  k2  ;a,b;v0n)g(yjt2  +k2;a, 

v0n  -k2t-  k2  j 

(v 

^o«2  +  * j 

1 

<N 

1 _ 1 

{y0m  ~k2 

); 

dt 


1 


1 


(n>ma)7o  (v0fla)  bY0  (v0mb)Y0  (v0nb) 


110 


(7.13) 


NAWCADPAX/TR-20 13/115 


APPENDIX:  EVALUATION  OF  INTEGRALS  FOR  SECTION  3 


We  first  observe  that,  from  (2.5)  and  (2.7) 

\‘pJ.Wp  =  pM\'.=  °- 


(A.l) 


We  next  evaluate  the  integral 


bdP0 m(p)  dJ 0  (ap) 


dp  dp  dp 


pdp  =  \ 

J  a 


b  d 


J0  (y0mp)~ 


dp  (u»ia) 
Yo  (voma) 


YP  (VomP) 


aJ{)  (ap) 
dp 


pd  p 


-  av, 


(vomP)~ 


dp  (V0ma) 

Y0  (YomO) 


Y\(ypmP) 


Jx  (ap)  pd  p 


=  |  £  di  (v0 mp)Jx  (ap)pdp 


^oMr‘ 
Yo(y0ma) 


£  Yi(yomp)JMp)pdp\  ■  (A.2) 


Using  standard  formulas  (reference  3,  p.  87),  we  proceed  to  evaluate  the  last  two  integrals 
pdP{)m(p)  dJp(ap)  av0m  , 

I  - , - ~ - Pdp  =  - j - 2  P  {  ^0,^2  (U»  P)^l  («P)  -  (Um  P)*/2  («P) 

Ja  dp  v0  -  a2 


Y0  (VQma) 

Yo(VOma) 


[ vo,„  Y2  (y0 m P)d , 1  (ap)  -aYl  (v0m  p)J 2  (ap)]  l  * 


gyp,, 

2  2 
Um  -a 


P  IE 


0  m 


Yo(v0ma) 


J i  (ap) 


a 


7oOom«) 


d 2  (ap) 


avQmp 


vn 


■a~  P 


a  v. 


0  m 


2  2 
Vpm  -a 


p 


a  v, 


0  m 


2  2 
~a 


P 


Y0  (VPma) 

dx  (y0mp)  -Jy^a\Yi(vpmp) 

Yo(V0  ma) 

d \  (y Pm P )  Yl  (V0, m P) 


~Vn 


dp(VpmP)~ 


dp(Vpma) 

Yo(yoma) 

b 


Yp(yP,nP) 


Jx  (ap)  -  J 0  (ap) 


>  J i  (ap) 


Yp(yPnP) 


ap 


dp(ap)  ba 


a  Vn 


2  2 

y0m  ~a 


Yp(yp,„b) 


J0  (ab)  -  a 


J°{d"a\ 

Yp(yp„P) 


J0  (a  a)  | 
(A.  3) 


or 
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tbdpo«,(p)  dJQ{ap) 

dp  dp 

Finally,  we  compute 

f bdPpm(p)  ^Jap) 
dp  dp 


pdp  = 


2a 


(  2  2 

~vom 


J0(aa )  J0(ab) 


ro(V«ma)  Y0(V0mb) 


(A.4) 


^Om  ^ On 


pd  p 
Jo  (V0ma) 


J^omP)-^r^Yx{y0mp) 

Yo(Yoma) 


J\  (V0nP)  ~  l°(l/°"a)  YMo nP) 


pd  p 


=  V0,nV0n  \  \h  Jx<P'QmP)Jx(VQnp)pd  p-Jo{V°mCl)  ^ 


Jq  (y On  J)  p 


v°.VK°"^f  J\(y0ap)Yl(y0np)pdP  +  — - r— - - 

Yo(VOna)la  Yo(VOma)  Yo(VOna) 


Y0^0,P) 

\'>Yi(y0mp)Jl(y0np)pdp 

J  a 

\\(V0mPWV0nP)Pd^ 


Yo(VOma) 

J0(V0n,a)  Jo(VOna)  f6 


(A.5) 


If  we  proceed  to  complete  this  calculation,  we  will  find  that  the  result  is  zero.  On  the  other  hand, 
if  m  =  n  ,  then 


JPon  (aP) 


dp 


pdp  =  v  J~\l 

J  a 


JSV0nP)~ 


JpQ'onJ) 

Yo(VOna) 


YX^0  nP) 


pd  p 


=  v, 


On 


[J i(v o„. P)\~  pdp -2 


Jq  (V0na) 


“1 2 


JoO/olla)  r h 

Yo(VOna) 

[[Y^o„P)f  Pdp 


J  Jl(VonP)Yl(VOnP)pdp 

J  a 


Jo  (V0na). 

Proceeding  with  the  same  kind  of  formulas  as  above,  we  find 

JPon^pV 


(A.6) 


dp 


=  '/0„  < 


p~ 


pd  P 


[Jx  (V0n  P)  )2  “  J 2  (Yon  P)J*  Oo„  P)  + 


Jo0/On(Y) 

Yo(VOna) 


-|2 

{( ■ YX  K„  P)  f  ~  Y2  O' On  P)Y0  0/O  n  P)) 


'  j"^(h'a\  ^[2JxO'onP)YxO'onP)-J20'onP)YoO'o„P)-JoO'onP)Y20'onP)^ 
YoO'onO)  2 


2  2 

Vo„  P 


IjOv«,p)-'yP^J<Oo,p) 

JxO'onP)- 

LL  Yoiy0„a)  \ 

L  YoO'onO)  \ 

nP) 
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+ 


^o(^o„a) 

Y0  (V0na). 
J0(V0n<*) 

Yo(vo„a) 


Jl  (Yon  P)  ~  YX  (Yon  P) 

Yo(vona) 


MVOnP)~ 


Mvo«a) 

Yo  (v0«a) 


Y0(^0nP) 


Yx(y*nP) 


2  2 

Von  P 


Jx  iy0np)~ 


Jq(vo,<*) 

Y0  iy0na) 


Yl(VonP) 


M  VPnP) 


Jp  (R)  »^0 
Yo(vo„a) 


MVpnP)- 


Jp(.Vpna) 

Yp  (y0na) 


Yx^PnP) 


Yx(Vp  nP) 


b 

a 


2  2 

Vq«  P 

2 


JA^PnP)- 


Jp  (yona) 
Yo  (Yo*a) 


Yx(Vo„P) 


|2 


b 

a 


2 

V0n  < 

b  2 

2 

a  2 

2  ' 

2 

JpiVpnb)  KVpnb_ 

Y0  (y0na)  xv0 na_ 

(A.7) 


0PqMp) 

dp 


pdp  = 


n 


1  _  1 

[r0(t-0„i»)]=  [r„(v0l,a)f  j 


(A.8) 
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PART  4 

NOTES  ON  COMPUTING  THE  INFINITE  SYSTEM  OF  EQUATIONS 


ABSTRACT 

This  is  the  fourth,  and  final,  part  of  the  report  on  the  formulation  of  the  problem  of  radiation  of 
a  coaxial  line  into  a  half  space  in  terms  of  BIEs.  In  it,  we  give  a  detailed  presentation  of  how  we 
compute  the  infinite  system  of  equations.  Specifically,  we  truncate  the  infinite  system  and 
consider  finite  systems  of  dimensions  from  lxl  to  11x11.  We  examine  the  convergence  of  the 
coefficients  of  these  systems  for  four  different  coaxial  lines  and  use  them  in  computing  far-field 
quantities,  such  as  directivity  and  gain. 
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1.  INTRODUCTION:  WHAT  IS  TO  BE  COMPUTED 

We  use  here  the  results  of  Part  3  to  compute  the  radiated  fields  of  a  coaxial  transmission  line 
that  opens  into  a  half  space.  We  use  Mathematical  (reference  1)  in  all  computations.  In  this 
section  we  assign  values  to  all  relevant  parameters. 

We  set  the  characteristic  impedance  of  the  coaxial  line  at  50  Q;  thus, 

50  =  Z  =  z'^±  =  (U) 

C  U  nt  4  /  nt  '  ' 

2  x  s0  2  n 


From  this 

Z  =  -  =  e100'r'/i“7^  =  2.3023 
a 


(1.2) 


or 


p  =  -  =  0.4343  . 
b 


For  all  modes  (except  the  TEM)  to  be  suppressed 


—  <1 
2 


(1.3) 


(1.4) 


or 


kb  <  1.3944  or  -<0.2219. 

A 


We  set  A  - 1  (300  MHz)  and  give  b  the  values  shown  in  Table  1.1. 


(1.5) 


Table  1.1.  Values  of  the  parameters  for  the  numerical  computation  ( A  -  1 ). 


Case 

b( m) 

a  (m) 

kb 

ka 

1 

0.05 

0.021715 

0.314159 

0.136439 

2 

0.10 

0.043430 

0.628319 

0.272879 

3 

0.15 

0.065145 

0.942478 

0.409318 

4 

0.20 

0.086860 

1.256637 

0.545757 
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2.  FIRST  INTEGRAL 

In  this  section,  we  compute  the  integral  in  (4.22)  of  Part  3.  The  imaginary  part  is  the  sum  of 
two  integrals.  Both  were  computed  in  Mathematical  (reference  1).  Their  sum  is  given  in  Table 
2.1  for  the  four  different  coaxial  lines  in  Table  1.1. 

Table  2.1.  Numerical  results  for  integral  In  ,  the  imaginary  part  of  (4.22). 


Case 

In 

1 

-  0.0000420305  i 

2 

-0.000649198  i 

3 

-  0.00309862  i 

4 

-0.00901561  i 

The  result  for  the  second  line  in  (4.22)  is  shown  in  Figure  2.2. 

Table  2.2.  Numerical  results  for  integral  In  ,  the  real  part  of  (4.22). 


Case 

In 

1 

0.0136816 

2 

0.0285517 

3 

0.0452372 

4 

0.0634187 

The  entire  integral  in  (4.22)  is  shown  in  Table  2.3. 

Table  2.3.  Numerical  results  for  integral  /, ,  the  entire  expression  in  (4.22). 


Case 

h 

1 

0.0136816-  0.0000420305  i 

2 

0.0285517-  0.000649198  i 

3 

0.0452372  -  0.00309862  i 

4 

0.0634187-  0.00901561  i 
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3.  SECOND  INTEGRAL 

The  second  integral  appears  in  (5.20)  of  Part  3.  We  compute  this  integral  for  the  four  cases  in 
Table  1.1  and  the  first  10  roots  of  the  transcendental  equation 

J,(y0na)Y0(y,„b)~  J0(v0nb)Y0(v0na)  =  0  .  (3.1) 

The  roots  for  the  four  cases  have  been  computed  in  Mathematica  (reference  1).  We  observe  that 
each  successive  case  (column)  is  equal  to  the  first  one  divided  by  the  number  of  the  column.  We 
note  that  Mathematica®  has  a  subroutine  that  finds  the  roots  of  this  equation  but  we  did  not  take 
advantage  of  it. 


Table  3.1.  First  ten  roots  of  the  transcendental  equation  (3.1) 


Root  No. 

Case  1 

Case  2 

Case  3 

Case  4 

1 

110.131 

55.0656 

36.7104 

27.5328 

2 

221.636 

110.818 

73.8786 

55.4089 

3 

332.867 

166.434 

110.956 

83.2168 

4 

444.02 

222.01 

148.007 

111.005 

5 

555.14 

277.57 

185.047 

138.785 

6 

666.243 

333.122 

222.081 

166.561 

7 

777.337 

388.668 

259.112 

194.334 

8 

888.424 

444.212 

296.141 

222.106 

9 

999.508 

499.754 

333.169 

249.877 

10 

1110.59 

555.294 

370.196 

277.647 

The  entire  expression  (5.20)  is  computed  in  Mathematica®.  The  results  are  displayed  in  Tables 
3.2.1  through  3.2.5. 

Table  3.2.1  Numerical  results  for  integral  I2,  the  expression  (5.20)  (roots  1  and  2) 


Case 

1 

Toot 

1 

2 

1 

1.26328  + 0.0036659  i 

-  2.07424  +  0.000668069  i 

2 

1.22771  +  0.0277552  i 

-  2.06343  +  0.00443626  i 

3 

1.18825  +  0.0853679  i 

-2.03808  + 0.0109216  i 

4 

1.16772  + 0.177285  i 

-  1.99193  + 0.0109216  i 
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Table  3.2.2  Numerical  results  for  integral  I2,  the  expression  (5.20)  (roots  3  and  4) 


Case 

Root 

3 

4 

1 

17.9066  + 0.0147716  i 

4.7176-0.000653382  i 

2 

17.7639  +  0.111543  i 

4.70708- 0.00433879  i 

3 

17.6066  + 0.341594  i 

4.68246-0.0100815  i 

4 

17.5259  + 0.705258  i 

4.63767-0.0107471  i 

Table  3.2.3  Numerical  results  for  integral  I2,  the  expression  (5.20)  (roots  5  and  6) 


Case 

Root 

5 

6 

1 

-3.52255  -0.00168172  i 

3.76781  -  0.00032758  i 

2 

-3.50631  -0.0126963  i 

3.76253  -  0.00217531  i 

3 

-  3.48842  -  0.0388686  i 

3.7502  -  0.00505523  i 

4 

-3.47925  -0.0802122  i 

3.72778  -  0.00539464  i 

Table  3.2.4  Numerical  results  for  integral  I2,  the  expression  (5.20)  (roots  7  and  8) 


Case 

] 

loot 

7 

8 

1 

-  11.05 12  -  0.00370696  i 

-  11.7554  +  0.00074151  i 

2 

-  11.0154-  0.0279843  i 

-  11.7434  + 0.00492407  i 

3 

-  10.9759-  0.0856635  i 

-  11.7155  + 0.0114437  i 

4 

-  10.9557-  0.17676  i 

-  11.6648  + 0.0122167  i 

Table  3.2.5  Numerical  results  for  integral  I2,  the  expression  (5.20)  (roots  9  and  10) 


Case 

Root 

9 

10 

1 

5.31419  +  0.00137324  i 

-  4.72606  +  0.000233462  i 

2 

5.30092  +  0.0103665  i 

-4.72226  +  0.00155033  i 

3 

5.28632  + 0.031732  i 

-4.71348  + 0.00360309  i 

4 

5.27883  +  0.0654731  i 

-4.69751  +  0.0038471  i 
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4.  THE  THIRD  INTEGRAL 


The  third  integral  is  given  by  (6.3)  in  Part  3.  The  entire  expression  is  computed  in 
Mathematica  .  The  results  are  displayed  in  Tables  4.1.1  through  4.1.5. 


Table  4.1.1  Numerical  results  for  integral  /3 ,  the  expression  (6.3)  (roots  1  and  2) 


Case 

Root 

1 

2 

1 

-  0.00354536  +  0.0000333409  i 

0.00433068  + 3.01548*  10'bi 

2 

-  0.00772607  +  0.000507352  i 

0.00875998  +  0.0000400965  i 

3 

-0.0141551  +  0.00236027  i 

0.0136065  + 0.000139937  i 

4 

-0.0219913  + 0.00661357  i 

0.0191409  + 0.000198389  i 

Table  4.1.2  Numerical  results  for  integral  /3 ,  the  expression  (6.3)  (roots  3  and  4) 


Case 

Root 

3 

4 

1 

-  0.0057484  +  0.0000443849  i 

-0.00269068  -  1.47166*10'b  i 

2 

-0.01 19776  + 0.000670672  i 

-0.00543566-  0.0000195511  i 

3 

-0.0199091  +  0.0030836  i 

-0.00797123  -0.0000681768  i 

4 

-0.026884  +  0.00849921  i 

-0.011534  -0.0000969721  i 

Table  4.1.3  Numerical  results  for  integral  /3 ,  the  expression  (6.3)  (roots  5  and  6) 


Case 

Root 

5 

6 

1 

0.000420177  - 3.02957*  10'6i 

-  0.00102198  -  4.91704*10'7  i 

2 

0.000916501  -0.0000457526  i 

-0.00206025  -6.53 125*  10'bi 

3 

0.0014815  -0.000210169  i 

-0.00314133  -  0.0000227721  i 

4 

0.00234308  -  0.000578554  i 

-0.00434858  -  0.0000324115  i 

Table  4.1.4  Numerical  results  for  integral  /3 ,  the  expression  (6.3)  (roots  7  and  8) 


Case 

Root 

7 

8 

1 

0.000739879  - 4.76895*  10'bi 

0.00185938  +  8.34656*10'vi 

2 

0.00154382  -0.0000720099  i 

0.00374666 +  0.00001 1086  i 

3 

0.00231508  -0.000330701  i 

0.00562887  +  0.0000386513  i 

4 

0.00340646  - 0.000910044  i 

0.00785099  +  0.0000550258  i 
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Table  4.1.5  Numerical  results  for  integral  I3 ,  the  expression  (6.3)  (roots  9  and  10) 


Case 

Root 

9 

10 

1 

-  0.00020 1 699  +  1 .37394*  1 0‘6  i 

0.000498893  +  2. 1 02 18*1 0"v  i 

2 

-  0.000431858  +  0.0000207448  i 

0.00100607  +  2.79208*10'6  i 

3 

-  0.000688448  +  0.0000952597  i 

0.0015313  +  9.73432*  10'6i 

4 

-0.00108081  +  0.000262104  i 

0.00210175  + 0.0000138596  i 
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5.  THE  FOURTH  INTEGRAL 


The  fourth  integral  is  given  by  (7.1)  in  Part  3,  and  in  computable  form  by  (7.13).  The  result  in 
(7.13)  depends  not  only  on  the  index  m  but  also  on  the  index  n.  Since  n  runs  from  1  to  10,  we 
have  10  times  as  much  data  as  in  the  previous  cases.  For  this  reason,  we  do  not  form  tables  but, 
rather  we  display  the  raw  data  from  Mathematical  For  the  first  10  roots,  we  have 

{{{-0.566603+0.002908  i, 0.00857833  +0.000529959  i,- 
5.38481  +  0.0117177  i , - 0 . 01 072 52  -  0 . 0 0 05183 08  i, 0.81978  -0.00133404 
i, -0 . 00604899-0 . 000259859  i, 2. 12689  -0.00294057  i , 0 . 0147767 
+  0.000588217  i , - 0 . 8 7 91 04  +  0 . 0 01 08 933  i,0. 00491084  +0.000185199 
i} , {0 . 00425818  +0.000263014  i , - 0 . 84 0763  +  0 . 0 0 0 04 7 93 8 1  i,0. 0279246 
+  0.00105982  i, 1.54291  -0.0000468842  i , - 0 . 0 03 6 74 76  -  0 . 0 0 012 06 58 
i, 1.03996  -0.0000235059  i , - 0 . 0 08 8 0162  -  0 . 0 0 02 6 5 962  i,- 
2.82511+0.000053208  i,0. 00344896  +0.0000985254  i,- 
1.01075+0.0000167524  i },{ -1 . 77903+0 . 00387126  i, 0.0185789 
+  0.00070551  i, -25 .442  +  0 . 0155992  i , - 0 . 022224 5  -  0 . 0 0 06 9  i,4.5937  - 
0.00177594  i, -0 . 0122724-0 . 000345938  i, 13. 1365  -0.00391463 
i,0. 0295536  +0.000783066  i , -5 . 7899+0 . 00145017  i, 0.00974883 
+  0.000246546  i } ,  { - 0 . 0 02 6 556 1 - 0 . 0 0 012 83 6 1  i, 0.769935  -0.0000233956 
i, -0 . 0166423-0 . 00051723  i , -1 . 64689+0 . 0000228813  i , 0 . 00215396 
+0.0000588857  i , -1 . 20977+0 . 0000114717  i, 0.00510935  +0.000129799 
i,3. 47599  -0.0000259675  i , - 0 . 0 01 98 93 9  -  0 . 0 0 0 04 8 084  i, 1.29426  - 
8 . 17578x10  6  i},{0. 162381  -0.000264239  i , - 0 . 0 014 64 57  -  0 . 0 0 0 04 8 1557 
i,2. 75418  -0.00106475  i, 0.00172414  +0.0000470971  i ,  - 
0.543394+0.000121219  i , 0 . 000942938  +0.0000236126  i,- 
1.64399  +  0.000267199  i , - 0 . 0 022 62 52  -  0 . 0 0 0 05344 95  i, 0.753789  - 
0.0000989835  i ,- 0 . 000741743  -  0 . 0000168284  i } , { - 0 . 0 0 0 9 984 6 5  - 
0.0000428871  i, 0.345853  -7 . 81679xl0"6  i  ,  -  0 . 0  06 13  222  -  0 . 0  0  0172  8 14 
i, -0.806248  +  7. 64495x10 ~6  i , 0 . 000787625  +0.0000196745  i,- 
0 . 625213+3 . 83287xl0"6  i,0. 00186647  +0.0000433678  i, 1.86557  - 
8 . 6761x10  6  i  , -0 . 000725125-0 . 0000160655  i, 0.71443  -2 . 73165xl0“6 
i},{0. 300869  -0.000415949  i , - 0 . 0 02 5 054 5  -  0 . 0 0 0 0758 03 7  i , 5 . 62488  - 
0.00167606  i, 0.00292745  +0.0000741373  i ,- 1 . 1741+0 . 000190816 
i,0. 00159898  +0.0000371694  i , -3 . 69107  +  0 . 000420608  i ,- 0 . 00382792 - 
0.0000841368  i, 1.74081  -0.000155814  i , - 0 . 0 012 5 978  -  0 . 0 0 0 02 64 9 02 
i} ,  {0 . 00182884  +0.0000727999  i , - 0 . 7 04 6 02  +  0 . 0 0 0 013 2 6 8 8  i, 0.0110647 
+  0.000293348  i, 1.73733  -0.0000129771  i , - 0 . 0 014 15 93  -  0 . 0 0 0 0333 971 
i, 1.39912  -6 . 50621x10  6  i , - 0 . 0 0334 1 1 9  -  0 . 0 0 0 073 6 15 9  i,~ 
4.29179+0.0000147275  i,0. 00129871  +0.0000272709  i,- 
1 . 67888+4 . 63691xl0"6  i  },{- 0 . 0967212  +  0 . 000119835  i  ,  0 . 000762736 
+  0.0000218391  i, -1 . 92827  +  0 . 000482874  i , - 0 . 0 0 08 8 56 03  -  0 . 0 0 0 02 13 5 9 
i,0. 418727  -0.0000549743  i , - 0 . 0 0 04 8 13 06  -  0 . 0 0 0 01 07 08 6  i, 1.35405  - 
0.000121178  i, 0.00115215  +0.0000242399  i , - 0 . 6 524 57  +  0 . 0 0 0 044 8 9 01 
i,  0 . 000377498  +7 . 63187xl0~6  i  },{  0 . 000486278  +0.0000183356  i,- 
0 . 201679  +  3 . 34192xl0"6  i,0. 00292129  +0.0000738832  i  ,  0 . 51754  - 
3 . 26845x10  6  i  ,  -  0 . 0  0  03  72  04  8  -  8 . 4 114  7x1  O'6  i, 0.428674  -  1 . 63  8  6  7x1 0"6 
i,  -  0 . 0008809  94  -0.000018  54 11  i,  -1 . 343  22  +  3 . 7  0  93  xlO"6  i,  0 . 0  0  0341725 
+  6 . 86851xl0’6  i, -0 . 53434  +  1 . 16786xl0"6  i  }},{{- 0 . 544 13  5  +  0 . 02 16  915 
i,0. 00803394  +0.00346806  i , -5 . 29309  +  0 . 0871752  i , - 0 . 01 02 082 - 
0.00339188  i, 0.809014  -0.00992268  i , - 0 . 0 057 922 5  -  0 . 0 017 0 056 
i,2. 10209  -0.0218709  i,0. 0141846  +0.00384942  i,- 
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0.869407+0.00810185  i,0. 00472671  +0.00121198  i },{ 0 . 00397402 
+  0.00171479  i, -0 . 835292  +  0 . 000274929  i,0. 0268598  +0.00689244 
i, 1.53667  -0.000268894  i 0 . 00355686  -  0 . 000784538  i, 1.0361  - 
0.000134814  i, -0 . 00853931-0 . 00172923  i , -2 . 81408+0 . 000305168 
i,0. 00335349  +0.000640577  i , -1 . 00638+0 . 0000960814  i},{- 
1.74106+0.0286746  i,0. 0178479  +0.00458514  i , -25 . 2715+0 . 11524  i,- 
0.0215319-0.00448442  i, 4. 57028  -0.0131171  i ,- 0 . 011922  -  0 . 00224833 
i , 13 . 072  -0.0289119  i,0. 0287416  +0.00508934  i , -5 . 75974+0 . 0107101 
i,0. 00948526  +0.00160237  i } , { - 0 . 0 02 5154  -  0 . 0 0 083 6 13 5  i,0. 766117  - 
0.000134058  i, -0 . 0161017-0 . 00336077  i , -1 . 6413 9+0 . 000131116 
i,0. 00209031  +0.000382542  i , - 1 . 2 0557+0 . 0 0 0 06 573 6 9  i,0. 00495753 
+  0.000843177  i, 3. 46201  -0.000148803  i , - 0 . 0 01 92 96 9  -  0 . 0 0 03 1234 7 
i, 1.28809  -0.0000468503  i  }  ,  { 0 . 159489  -0.00195615  i, -0.001414- 
0.000312798  i, 2. 73914  -0.00786156  i,0. 00167437  +0.000305926  i ,  - 
0.540922  +  0.000894839  i , 0 . 000917494  +0.00015338  i,- 
1.63608  +  0.00197235  i , - 0 . 0 02 1 96 93  -  0 . 0 0 034 71 94  i, 0.749676  - 
0.000730635  i , - 0 . 0 0 072 0 98 9  -  0 . 0 0 01 0 93 13  i },{- 0 . 000951622  - 
0.00027932  i, 0.344182  -0.0000447838  i ,- 0 . 00595562  -  0 . 0011227  i ,  - 
0.803284+0.0000438008  i , 0 . 000768246  +0.000127793  i,- 
0.622639+0.0000219602  i, 0.001819  +0.000281673  i, 1.85637  - 
0.0000497095  i ,- 0 . 0 0 07 0751 - 0 . 0 0 01 04343  i, 0.710204  -0.0000156509 
i},{0. 295914  -0.00307878  i , - 0 . 0 024  24 54  -  0 . 0 0 04 92  3 13  i,5.5944  - 
0.0123733  i, 0.0028485  +0.000481498  i , -1 . 16828+0 . 00140839 
i,0. 00155717  +0.000241406  i , -3 . 67066  +  0 . 00310427  i , - 0 . 0 03 723 1 - 
0.000546449  i, 1.72959  -0.00114995  i , - 0 . 0 01223 5  -  0 . 0 0 0172 04 8 
i} , {0 . 00174669  +0.000474114  i , - 0 . 7 0 0 98  +  0 . 0 0 0 076 0 155  i,0. 0107399 
+  0.00190566  i, 1.72976  -0.000074347  i , - 0 . 0 013 74 8  -  0 . 0 0 02 16 913 
i, 1.39203  -0.000037275  i , - 0 . 0 03234 8 5  -  0 . 0 0 04 78 1 07  i,- 
4.26555+0.0000843763  i, 0.0012527  +0.00017711  i,- 
1.66659  +  0.0000265657  i } , { - 0 . 0 9517 97  +  0 . 0 0 08 8 6 94 7  i , 0 . 000738688 
+  0.000141827  i, -1 . 91701  +  0 . 00356454  i ,- 0 . 000859986  -  0 . 000138712 
i,0. 416318  -0.000405732  i , - 0 . 0 0 04 6 6 8 94  -  0 . 0 0 0 06 954 51  i, 1.3451  - 
0.000894289  i, 0.00111042  +0.000157423  i , - 0 . 64 73 9  +  0 . 0 0 033 12 8 
i,0. 00036287  +0.0000495643  i },{ 0 . 000465729  +0.000119408  i,- 
0.200534+0.000019145  i, 0.0028385  +0.000479952  i, 0.514829  - 
0.0000187247  i  ,  -  0 . 0  0  03  624  02  -  0 . 0  0  0  054  63  0  9  i,  0.426016  -  9 . 38793xl0~6 
i, -0 . 000853457-0 . 000120414  i, -1.33318+0. 0000212507  i, 0 . 00 0330501 
+  0.0000446  062  i,  -  0 . 52  95  91  +  6 . 6  9  072  x  1  O'6  i  }  }  ,  {{-0.52  370  9  +  0.065  03  68 
i, 0.0022977  +0.00786479  i , -5 . 21129  +  0 . 26026  i ,- 0 . 00464348 - 
0.00769733  i, 0.799647  -0.0296142  i , - 0 . 0 03 0 06 96  -  0 . 0 03 8 5 971 
i,2. 08124  -0.0652676  i, 0.00788145  +0.00873736  i,- 
0.861588+0.0241768  i,0. 00274315  +0.00275099  i },{ 0 . 00113137 
+  0.00386437  i , - 0 . 82 6442  +  0 . 0 0 04 8 063 5  i,0. 0156485  +0.0154798 
i, 1.5279  -0.000470437  i ,- 0 . 00228503  -  0 . 00176153  i, 1.03156  - 
0.000235897  i , - 0 . 0 0574 0 97  -  0 . 0 03 8 823 8  i , -2 . 80337  +  0 . 000534009 
i,0. 002319  +0.00143815  i ,- 1 . 00283+0 . 000168135  i},{- 
1.70148+0.0849746  i,0. 0103677  +0.0102862  i , -25 . 1102+0 . 340059  i ,  - 
0.0142872-0.0100672  i, 4. 55118  -0.0386942  i , - 0 . 0 082 97 91 - 0 . 0 05 04 8 06 
i , 13 . 0274  -0.0852795  i, 0.0205462  +0.0114275  i ,- 5 . 742 04+0 . 03 158 97 
i,0. 00690648  +0.00359798  i },{- 0 . 00113422  -  0 . 00188271  i,0. 760584  - 
0.000234183  i , - 0 . 01 06 512  -  0 . 0 0754 175  i , -1 . 63568  +  0 . 000229214 
i,0. 00146939  +0.000858217  i , - 1 . 2 0244+0 . 0 0 01 14 93 7  i,0. 00358463 
+  0.00189149  i, 3. 45405  -0.000260189  i , - 0 . 0 014 1 951 - 0 . 0 0 07 0 06 64 
i, 1.28524  -0.0000819216  i},{0. 156387  -0.00579165  i ,-0.000903337- 
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0.000701136  i, 2. 72612  -0.0231776  i,0. 00117898  +0.000686209  i,- 
0.539298+0.0026373  i , 0 . 000669242  +0.000344089  i,- 
1.63204  +  0.00581244  i , - 0 . 0 016333 6  -  0 . 0 0 0778 92 6  i, 0.74795  - 
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6.  THE  SYSTEM  OF  EQUATIONS 


The  infinite  system  of  equations  is  given  by  (3.13)  and  (3.14)  in  Part  3.  We  can  solve  it  by 
truncating  it  to  a  finite  system  of  equations.  What  we  have  done  is  to  work  with  a  lxl  system 
(returning  TEM  mode  only)  up  to  a  1  lxl  1  system.  For  systems  of  order  higher  than  this,  we  need 
to  compute  the  integrals  above  for  more  roots  of  the  transcendental  equation  (3.1). 

Since  we  stop  with  a  11x11  system,  we  examine  the  convergence  of  the  first  four  coefficients 
as  a  function  of  the  order  of  the  system  for  each  of  the  four  coaxial  lines  in  Table  1.1.  In  Figures 
6.1  and  6.2,  we  display  the  real  and  imaginary  part,  respectively,  of  the  coefficient  A  as  a 
function  of  the  order  of  the  system.  We  see  that,  for  engineering  purposes,  A  has  stabilized  for 
the  11x11  system  and  for  all  four  coaxials.  If  we  put  these  results  under  the  microscope, 
however,  we  may  conclude  that  the  data  has  converged  to  at  most  two  significant  digits,  as  seen 
from  Figures  6.3  and  6.4. 

We  see  that,  as  the  radii  of  the  coax  get  larger,  the  reflection  coefficient  becomes  smaller  in 
magnitude.  This  means  that  more  energy  escapes  into  the  upper-half  space.  This  will  become 
evident  when  we  compute  the  gain  of  the  coaxial  line  as  an  antenna. 
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Figure  6. 1 .  The  real  part  of  the  reflection  coefficient  A  as  a  function  of  the  order  of  the 
system  of  equations.  The  smallest  system  is  lxl  and  the  largest  11x11. 
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Figure  6.2.  The  imaginary  part  of  the  reflection  coefficients  as  a  function  of  the  order 
of  the  system  of  equations.  The  smallest  system  is  lxl  and  the  largest  11x11. 


b  =  5  cm 


b  =  10  cm 


O □ □ 


b  =  1 5  cm 


<*»■  — — - <s>  — * -  ❖ - <> 


b  =  20  cm 


1 - ; - m 


127 


NAWCADPAX/TR-20 13/115 


In  Figures  6.5  through  6.8,  we  present  the  convergence  of  coefficient  B\.  Again,  we  have 
convergence  to  two  significant  digits. 
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Figure  6.5.  The  real  part  of  coefficient  B\  as  a  function  of  the  order  of  the  system  of  equations. 
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Figure  6.6.  The  imaginary  part  of  coefficient  B\  as  a  function  of  the  order  of  the  system  of 
equations.  The  top  curve  is  for  Case  1  (smallest  coaxial),  the  one  below  it  for  Case  2,  and  so  on. 

The  smallest  system  in  this  case  is  the  2x2. 
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Figure  6.7.  Detailed  version  of  Figure  6.5.  Case  1  is  the  top  left  graph.  Case  2  is  the  top  right, 
Case  3  the  bottom  left  and  Case  4  the  bottom  right. 
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Figure  6.8.  Detailed  version  of  Figure  6.6.  Case  1  is  the  top  left  graph.  Case  2  is  the  top  right, 
Case  3  the  bottom  left  and  Case  4  the  bottom  right. 
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In  Figures  6.9  through  6.12,  we  examine  the  convergence  of  coefficient  /A  The  results  are 
displayed  in  Figures  6.9-6.12  and,  again,  we  get  convergence  to  two  significant  digits. 
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Figure  6.9.  The  real  part  of  coefficient  B2  as  a  function  of  the  order  of  the  system  of  equations. 
The  smallest  system  is  3x3  and  the  largest  is  1  lxll. 
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Figure  6.10.  The  imaginary  part  of  coefficient  B2  as  a  function  of  the  order  of  the  system  of 
equations.  The  smallest  system  is  3x3  and  the  largest  11x11. 
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Figure  6.1 1.  Detailed  version  of  Figure  6.9.  Case  1  is  the  top  left  graph.  Case  2  is  the  top  right, 
Case  3  the  bottom  left  and  Case  4  the  bottom  right. 
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Figure  6.12.  Detailed  version  of  Figure  6.10.  Case  1  is  the  top  left  graph.  Case  2  is  the  top  right. 
Case  3  the  bottom  left  and  Case  4  the  bottom  right. 
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In  Figures  6.13  through  6.16,  we  examine  the  convergence  of  coefficient  B 3.  Convergence 
appears  to  be  to  two  significant  digits. 
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Figure  6.13.  The  real  part  of  coefficient  B  3  as  a  function  of  the  order  of  the  system  of  equations. 

The  smallest  system  is  4x4  and  the  largest  1  lxl  1. 
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Figure  6.14.  The  imaginary  part  of  coefficient  f?3  as  a  function  of  the  order  of  the  system  of 
equations.  The  top  curve  is  for  Case  1  (smallest  coaxial),  the  one  below  it  for  Case  2,  and  so  on. 

The  smallest  system  in  this  case  is  the  4x4. 
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Figure  6.15.  Detailed  version  of  Figure  6.13.  Case  1  is  the  top  left  graph.  Case  2  is  the  top  right, 
Case  3  the  bottom  left  and  Case  4  the  bottom  right. 
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Figure  6.16.  Detailed  version  of  Figure  6.14.  Case  1  is  the  top  left  graph.  Case  2  is  the  top  right, 
Case  3  the  bottom  left  and  Case  4  the  bottom  right. 


137 


NAWCADPAX/TR-20 1 3/ 1 15 


In  Figure  6.17,  we  have  plotted  the  magnitude  of  the  coefficients  B\  to  B io,  normalized  to  the 
magnitude  of  the  reflection  coefficient  A.  The  numbers  for  the  coefficients  came  out  of  the  11x11 
system.  As  we  saw  above,  the  accuracy  of  these  coefficients  does  not  go  beyond  the  first  couple 
of  digits.  This  graph,  however,  is  useful  in  displaying  their  order  of  magnitude  relative  to  that  of 
A. 


Figure  6.17.  Magnitude  of  coefficients  B\  to  B  i0  relative  to  that  of  reflection  coefficient  A. 

The  numbers  are  those  of  the  11x11  system. 


In  Figure  6.18  we  plot  the  real  part  of  A  on  the  horizontal  axis  and  the  imaginary  on  the 
vertical.  The  values  of  the  four  coax  cases  are  plotted  and  connected  in  sequence  using  straight- 
line  segments.  It  is  clear  that,  as  the  dimensions  of  the  coax  increase,  the  real  part  of  A  decreases 
while  the  absolute  value  of  the  imaginary  part  increases. 
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Real  Part  of  Coefficient  A 


Figure  6.18.  The  real  and  imaginary  parts  of  the  reflection  coefficients  for  the  four  coax  sizes. 
The  four  computed  points  are  connected  using  straight-line  segments. 

From  Figures  6.19  and  6.21,  for  B\  and  B2,  both  the  real  part  and  the  absolute  value  of  the 
imaginary  part  increase. 

From  Figure  6.20,  both  the  absolute  value  of  the  real  part  and  the  imaginary  part  of  B2 
increase. 
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Figure  6.19.  The  real  and  imaginary  parts  of  the  reflection  coefficient  B\  for  the  four  coax  sizes. 
The  four  computed  points  are  connected  using  straight-line  segments. 


Real  Part  of 

Figure  6.20.  The  real  and  imaginary  parts  of  the  reflection  coefficient  Bj  for  the  four  coax  sizes. 
The  four  computed  points  are  connected  using  straight-line  segments. 
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Real  Part  of  B ^ 

Figure  6.21.  The  real  and  imaginary  parts  of  the  reflection  coefficient  Bj  for  the  four  coax  sizes. 
The  four  computed  points  are  connected  using  straight-line  segments. 
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7.  FAR  FIELD 


In  this  section  we  compute  quantities  related  to  the  far  field.  The  far  field  is  given  by  (7.22) 
and  (7.23)  of  Part  2.  For  the  calculations,  we  use  the  coefficients  A,  B\,  B2,  and  B2.  We  compute 
far-field  amplitude  and  phase,  far-field  intensity  (average  radiated  power  density),  total  radiated 
average  power,  directivity,  and  gain.  The  definitions  we  used  are  the  following. 

Far-field  amplitude  and  phase: 


FarFieldAmplitude  =  F(3) 


(7.1) 


FarFieldPhase  =  tan 


Im{F(i9)} 

Re{F(,7)} 


(7.2) 


where,  from  (7.22)  of  Part  2, 

1  +  A  J0  ( ka  sin  3)  -  JQ  (kb  sin  9) 


F{9)  = 


Hz) 

2k2  sin.9  3 


sin  .9 


n 


I*. 


\K,2  -k2 

.7,,  (ka  sin  $) 

J0(A:hsin  <9) 

2  -k2  sin2 19 

7o(R)„6) 

Far-field  intensity  (reference  2,  p.  38): 


FarFieldlntensity  = 


F(»9) 

2Z„ 


(7.3) 


(7.4) 


Total  radiated  average  power  (reference  2,  p.  38): 


fin  fitl  2  _z '  i  i  rji/i  2 

TotalRadiatedAveragePower  =  dcp^  ^  d i9sini9  =  kY0  d9sin9  F(9)~  . 


F(3Y 


*n!  2 


2Zn 


(7.5) 


Directivity  (reference  2,  p.  39): 


Directivity 


FarFieldlntensity 
T  otalRadiatedAveragePower 
2k 


^  FarFieldlntensity 

T  otalRadiatedAveragePower 


(7.6) 


Average  input  power: 
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AveragelnputPower  = - = - 

2  Zc  100 


(7.7) 


since  the  characteristic  impedance  of  the  line  is  50  Q  and  the  input  is  one  Volt. 
Gain  (reference  2,  p.  58): 


Gain 


FarFieldlntensity 
AveragelnputPower 
2  n 


2 k  x  FarFieldlntensityx  100  . 


(7.8) 


With  these  definitions,  we  proceed  to  display  some  of  the  calculations  we  performed.  We 
begin  with  the  far- Held  amplitude.  The  results  are  shown  in  Figure  7.1.  The  polar  angle  6  is 
measured  from  the  perpendicular  to  the  infinite  plane  toward  the  plane  (i.e.,  it  is  equal  to  90  deg 
minus  elevation).  We  see  that,  as  the  dimensions  of  the  line  increase,  so  is  the  energy  that 
escapes  into  the  upper-half  space.  Equivalently,  as  the  operating  frequency  approaches  the  cut¬ 
off  frequency,  more  and  more  energy  escapes  into  the  upper  half  space. 


0  20  40  60  80 

Polar  Angle  6  (deg) 

Figure  7.1.  Far-field  amplitude  for  the  four  coaxial  cases. 
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Figure  7.2  shows  the  phase  of  the  far  field.  The  four  different  cases  are  shown  separately 
because  of  the  slow  variation  of  phase.  As  the  size  of  the  coax  increases,  so  does  the  variation. 


174.51675 

174.5167 

174.51665 

174.5166 

174.51655 

174.5165 

174.51645 

174.5164 

174.51635 


Polar  Angle  6  (deg) 


Polar  Angle  0  (deg) 


40  60 

Polar  Angle  0  (deg) 


Polar  Angle  6  (deg) 


Figure  7.2.  Far-field  phase  for  the  four  coaxial  cases. 
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The  results  for  directivity  are  shown  in  Figure  7.3.  The  variation  in  the  four  cases  is  small  but 
there  is  a  clear  tendency  for  the  directivity  to  rise  as  the  coaxial  gets  bigger. 


Polar  Angle  0  (deg) 

Figure  7.3.  Directivity  (dB)  of  four  coaxials.  Case  1  is  the  bottom  figure  while  Case  4  is  the 
top  one.  The  rest  follow  in  a  clockwise  fashion. 
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The  gain  results  are  shown  in  Figure  7.4.  Clearly,  the  gain  improves  toward  the  horizon  and 
with  increasing  coaxial  size. 


Polar  Angle  6  (deg) 

Figure  7.4.  Variation  of  gain  with  polar  angle  for  the  four  coaxials. 


We  have  performed  additional  calculations  that  appear  in  the  paper  that  came  out  of  this  study 
(reference  3).  There,  we  compared  the  coefficient  A  to  that  obtained  by  Bird  (reference  4)  and 
found  it  to  be  in  agreement  to  two  significant  digits  in  all  four  cases.  We  also  computed  gain 
relative  error  when  using  only  the  lxl  system  rather  than  the  11x11,  and  we  did  the  same  for  the 
normalized  admittance  at  the  opening  to  the  half-space. 
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CONCLUSIONS 

We  have  solved  the  classical  problem  of  radiation  by  a  coaxial  line  using  BIEs.  Because  of  the 
circular  symmetry  in  the  geometry,  we  can  convert  the  vector  integral  equations  to  three  scalar 
equations  and  we  show  that,  the  solution  of  any  of  these  will  determine  the  unknown  current 
density.  We  express  the  latter  as  an  infinite  series  in  the  modes  of  the  coaxial  line  and  we 
determine  the  unknown  coefficients  by  taking  advantage  of  the  orthogonality  of  the  modes.  This 
results  in  an  infinite  system  of  linear  equations  with  the  unknowns  being  the  coefficients  of  the 
infinite  series.  In  the  last  part  of  the  report,  we  have  presented  numerical  examples.  The  principal 
asset  of  this  method  is  that  we  can  provide  an  engineering  solution  to  the  problem  by  truncating 
the  infinite  system  to  a  small  system  of  linear  equations.  In  Part  4,  the  largest  system  we  used  has 
11  unknowns  (11x11).  By  contrast,  use  of  other,  well  known  numerical  methods,  such  as  finite 
elements  or  the  method  of  moments  or  a  hybrid  of  the  two,  would  require  solution  of  systems 
with  hundreds  if  not  thousands  of  unknowns.  In  closing,  we  wish  to  mention,  that  the  precise 
same  method  can  be  applied  to  radiation  from  a  circular  or  rectangular  waveguide  or  any  other 
guide  whose  natural  modes  are  known  analytically. 
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